findpeaks: trim a dataset to only peak(s)

In coloc: Colocalisation Tests of Two Genetic Traits

trim a dataset to only peak(s)

Description

tries to be smart about detecting the interesting subregion to finemap/coloc.

Usage

findpeaks(d, maxz = 4, maxr2 = 0.1, do.plot = FALSE)

Arguments

d	a coloc dataset
maxz	keep all snps between the leftmost and rightmost snp with |z| > maxz
maxr2	expand window to keep all snps between snps with r2 > maxr2 with the left/rightmost snps defined by the maxz threshold
do.plot	if TRUE, plot dataset + boundaries

Details

Differs from findends by finding multiple separate regions if there are multiple peaks

Value

logical vector of length d$position indicating which snps to keep

Author(s)

Chris Wallace

See Also

findends

coloc documentation built on Oct. 3, 2023, 5:07 p.m.