Changes

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 />result.adj = glm(affection ~ sex + rs1112, family=binomial("logit"), data=snp.data)<br />summary(result.adj)