pdf("donner.pdf")

don1 <- read.fwf("C:\\data\\donner.dat", width=c(27,3,1,15,12,9,14,9,2),
col.names=c("name", "age", "sex", "survive", "deathdt", "rescue", "joindt", "trapped", "camp"))



male <- as.numeric(don1$sex)

male = male-1

agemale = don1$age*male

donner <- data.frame(age=don1$age, survive=don1$survive, male=male, agemale=agemale)

attach(donner)

donner

donner.mod1 <- glm(formula = survive ~ age + male + agemale, family=binomial("logit"))

summary(donner.mod1)

anova(donner.mod1, test="Chisq")

donner.mod2 <- glm(formula = survive ~ age + male, family=binomial("logit"))

summary(donner.mod2)

anova(donner.mod2, test="Chisq")

dev.off()