pdf("race7579.pdf") race1 <- read.fwf("C:\\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") dev.off()