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"                     
#>  [3] "tauSe"                    "enrichment"              
#>  [5] "enrichmentSe"             "enrichmentP"             
#>  [7] "propH2"                   "propSnps"                
#>  [9] "h2g"                      "annotationJackknifeCoefs"
#> [11] "nBlocks"