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)
plot(run1$cycle, run1$flour, xlab="Cycle", ylab="Fluorescence")
lines(run1$cycle, predict(model))
Run the code above in your browser using DataLab