com = c(0,5,3,6,1,9,0,0,4,1,0,0,4,0,2,0,0,0,2)
sem = c(9,4,8,3,8,2,9,9,8,10,1,2,0,2,0,2,4,2,0)
n.tot = c(340.52,344.07,147.43,411.64,264.75,263.24,1087.00,1028.32,697.57,
428.56,414.46,538.18,519.69,590.99,472.54,381.82,465.55,528.43,334.95)
D.L70 = c(4.44,3.92,6.02,3.80,3.23,4.56,8.53,9.56,2.29,3.61,14.99,5.85,
3.96,15.30,3.96,1.01,1.32,1.79,4.51)
GLM = glm(cbind(com,sem)~n.tot+D.L70,family = binomial)
summary(GLM)