Description Usage Arguments Examples
Find the probability of detection.
1 | probDetEqual3(MAF, Xmean, T, N.p, error)
|
MAF |
|
Xmean |
|
T |
The threshold, usually calculated from diffVariantError() |
N.p |
|
error |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function(MAF, Xmean, T, N.p, error)
{
pDetEqual <- 0;
for(n in 0:(2*N.p))
{
pConf <- dbinom(n,2*N.p,MAF);
pDet.not <- 0;
theta <- (1-error)*n/(2*N.p) + error*(1-n/(2*N.p))
pDet <- pbinom(T-1, Xmean, theta, lower.tail=FALSE)
pDetEqual <- pDetEqual + pConf*pDet;
}
pDetEqual
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.