pdf ("ginkgoams.pdf")

gams0 <- read.fwf("C:\\data\\ginkgoams.dat", width=c(8,8,8,8,8,8),
 col.names=c("trt", "acet", "gb", "illness", "amscount", "hacount"))

gams1 <- data.frame(acet=gams0$acet, gb=gams0$gb, illness=gams0$illness, amscount=gams0$amscount)

attach(gams1)

gams1

gams.mod1 <- glm(illness ~ 1 + acet + gb + acet:gb, binomial("logit"), weights=amscount)

summary(gams.mod1)

dev.off()