Skip to contents

QTL Enrichment Pipeline (Genome-Wide)

Source: R/qtlEnrichmentPipeline.R
qtlEnrichmentPipeline.Rd

Genome-wide pipeline that computes per-pair (GWAS study, QTL context) enrichment estimates by passing the GWAS PIP vector and the QTL credible-set posteriors to qtlEnrichment. The returned table feeds colocPipeline via its enrichment argument.

Not gene-parallelisable: the enrichment estimator runs over the full genome of GWAS PIPs and the full collection of QTL fits at once.

Usage

qtlEnrichmentPipeline(
  gwasFineMappingResult,
  qtlFineMappingResult,
  numGwas = NULL,
  piQtl = NULL,
  lambda = 1,
  impN = 25,
  numThreads = 1L,
  ...
)

Arguments

gwasFineMappingResult

See above.

qtlFineMappingResult

See above.

numGwas

Number of GWAS variants used to estimate piGwas. When NULL (default) it is estimated from the data — bias warning applies if the input PIP vector is not genome-wide.

piQtl

Per-variant prior of being a QTL causal variant. NULL (default) estimates from the data.

lambda

Shrinkage parameter for the enrichment estimator. Default 1.0.

impN

Number of imputed samples used by the estimator. Default 25.

numThreads

Number of threads used by qtlEnrichment. Default 1.

...

Additional arguments forwarded to qtlEnrichment.

Value

A data frame with one row per (gwasStudy, qtlStudy, qtlContext) triple and columns gwasStudy, qtlStudy, qtlContext, enrichment, enrichmentSe, enrichmentLogOdds, plus any extras the underlying estimator emits. Suitable as the enrichment argument to colocPipeline (which joins on the same triple).

Inputs

  • gwasFineMappingResult: a genome-wide GwasFineMappingResult (one row per (study, LD block) tuple). Each entry's FineMappingEntry$trimmedFit must carry a pip vector.

  • qtlFineMappingResult: the genome-wide QtlFineMappingResult. Each entry's trimmedFit must carry alpha, pip, and prior-variance fields (V).

LD-sketch identity check

The GWAS FineMappingResultBase must have a non-NULL ldSketch (RSS-derived). If the QTL FMR also has a non-NULL ldSketch, the two must match exactly. When the QTL FMR's ldSketch is NULL (individual-level QTL fit), validation is skipped on the QTL side.