data <- matrix(scan(),ncol=3,byrow=TRUE) 2.52 38 0 2.46 36 0 2.29 36 0 3.15 36 0 2.88 30 0 2.29 31 0 2.99 36 0 2.38 37 1 2.56 31 0 3.22 38 0 2.35 29 0 3.53 46 1 2.65 46 0 2.15 31 0 3.32 35 0 2.23 37 0 2.19 33 0 2.21 37 0 5.06 37 1 2.68 34 0 2.09 44 1 2.54 28 0 2.18 31 0 2.67 39 0 3.41 37 0 3.15 39 0 3.34 32 1 2.60 38 0 2.28 36 0 3.93 32 1 2.60 41 0 3.34 30 0 data Example1223 <- data.frame(data) colnames(Example1223) <- c("fib","gamma","health") rm(data) Example1223 Example1223.Fit <- glm(health ~ fib + gamma, family=binomial, data=Example1223) anova(Example1223.Fit) summary(Example1223.Fit)