pdf("race7579.pdf") race1 <- read.fwf("http://www.stat.ufl.edu/~winner/data/race7579.dat", width=c(8,8,8,8,8,8,8,8,8,12,40), col.names=c('srace', 'yr', 'yrace', 'drivers', 'trklen', 'laps', 'roadtrk', 'cautions', 'leadchng', 'trkid', 'track')) race <- data.frame(drivers=race1$drivers, trklen=race1$trklen, laps=race1$laps, cautions=race1$cautions) attach(race) race.mod <- glm(formula = cautions ~ drivers + trklen + laps, family=poisson("log")) summary(race.mod) anova(race.mod, test="Chisq") muhat <- predict(race.mod, type="response") print(cbind(cautions, muhat)) (pearson.x2 <- sum((cautions - muhat)^2/muhat)) dev.off()