Standardize tau and compute EnrichStat for one polyfun run — standardize_sldsc

Applies the Gazal standardization $\tau^*_C = \tau_C \cdot sd_C \cdot M_{ref} / h^2_g$ to the point and to each jackknife block. For `mode = "single"`, additionally computes EnrichStat and back-solves its standard error from polyfun's reported `Enrichment_p` using $|Z| = \Phi^{-1}(1 - p/2)$.

Usage

standardize_sldsc_trait(
  trait_data,
  sd_annot,
  M_ref,
  target_categories = NULL,
  mode = c("single", "joint")
)

Arguments

trait_data: List from read_sldsc_trait.
sd_annot: Named numeric vector from compute_sldsc_annot_sd.
M_ref: Scalar from compute_sldsc_M_ref.
target_categories: Character vector or NULL. If NULL, intersects `trait_data$categories` with `names(sd_annot)`.
mode: Character: `"single"` or `"joint"`.

Value

A list with `summary` (data frame), `tau_star_blocks` (matrix), `h2g`, `n_blocks`, `mode`.