Skip to contents

Reads the regression outputs produced by `polyfun/ldsc.py` for a single polyfun run (one trait, one annotation set) and returns them as a tidy list ready for downstream standardization. Hides the underlying file formats; downstream code consumes only modeling quantities.

Usage

readSldscTrait(prefix)

Arguments

prefix

Character. Path prefix to the polyfun outputs for one trait/run. The function appends `.results`, `.log`, and `.part_delete` to this prefix. Example: `"/path/to/cwd/CAD_META.filtered.sumstats.gz"`.

Value

A named list. See `sldscPostprocessingPipeline` for components.

Examples

prefix <- file.path(
  system.file("extdata", "sldsc_trait", package = "pecotmr"),
  "sumstats.parquet")
run <- readSldscTrait(prefix)
names(run)
#>  [1] "categories"   "tau"          "tauSe"        "enrichment"   "enrichmentSe"
#>  [6] "enrichmentP"  "propH2"       "propSnps"     "h2g"          "tauBlocks"   
#> [11] "nBlocks"