我想将对照组均值添加到我的回归表中Stargazer
。我目前的想法是用对照组平均值替换截距并重命名该术语。我知道如何替换拦截的值,但是有没有办法重命名拦截项?用对照组平均值替换常数项通常会提供信息。
library(stargazer)
attitude$treatment <- rbinom(n = 30, size = 1, prob = .5)
lm_model <- lm(rating ~ treatment + complaints + privileges + learning
+ raises + critical, data=attitude)
#add in control mean
control_mean <- mean(attitude[attitude$treatment == 0, "rating"])
lm_model[["coefficients"]][names(lm_model[["coefficients"]]) ==
"(Intercept)"] <- control_mean
stargazer(lm_model)