Wrapper Function for SuSiE RSS with Dynamic L Adjustment

Performs SuSiE RSS analysis with dynamic L adjustment. Supports both z+R (correlation matrix) and z+X (genotype matrix) interfaces.

Usage

susie_rss_wrapper(
  z,
  R = NULL,
  X = NULL,
  n = NULL,
  L = 10,
  max_L = 30,
  l_step = 5,
  coverage = 0.95,
  stochastic_ld_sample = NULL,
  ...
)

Arguments

z

Z score vector.

R

LD correlation matrix. Mutually exclusive with X.

X

Genotype matrix (samples x variants). Mutually exclusive with R. When provided, susie_rss uses the low-rank X interface.

n

Sample size.

L

Initial number of causal configurations.

max_L

Maximum number of causal configurations.

l_step

Step size for increasing L when the limit is reached.

stochastic_ld_sample

Stochastic LD parameter passed to susie_rss. NULL (default): no variance inflation. TRUE: infer sketch size from X (requires X, not R). Integer: explicit sketch size (for R only).

...

Extra parameters passed to susie_rss (e.g., var_y, coverage).

Value

SuSiE RSS fit object after dynamic L adjustment