## alldata from session 3 table(alldata$snp3) alldata$add <- with(alldata, (snp3=="TC")+2*(snp3=="TT")) alldata$dom <- with(alldata, (snp3=="TC") | (snp3=="TT")) alldata$rec <- with(alldata, snp3=="TT") alldata$df2 <- factor(alldata$snp3) names(alldata) addlm<-lm(sbp~add, data=alldata) addlm coef(addlm) confint.default(addlm) summary(addlm) domlm<-lm(sbp~dom, data=alldata) coef(domlm) confint.default(domlm) summary(domlm) adjlm<-lm(sbp~add+dbp+bmi, data=alldata) summary(adjlm) confint.default(adjlm) alldata$hibp<-alldata$sbp>140 addglm<-glm(hibp~add, data=alldata,family=binomial) coef(addglm) confint.default(addglm) summary(addglm) exp(coef(addglm)) exp(confint.default(addglm)) adjglm<-glm(hibp~add+dbp+bmi, data=alldata, family=binomial) summary(adjglm) confint.default(adjglm)