q*********i 发帖数: 696 | 1 有一组数据,k是trial的次数,y是success的次数,pred是唯一的predictor
每个observation的trial次数都不一样。
R代码写成下面这样可以吗?
glm (y/k~pred,binomial(link = "logit"), weights =k)
或
glm (y/k~pred,binomial(link = "logit"))
看了半天help也没看明白weights参数在这里的意思,试着fit了一下两个命令得到的参
数差不多。deviance和AIC都是有weight的那个小一些。 |
|