## single plot
m1 <- pcrfit(reps, 1, 2, l5)
pcrplot(m1)
## add another plot in blue
## with 99\% prediction interval
m2 <- pcrfit(reps, 1, 12, l5)
pcrplot(m2, add = TRUE, col = 4, confband = "confidence", level = 0.99)
## plot a 'modlist' batch with coloring of replicates
ml <- modlist(reps, 1, 2:13, model = l4)
pcrplot(ml, colvec = gl(3,4))
Run the code above in your browser using DataLab