R中使用relevel、fct_relevel后因子参考水平未变更问题求解

• 函数未赋值回原数据集：relevel()、fct_relevel()均为返回新对象的函数，不会原地修改原始数据列。你仅执行了因子修改操作但没有将结果赋值回chickdiet1_3_wide$Diet，建模时读取的仍是未修改的原始因子，因此参考水平保留为默认的1。
• 因子存在冗余水平：你筛选后的数据仅保留了Diet为1、3的观测，但原始因子的水平仍包含["1","2","3","4"]，冗余水平虽然不直接影响参考水平设置，但可能引发其他后续运算问题，建议提前清理。

先清理冗余因子水平，再修改参考水平，且必须将修改结果赋值回原数据列，参考代码如下：

chickdiet1 <- ChickWeight %>% subset(Diet == 1) 
chickdiet1_3 <- rbind(chickdiet1, chickdiet3)
chickdiet1_3_wide <- pivot_wider(chickdiet1_3)

# 核心修改：清理冗余水平+重设参考水平+赋值回原列
chickdiet1_3_wide$Diet <- chickdiet1_3_wide$Diet %>% 
  droplevels() %>% 
  relevel(ref = "3")

summary(lm(weight ~ Diet, data = chickdiet1_3_wide))

修改后拟合的模型截距项对应Diet3的平均体重，系数项为Diet1，代表Diet1相对Diet3的体重差异，符合设置参考水平为3的预期。

如果使用fct_relevel操作，写法逻辑一致，同样需要赋值：

chickdiet1_3_wide$Diet <- fct_relevel(chickdiet1_3_wide$Diet, "3")

内容的提问来源于stack exchange，提问作者chickennugget