## create a 'modlist'
ml1 <- modlist(reps, 1, 2:5, l4)
## convert into 'replist'
rl1 <- replist(ml1, group = rep(1, 4))
par(mfrow = c(2, 1))
plot(rl1, main = "A 'replist'")
## convert back to a 'modlist'
ml2 <- rep2mod(rl1)
plot(ml2, main = "A 'modlist'")
Run the code above in your browser using DataLab