skin = read.table("malignant.txt", header=TRUE) skin vaccine = read.table("vaccine.txt", header=TRUE) vaccine ulcer= read.table("ulcer.txt", header=TRUE) ulcer #skin res.min.skin = glm(freq~1, family=poisson, data=skin) summary(res.min.skin) res.add.skin = glm(freq~site+type, family=poisson, data=skin) summary(res.add.skin) res.sat.skin = glm(freq~site*type, family=poisson, data=skin) summary(res.sat.skin) #vaccine #for placebo vaccine$freq[1:3]/sum(vaccine$freq[1:3]) #for vaccine vaccine$freq[4:6]/sum(vaccine$freq[4:6]) #Do not do this, should include offset #res1.vaccine = glm(freq~1, family=poisson, data=vaccine) #summary(res1.vaccine) res2.vaccine = glm(freq~treatment*response,family=poisson, data=vaccine) summary(res2.vaccine) res3.vaccine = glm(freq~treatment+response,family=poisson, data=vaccine) summary(res3.vaccine) ulcer ulcer[1:4,] res1.gulcer = glm(freq ~ case.control + aspirin, family=poisson, data=ulcer[1:4,]) summary(res1.gulcer) res2.gulcer = glm(freq ~ case.control*aspirin, family=poisson, data=ulcer[1:4,]) summary(res2.gulcer) res1.dulcer = glm(freq ~ case.control + aspirin, family=poisson, data=ulcer[5:8,]) summary(res1.dulcer) res2.dulcer = glm(freq ~ case.control*aspirin, family=poisson, data=ulcer[5:8,]) summary(res2.dulcer) res1.dlcer = glm(freq ~ case.control*ulcer, family=poisson, data=ulcer) summary(res1.dlcer) res2.ulcer = glm(freq ~ case.control*ulcer + aspirin, family=poisson, data=ulcer) summary(res2.ulcer) res3.ulcer = glm(freq ~ case.control*ulcer + aspirin + aspirin*case.control, family=poisson, data=ulcer) summary(res3.ulcer) res4.ulcer = glm(freq ~ case.control*ulcer + aspirin + aspirin*case.control+aspirin*ulcer, family=poisson, data=ulcer) summary(res4.ulcer) res5.ulcer = glm(freq ~ case.control*ulcer*aspirin, family=poisson, data=ulcer) summary(res5.ulcer) res6.ulcer = glm(freq ~ 1, family=poisson, data=ulcer) summary(res6.ulcer) 1/8*sum((fitted(res4.ulcer)-ulcer$freq)^2/sqrt(fitted(res4.ulcer)))