if (FALSE) {
# Compare different types of models
model1 <- cnorm(group = elfe$group, raw = elfe$raw)
model2 <- cnorm.betabinomial(elfe$group, elfe$raw)
model3 <- cnorm.shash(elfe$group, elfe$raw)
# Compare traditional cnorm with ShaSh
compare(model1, model3, age = elfe$group, score = elfe$raw)
# Compare beta-binomial with ShaSh
compare(model2, model3, age = elfe$group, score = elfe$raw)
}
Run the code above in your browser using DataLab