Probability of sharing of rare variants in a subset of families
Computing probability of sharing of rare variants in a subset of families where rare variants are seen based on precomputed family-specific rare variant sharing probabilities.
a vector of names of all families where a variant is seen.
a vector of names of families where not all affected subjects share the rare variant (RV).
a data frame with at least the two following columns:
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 family-specific 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.
P-value of the exact rare variant sharing test requiring sharing by all affected subjects.
Alexandre Bureau <firstname.lastname@example.org>
Bureau, A., Younkin, S., Parker, M.M., Bailey-Wilson, J.E., Marazita, M.L., Murray, J.C., Mangold, E., Albacha-Hejazi, 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): 2189-96.
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)) vec = c("15008","15053","15157") get.psubset(vec,not="15008",pshare.data) get.psubset(vec,not=c("15008","15053"),pshare.data)