pcropt1 {qpcR}R Documentation

Combinatorial elimination of plateau and ground phase cycles

Description

The estimation of PCR efficiency and calculation of initial fluorescence (F0) is analyzed by refitting the (optimized) model on subsets of the data, thereby using all possible combinations of datapoints. The estimated parameters are then collated in a dataframe. This is intended to be the prerequisite for finding the optimal datapoints that minimize the fit or exhibit the best correlation to a calibration curve.

Usage

pcropt1(object, fact = 3, opt = FALSE, plot = TRUE, ...)

Arguments

object an object of class 'pcrfit'.
fact numeric. The multiplier for the scan border. See 'Details'.
opt logical. If true, model selection is applied for each combination of cycles. Beware: Slow!!
plot logical. If TRUE, boxplots of the resulting values are displayed.
... other parameters to be passed on to efficiency, i.e. where to take the PCR efficiency from.

Details

It has been shown by Rutledge (2004) that the estimation of PCR efficiency gives more realistic values when the number of plateau cycles are decreased. This paradigm is the basis for this function, but we also consider the cycles in the ground phase and all combinations between ground/plateau cycles. All datapoints between the lower border cpD1 - fact * (cpD1 - cpD2) and upper border cpD1 + fact * (cpD1 - cpD2) are cycled through.

Value

A dataframe with the border values, AIC, AICc, residual variance, efficiency and the estimated F(0) from the exponential and sigmoidal model.

Author(s)

Andrej-Nikolai Spiess

References

Sigmoidal curve fitting redefines quantitative real-time PCR with the prospective of developing automated high-throughput applications. Rutledge RG, Nucleic Acids Research, 2004, e178.

See Also

The function efficiency that is called by this function.

Examples

## Using one model throughout
## Not run: 
m <- pcrfit(reps, 1, 2, l4)
pcropt1(m)
## End(Not run)

## Selecting the best model for each combination
## by Akaike weights
## Not run: 
m <- pcrfit(reps, 1, 2, l4)
pcropt1(m, opt = TRUE, crit = "weights") 
## End(Not run)

[Package qpcR version 1.2-4 Index]