使用{metafor}进行meta-analisis中的模糊限制器估计和随机效应

时间:2015-07-28 09:49:10

标签: r

我使用对数响应率(“ROM”)作为效果大小,使用Hedges Estimator(“HE”)进行metafor:

rma(es, var, mods=~ Myc * N , method="HE", data=dat, weights=weightsTime)

现在我想要包含随机效果,但rma.mv仅允许method=REMLmethod=ML。这将是我想要的模型:

rma.mv(es, var, mods=~ Myc * N, random = list(~ Myc|SPS.in.SITE, ~N|SPS.in.SITE), method="HE", data=dat, W=weightsTime)

method="HE"不是有效选项。 如何使用Hedges Estimator包含随机效果? 感谢

1 个答案:

答案 0 :(得分:0)

你做不到。 rma.mv()只进行ML / REML估算。其他估算器,例如HE,DL等(可用于rma.uni())尚未推广到rma.mv()函数处理的所有复杂数据结构。关于(受限制的)最大似然估计的好处是它是一个通用框架,只要你可以记下似然函数就可以工作。另外,ML / REML估计器具有良好且经过充分研究的统计特性。