findpeaks {coloc}

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