将最后一级设置为所有回归分析的参考类别

Question

默认情况下，R将因子变量的 first 级别设置为回归中的参考类别。我知道如何使用relevel来指定另一个级别，但我想知道是否有设置或其他东西所以我可以全局强制R将 last 级别的因子设置为参考类别？

Answer 1

是的，因为那是SAS默认值，所以contr.treatment有一个包装器将最后一个级别设置为参考级别：

options(contrasts = c(unordered = "contr.SAS", ordered = "contr.poly"))
fit <-lm(Sepal.Length ~ Species, data = iris)
coef(fit)
#(Intercept)     Speciessetosa Speciesversicolor 
#      6.588            -1.582            -0.652 
levels(iris$Species)
#[1] "setosa"     "versicolor" "virginica"