ls()
dbp[1:5,]
result.snp12 = glm (affection ~ rs1112, family=binomial("logit"), data=dpbdbp)
print (result.snp12)
print ( class (result.snp12) )
print (ci)
print ( exp(ci) )
snp.data = dpbdbp[,c("affection", "rs1112")]
summary(snp.data)
snp.data[,"rs1112"] summary(snp.data)
summary(result.all)
print(dev.all)
snp.data = dpbdbp[,c("affection", "trait","sex", "age", "rs1112", "rs1117")]
summary(snp.data)
snp.data[,"rs1112"] snp.data[,"rs1117"] result.adj = glm (affection ~ sex + rs1112 , family=binomial("logit"), data=snp.data)