data(HairEye)
haireye.crg <- corregp(Eye ~ Hair * Sex, data = HairEye, b = 3000)
haireye.crg
summary(haireye.crg, parm = "b", add_ci = TRUE)
screeplot(haireye.crg, add_ci = TRUE)
plot(haireye.crg, x_ell = TRUE, xsub = c("Hair", "Sex"))Run the code above in your browser using DataLab