Probability of sharing of rare variants in a subset of families
Description
Computing probability of sharing of rare variants in a subset of families where rare variants are seen based on precomputed familyspecific rare variant sharing probabilities.
Usage
1  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 (RV). 
pshare.data 
a data frame with at least the two following columns:

Details
All the subsets of families of size equal or inferior to the lenght of not
are created, and the joint probability of each such subset not sharing a rare variant and the remaining families sharing a rare variant is obtained as the product of the familyspecific rare variant sharing probabilities or its complement. The function then sums the pattern probabilities inferior or equal to the probability of the observed pattern of the not
families not sharing a rare variant and the remaining families sharing a rare variant.
Value
Pvalue of the exact rare variant sharing test requiring sharing by all affected subjects.
Author(s)
Alexandre Bureau <alexandre.bureau@msp.ulaval.ca>
References
Bureau, A., Younkin, S., Parker, M.M., BaileyWilson, J.E., Marazita, M.L., Murray, J.C., Mangold, E., AlbachaHejazi, H., Beaty, T.H. and Ruczinski, I. (2014) Inferring rare disease risk variants based on exact probabilities of sharing by multiple affected relatives. Bioinformatics, 30(15): 218996.
See Also
RVsharing,RVgene
Examples
1 2 3 4 5 6 7 8  data(ped.list)
# Setting up pshare.data data frame
pshare.data=data.frame(pshare=sapply(ped.list[1:39],function(ped)
RVsharing(ped)@pshare),ped.tocompute.vec=sapply(ped.list[1:39],function(ped) ped$famid[1]))
vec = c("15008","15053","15157")
get.psubset(vec,not="15008",pshare.data)
get.psubset(vec,not=c("15008","15053"),pshare.data)
