### R code from vignette source '4-logistic.Rnw' ################################################### ### code chunk number 1: readData ################################################### X <- data.frame(y=c(178, 138, 108, 570, 648, 442, 138, 252, 252), belief=rep(c("1-Fundam", "2-Moder", "3-Liber"), 3), degree=rep(c("1-