pull.markers: Extract markers with allele frequencies in specified range

In kbroman/fingers: Identifying Clusters of Related Individuals

Extract markers with allele frequencies in specified range

Description

Extract markers from a RAPD data set that have allele frequencies within a specified range.

Usage

pull.markers(dat,lo=0.1,hi=0.6,f=freq(dat))

Arguments

dat	A matrix of size (n.ind x n.mar) containing RAPD phenotypes, with 1 indicating the presence of a band and 0 indicating absence.
lo	Lower bound for band allele frequency.
hi	Upper bound for band allele frequency.
f	Vector of band allele frequencies (included in order to avoid recalculating it, if possible).

Value

A matrix, like the argument dat, but containing only those markers with band allele frequency between lo and hi.

Author(s)

Karl W Broman [e-mail]

See Also

freq

Examples

data(shiff1)
f <- freq(shiff1)
subset <- pull.markers(shiff1, 0.1, 0.6, f)

kbroman/fingers documentation built on May 17, 2023, 11:50 p.m.