Gene expression from real-time quantitative PCR

Gene expression levels from real-time quantitative polymerase chain reaction (qPCR) experiments on two different plant lines. Each line was used for 7 experiments each with 45 cycles.

Format

A data frame with 630 observations on the following 4 variables.

flour	numeric	Fluorescence level
line	factor	Plant lines rnt (mutant) and wt (wildtype)
cycle	numeric	Cycle number for the experiment
transcript	factor	Transcript used for the different runs

Source

Data provided by Kirsten Jorgensen <kij@life.ku.dk>. Added by Claus Ekstrom <ekstrom@life.ku.dk>

References

Morant, M. et al. (2010). Metabolomic, Transcriptional, Hormonal and Signaling Cross-Talk in Superroot2. Molecular Plant. 3, p.192--211.

Examples

data(qpcr)

#
# Analyze a single run for the wt line, transcript 1
#
run1 <- subset(qpcr, transcript==1 & line=="wt")

model <- nls(flour ~ fmax/(1+exp(-(cycle-c)/b))+fb,
             start=list(c=25, b=1, fmax=100, fb=0), data=run1)

print(model)
#> Nonlinear regression model
#>   model: flour ~ fmax/(1 + exp(-(cycle - c)/b)) + fb
#>    data: run1
#>       c       b    fmax      fb 
#> 29.5932  0.8406 96.7559  3.7226 
#>  residual sum-of-squares: 53.79
#> 
#> Number of iterations to convergence: 9 
#> Achieved convergence tolerance: 7.746e-06

plot(run1$cycle, run1$flour, xlab="Cycle", ylab="Fluorescence")
lines(run1$cycle, predict(model))