get.psubset {RVS} | R Documentation |
This function is deprecated with version >= 2.0 and should not be used, instead use multipleFamilyPValue
get.psubset(vec, not, pshare.data)
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 |
P-value of the exact rare variant sharing test requiring sharing by all affected subjects.
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))