Skip to contents

Reconciles per-condition (univariate) SuSiE fine-mapping into a single set of merged credible sets. Each row of the supplied QtlFineMappingResult is treated as one condition (its topLoci carrying that condition's credible sets); credible sets that share variants across conditions are unioned via connected components, and every variant is reported with its merged credible-set label plus the maximum and median PIP across the conditions it appears in. A typical use is selecting a representative lead variant per merged credible set to assemble the "strong" input for mashPipeline.

Usage

mergeSusieCs(fineMappingResult, coverage = 0.95)

Arguments

fineMappingResult

A QtlFineMappingResult (or any FineMappingResult) produced by per-condition SuSiE fine-mapping. Each entry's topLoci must carry a credible-set column (cs_<coverage*100>, e.g. cs_95, with values such as "susie_1" where the trailing integer is the set index and _0 means "not in a credible set") and a PIP column.

coverage

Credible-set coverage level selecting the cs_* column (default 0.95 -> cs_95).

Value

A data.frame with one row per variant: variant_id, credibleSetNames (the merged credible-set label), maxPip and medianPip; or NULL when no credible sets are present.