Probability of Detection

Description

Find the probability of detection.

Usage

1
probDetEqual3(MAF, Xmean, T, N.p, error)

Arguments

MAF
Xmean
T

The threshold, usually calculated from diffVariantError()

N.p
error

Examples

 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
  }