## create l5 model and plot residuals
m1 <- pcrfit(reps, 1, 2, l5)
resplot(m1)
## compare to l4 model (added smaller bars)
m2 <- pcrfit(reps, 1, 2, l4)
resplot(m2, add = TRUE, width = 0.5, space = c(1.4, 0.9))
Run the code above in your browser using DataLab