### Using one model throughout
m <- pcrfit(reps, 1, 2, l5())
temp <- pcropt1(m)
edit(temp)
### Selecting the best model for each combination
m <- pcrfit(reps, 1, 2, l4())
temp <- pcropt1(m, opt = TRUE)
edit(temp)
Run the code above in your browser using DataLab