snp.data = dbp[,c("affection", "rs1112")]
summary(snp.data)
snp.data[,"rs1112"]<-as.numeric(snp.data[,"rs1112"])-1<br />summary(snp.data)<br />result.all=glm(affection ~ rs1112, family=binomial("logit"), data=snp.data)<br />dev.all = anova(result.all, test="Chi")<br />print(dev.all)<br />snp.data = dbp[,c("affection","trait","sex","age","rs1112","rs1117")]<br />summary(snp.data)<br />snp.data[,"rs1112"]<-as.numeric(snp.data[,"rs1112"])-1<br />snp.data[,"rs1117"]<-as.numeric(snp.data[,"rs1117"])-1<br /><br />
Multifactorial Part 2
result.reg = glm (affection ~ sex + age + rs1112 + rs1117, family=binomial("logit"), data=snp.data)
summary(result.reg)
modelchoice.result<-step(result.reg)<br />summary(modelchoice.result)
GWAS Control Substructure