logreg<-read.table("logistic1.txt",header=T)
attach(logreg)
out1<-glm(cbind(Yes, No)~dust+race+sex+smoking+Empleng,family=binomial)
out1
summary(out1)
anova(out1)

out2<-glm( cbind(Yes, No)~dust+smoking+Empleng, family=binomial)
anova(out2,out1)
pchisq(3.053,2)
summary(out2)

residuals(out2, type="d")
residuals(out2, type="pear")
plot(residuals(out2, type="d"), xlab="Index", ylab="Deviance Residuals")
abline(h=0)
plot(residuals(out2, type="pear"), xlab="Index", ylab="Pearson Residuals")
abline(h=0)

out3<-glm( cbind(Yes, No)~dust+smoking+Empleng,family=binomial(link=probit))
out3
summary(out3)