Screen cTWAS regions

Step 3 of the three-step ctwasPipeline: runs ctwas::screen_regions on the estCtwasParam result and returns the screened-region set. Use this entry point to substitute hand-tuned priors for the ones estimated in step 2 (e.g. when the accurate EM diverges to NaN and you want to recover the prefit values).

Usage

screenCtwasRegions(estResult, L = 5L, ncore = 1L, ...)

Arguments

estResult

A list returned by estCtwasParam.

L

Unused. Retained for call-site compatibility with ctwasPipeline; ctwas's screening always uses the single-effect (SER) model and ignores L. L is applied by finemapCtwasRegions downstream.

ncore

Number of cores.

...

Additional arguments forwarded to ctwas::screen_regions (e.g. min_nonSNP_PIP, min_snp_pval, min_var, min_gene).

Value

The estResult list augmented with screen_res (the full ctwas output) and screened_region_data.