pre = c(9,10,8,9,6,7,9,7,8,6,8,4,8,11,6,7,9,9,6,6,9,9,10,9,8)

post = c(13,10,10,12,9,11,12,9,7,9,7,3,7,9,6,10,10,12,8,9,8,9,10,10,7)


y <- c(pre, post)
group <- rep(0:1, each = 25)
summary(model <- glm(cbind(y, 15 - y) ~ group, binomial) )