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 />result.adj = glm(affection ~ sex + rs1112, family=binomial("logit"), data=snp.data)<br />summary(result.adj)
result.adj = glm (affection ~ age + rs1112 , family=binomial("logit"), data=snp.data)
summary(result.adj)
result.adj = glm (affection ~ sex + age + rs1112, family=binomial("logit"), data=snp.data)
summary(result.adj)
result.adj = glm (affection ~ rs1117 + rs1112, family=binomial("logit"), data=snp.data)
summary(result.adj)
anova (result.adj, test="Chi")
result.adj = glm (affection ~ rs1112 + rs1117, family=binomial("logit"), data=snp.data)
summary(result.adj)
anova (result.adj, test="Chi")
result.adj = lm (trait ~ rs1112, data=snp.data)
summary(result.adj)
result.adj = lm (trait ~ sex + rs1112, data=snp.data)
summary(result.adj)
Multifactorial Part 2