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 /> summary(result.all) 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)