Skip to contents

Combine QtlSumStats collections

Source: R/qtlSumStats.R
combineQtlSumStats.Rd

Row-bind two or more QtlSumStats collections into one – e.g. per-region or per-context association results written by separate pipeline steps. Per-element metadata columns are carried through; a collection lacking an optional column is NA-padded.

Usage

combineQtlSumStats(..., ldSketch = NULL)

Arguments

...

Two or more QtlSumStats objects, or a single list of them.

ldSketch

Optional genotype panel (see readGenotypes) to attach to the combined collection, overriding the unioned one. Default NULL.

Value

A single combined QtlSumStats.

Details

The three collection-level slots are merged rather than taken from the first input, on the same rules as combineGwasSumStats: genome and the summaryStatsQc options must agree, the per-element qcInfo$entryAudit concatenates in element order, and the LD sketches union over the shared genotype handle (all-NULL stays NULL, which is the usual case for individual-level QTL data).

See also

combineGwasSumStats, combineFineMappingResults, combineTwasWeights

Examples

data(qtlSumStatsExample)
combineQtlSumStats(qtlSumStatsExample)
#> QtlSumStats: 1 entries, genome build hg19
#>   1 studies, 1 contexts, 1 traits
#>   LD sketch: plink1 @ pecotmr://extdata/toy_canonical