get.psubset: deprecated function

In RVS: Computes estimates of the probability of related individuals sharing a rare variant

Description

This function is deprecated with version >= 2.0 and should not be used, instead use multipleFamilyPValue

Usage

get.psubset(vec, not, pshare.data)

Arguments

vec	a vector of names of all families where a variant is seen
not	a vector of names of families where not all affected subjects share the rare variant
pshare.data	a data frame with at least two of the following columns: pshare: vector of RV sharing probabilities ped.tocompute.vec: vector of names of the families whose sharing probability is contained in pshare. The names in the arguments vec and not must be found in ped.tocompute.vec

Value

P-value of the exact rare variant sharing test requiring sharing by all affected subjects.

Examples

data(samplePedigrees)
notSharedFams <- c(15159, 15053, 15157)
famids <- sapply(samplePedigrees, function(p) p$famid[1])
notShared <- famids %in% notSharedFams
probs <- sapply(samplePedigrees, RVsharing)
get.psubset(famids, notShared, data.frame(pshare=probs,
ped.tocompute.vec=famids))

RVS documentation built on Nov. 8, 2020, 6:57 p.m.