r - 在 ggplot 中绘制具有聚集标准误差的回归的置信区间
问题描述
我想绘制一个带有置信区间的回归模型ggplot
。在我的模型中,我想使用聚集在变量上的稳健标准误差。但是,我找不到可以找到聚类错误变量的位置。
我已经尝试过geom_smooth
使用方法的函数,lm_robust
但找不到在哪里添加集群错误变量。
p1 +geom_smooth(data=data,aes(y=y,x=x),method='lm_robust',se = T)
我需要情节中以下行的等价物ggplot
:
lm_robust(y~ x, data = data, clusters = z)
解决方案
似乎非标准评估会带来问题,但您可以像这样传递整个 data.frame 列:
library(ggplot2)
library(estimatr)
ggplot(mtcars, aes(hp, qsec)) +
geom_smooth(method = 'lm_robust', method.args = list(cluster = mtcars$cyl))
请注意,当绘制多条线(例如使用颜色)或使用刻面时,这将不起作用。