我知道
glm(formula = y ~ x1 + x2 + x3, family = binomial(logit), data = birds_poop)
将为x
{x1
,x2
,x3
}提供参数估算值。
如果x的级别超过2级,如何设置参考级别以及如何使glm
函数显示每个x级别的参数估计值。
例如,假设x1有4个级别,我如何获得 x1-level2 , x1-level3 和 x1-level4 假设 x1-level0 是参考水平。
答案 0 :(得分:0)
R默认使用治疗对比。您可以将这些更改为其他类型的对比,通过查看help(contrasts)
和help(C)
可以找到更多信息。对于默认设置,(Intercept)
项将是x1-level0和所有其他变量在其基本级别或0为数值变量的情况下的日志(赔率)等级的估计值。 R将报告不是参考水平的因子水平的参数估计。一个常见的问题是如何更改使用factor()函数轻松完成的参考级别。您也可以使用predict
。
您应该阅读glm
的帮助文件和所有链接。帮助页面的另一个奇妙特征是它们通常包含有用的示例,如果您有其他问题并且难以制定一个好的示例,那么您可能会发现修改其中一个将使讨论更具体。