我使用对数响应率(“ROM”)作为效果大小,使用Hedges Estimator(“HE”)进行metafor:
rma(es, var, mods=~ Myc * N , method="HE", data=dat, weights=weightsTime)
现在我想要包含随机效果,但rma.mv
仅允许method=REML
或method=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包含随机效果?
感谢
答案 0 :(得分:0)
你做不到。 rma.mv()
只进行ML / REML估算。其他估算器,例如HE,DL等(可用于rma.uni()
)尚未推广到rma.mv()
函数处理的所有复杂数据结构。关于(受限制的)最大似然估计的好处是它是一个通用框架,只要你可以记下似然函数就可以工作。另外,ML / REML估计器具有良好且经过充分研究的统计特性。