Question

是否有人能够解释如何为r中的不同级别的分类变量设置常量？

我已阅读以下内容：How to set the Coefficient Value in Regression; R并且它可以很好地解释如何为整个分类变量设置常量。我想知道如何为每个级别设置一个。

举个例子，让我们看看MTCARS数据集：

df <- as.data.frame(mtcars)

df$cyl <- as.factor(df$cyl)
set.seed(1)
glm(mpg ~ cyl + hp + gear, data = df)

这给了我以下输出：

Call:  glm(formula = mpg ~ cyl + hp + gear, data = df)

Coefficients:

(Intercept)         cyl6         cyl8           hp         gear
     19.80268     -4.07000     -2.29798     -0.05541      2.79645  

Degrees of Freedom: 31 Total (i.e. Null);  27 Residual

Null Deviance:      1126 

Residual Deviance: 219.5    AIC: 164.4

如果我想将cyl6设置为-.34并将cyl8设置为-1.4，然后重新运行以查看它如何影响其他变量，我该怎么做？

Answer 1

我认为这就是你能做的事情

df$mpgCyl=df$mpg
df$mpgCyl[df$cyl==6]=df$mpgCyl[df$cyl==6]-0.34
df$mpgCyl[df$cyl==8]=df$mpgCyl[df$cyl==8]-1.4
model2=glm(mpgCyl ~ hp + gear, data = df)
> model2

Call:  glm(formula = mpgCyl ~ hp + gear, data = df)

Coefficients:
(Intercept)           hp         gear  
   16.86483     -0.07146      3.53128

更新评论：

cyl是一个因素，因此默认情况下它会将glm作为偏移而不是斜率。实际上cyl==4是隐藏的＆＃39;但也存在于glm中。因此，在您的第一个模型中，模型所说的是：

1) for cyl==4: mpg=19.8-0.055*hp+2.79*gear
2) for cyl==6: mpg=(19.8-4.07)-0.055*hp+2.79*gear
3) for cyl==8: mpg=(19.8-2.29)-0.055*hp+2.79*gear

也许您也可以在这里查看https://stats.stackexchange.com/questions/213710/level-of-factor-taken-as-intercept和Is there any way to fit a `glm()` so that all levels are included (i.e. no reference level)?

希望这有帮助

在R中为分类变量设置不同的常量级别

1 个答案: