r - 如何使用summary()让R按级别名称标记因子？

Question

我使用汇总进行线性回归，通常汇总用因子水平的名称标记系数。但是，当我使用 options(contrasts=c("contr.sum","contr.poly")) 时,
总结出来的时候用变量名 1 和变量名 2 来标记系数，而不是级别的名称。我怎样才能让它显示级别名称？

如果我必须手动完成，我怎么知道哪个系数对应哪个级别？

Answer 1

我会为每个因素手动定义对比，如果需要，还可以定义标签，如下所示:

data(iris)
contrasts(iris$Species)
#>            versicolor virginica
#> setosa              0         0
#> versicolor          1         0
#> virginica           0         1
lm(Sepal.Length ~ Species, data = iris)
#> 
#> Call:
#> lm(formula = Sepal.Length ~ Species, data = iris)
#> 
#> Coefficients:
#>       (Intercept)  Speciesversicolor   Speciesvirginica  
#>             5.006              0.930              1.582
contrasts(iris$Species) <- contr.sum(3)
colnames(contrasts(iris$Species)) <- c("_setosa", "_versicolor")
contrasts(iris$Species)
#>            _setosa _versicolor
#> setosa           1           0
#> versicolor       0           1
#> virginica       -1          -1
lm(Sepal.Length ~ Species, data = iris)
#> 
#> Call:
#> lm(formula = Sepal.Length ~ Species, data = iris)
#> 
#> Coefficients:
#>        (Intercept)      Species_setosa  Species_versicolor  
#>            5.84333            -0.83733             0.09267