我有一个关于在森林图中显示稳健估计的快速问题。
我使用metafor
包中的rma.mv()
来拟合多层随机效应模型。为了进行稳健估计,我使用clubSandwich
包中的conf_int()
和coef_test()
函数。然而,metafor
的forest()
函数只接受rma
类的对象,这些对象不是由函数产生的稳健估计。
是否有办法使用forest()
函数来创建一个漂亮的森林图,但以某种方式整合clubSandwich
函数的稳健总体估计?
感谢任何建议或替代解决方案!
您可以使用forest()
绘制单个估计,然后使用addpoly()
添加多边形。对于后者,您只需要从coef_test()
返回的对象中提取池估计和相应的标准误差。