Question

如何打印将分类/名义变量输入线性模型时使用的参考类别。这是一个例子：

summary(lm(data = iris, Sepal.Length ~ Species))

Coefficients:
                  Estimate Std. Error t value Pr(>|t|)    
(Intercept)         5.0060     0.0728   68.76  < 2e-16 ***
Speciesversicolor   0.9300     0.1030    9.03  8.8e-16 ***
Speciesvirginica    1.5820     0.1030   15.37  < 2e-16 ***

这是我想要的：

                  Estimate Std. Error t value Pr(>|t|)    
(Intercept)         5.0060     0.0728   68.76  < 2e-16 ***
**Reference: Speciessetosa**
Speciesversicolor   0.9300     0.1030    9.03  8.8e-16 ***
Speciesvirginica    1.5820     0.1030   15.37  < 2e-16 ***

如果有一种方法可以使这项工作普遍存在（当有多个分类预测因子时，那么每个参考组都很容易识别），这将是最优秀的。如果有一种方法可以使格式化特别清晰，那将是非常好的（我不喜欢上面的示例格式）。

Answer 1

您可以指定截距为零。

summary(lm(Sepal.Length ~ Species + 0, data = iris))
#Coefficients:
#                  Estimate Std. Error t value Pr(>|t|)    
#Speciessetosa       5.0060     0.0728   68.76   <2e-16 ***
#Speciesversicolor   5.9360     0.0728   81.54   <2e-16 ***
#Speciesvirginica    6.5880     0.0728   90.49   <2e-16 ***

在R的lm摘要中打印（显示）参考类别？

1 个答案: