我正在使用该包检查两个连续预测变量之间的相互作用emmeans。我使用包lm_robust()中的estimatr内容来执行线性回归并获得集群稳健的标准误差。结果变量以 SD 单位方差为中心并缩放。例如:
fit <- lm_robust(scale(Y) ~ X1 * X2 + X3 + X4, data = mydata, cluster = school, se_type = 'CR2')
X2然后,我可以使用类似于以下的代码执行成对对比或可视化三个级别的线条:
emmip(fit, X2 ~ X1, CIs = TRUE, at = list(X2 = c(mean(X2) - sd(X2),
mean(X2),
mean(X2) + sd(X2))))
Run Code Online (Sandbox Code Playgroud)
我不想将结果变量转换为其原始规模。
我的问题是是否emmeans使用集群稳健标准误差来计算其报告的置信区间或 p 值,并且此行为是否取决于结果变量是处于其原始规模还是经过转换?estimatr包创建者网站上的一个简短示例表明lm_robust对象可以与 一起使用,但我在“emmeans 支持的模型”小插图页面或包文档中看emmeans不到将其列为受支持的模型。lm_robust