当x超过2个级别时,r glm函数用于二项式

时间:2014-02-25 05:55:18

标签: r glm

我知道

glm(formula = y ~ x1 + x2 + x3, family = binomial(logit), data = birds_poop)

将为x {x1x2x3}提供参数估算值。

如果x的级别超过2级,如何设置参考级别以及如何使glm函数显示每个x级别的参数估计值。

例如,假设x1有4个级别,我如何获得 x1-level2 x1-level3 x1-level4 假设 x1-level0 是参考水平。

1 个答案:

答案 0 :(得分:0)

R默认使用治疗对比。您可以将这些更改为其他类型的对比,通过查看help(contrasts)help(C)可以找到更多信息。对于默认设置,(Intercept)项将是x1-level0和所有其他变量在其基本级别或0为数值变量的情况下的日志(赔率)等级的估计值。 R将报告不是参考水平的因子水平的参数估计。一个常见的问题是如何更改使用factor()函数轻松完成的参考级别。您也可以使用predict

您应该阅读glm的帮助文件和所有链接。帮助页面的另一个奇妙特征是它们通常包含有用的示例,如果您有其他问题并且难以制定一个好的示例,那么您可能会发现修改其中一个将使讨论更具体。