pdf("westnilesc.pdf")

wnsc <- read.fwf("C:\\data\\westnilesc.dat", 
width=c(24,3,8,8,8,8,8,12,12),
col.names=c('cnty','wnbird', 'wnequine', 'farms', 'area', 
 'pop', 'humdens', 'posbrdrt', 'poseqnrt'))

wnsc1 <- data.frame(wnsc)

attach(wnsc1)

wnsc.mod1 <- glm(formula = wnequine ~ offset(log(farms)) + 
 posbrdrt + humdens + posbrdrt:humdens, family=poisson(log))

summary(wnsc.mod1)




dev.off()