Description Usage Arguments Details Value Author(s) Examples
The assay sensitivity is the minimum number of copies that the digital PCR assay can detect. Users provide serial dilution results in the format of counts of positive and total reaction wells. The output is the estimated assay sensitivity and the copy number per well in the initial dilute.
1 | digitalPCR(pos, neg, dilution, Nboot, single.copy = c("FALSE", "TRUE"), upper.copy=100)
|
pos |
a numeric vector of number of positive reactions in the order of dilutions |
neg |
a numeric vector of number of negative reactions in the order of dilutions |
dilution |
a numeric vector of folds of dilutions relative to the initial dilution |
Nboot |
number of bootstrapping |
single.copy |
character tells whether the assay is assumed single-copy sensitive. default "FALSE" |
upper.copy |
the upper bound of copy number in highest concentration |
the pos, neg and dilution vector must be in the same order from high to low concentrations.
The returned value is a list of four components:
mean copy number |
mean copy number |
sd copy number |
standard deviation of copy number |
copy numbers |
bootstrapped copy number estimates |
thresholds |
bootstrapped assay sensitivity estimates |
Xutao Deng
1 2 3 4 5 6 7 8 9 10 11 12 13 | #dilution at 1,2,4,8 fold
pos=c(221,97,39,14)
dilution=c(1,2,4,8)
neg=384-pos
#note in practice, set this 100 or larger
Nboot=10
result1=digitalPCR(pos, neg, dilution, Nboot)
print (paste(result1$"mean copy number", result1$"sd copy number"))
hist(result1$"thresholds",xlim=c(-2,10))
#the following example contains only 1 dilution
result2=digitalPCR(230, 355, 1, Nboot, "TRUE")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.