View source: R/pvalue_functions.R
get.psubset | R Documentation |
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.
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 (RV). |
pshare.data |
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 <alexandre.bureau@msp.ulaval.ca>
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.
RVsharing,RVgene
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.