estimatePrimerDimer {MassArray} | R Documentation |
Function to estimate level of signal due to primer dimers in a given spectrum
estimatePrimerDimer(fragments, peaks, method = c("ratio", "mann-whitney"))
fragments |
List of MassArrayFragment objects corresponding to the sample |
peaks |
List of MassArrayPeak objects comprising spectral data for a complete assay |
method |
Specifies which algorithm to use when estimating primer dimer levels (either |
Primer dimers are calculated by: 1) identifying fragments that occur within the expected primer sequence, 2) identifying which of these fragments is assayable, and 3) comparing the overall signal for primer peaks and peaks from the rest of the amplicon.
Returns a list containing primer dimer ratios or significance estimates (i.e. p-values) depending on the analytical method specified ("ratio"
or "mann-whitney"
, respectively). Returns "NA"
in cases where insufficient data is present to calculate primer dimer levels.
Reid F. Thompson (rthompso@aecom.yu.edu), John M. Greally (jgreally@aecom.yu.edu)
See Also MassArrayData
data(MassArray.example.data) primer.data <- estimatePrimerDimer(MassArray.example.data$fragments.T, MassArray.example.data$samples[[1]]$peaks, "ratio") mean(primer.data) primer.data <- estimatePrimerDimer(MassArray.example.data$fragments.T, MassArray.example.data$samples[[1]]$peaks, "mann-whitney") mean(primer.data)