问题描述
我正在对论文中的相关系数进行元分析。所以 lvl 1 = 论文,lvl 2 = 论文中的系数。
我想在论文和系数级别包含版主,但我不确定如何指定应应用版主的级别。例如,论文层面的不同方法如何影响结果,系数层面的不同参与者如何影响结果。
我目前的模型是:
m1 <- rma.mv(yi,vi,mods = ~ paper_factor + coefficient_factor,random = ~1|paper/id,data=dat)
但是,这并没有专门编码论文之间的 paper_factor 不同,以及论文中的系数之间的 coefficient_factor 不同。我该怎么做?
注意:Metafor
的语法与 lme4
相似,所以这里的适当语法也许也适用?
解决方法
主持人的“级别”取决于其编码方式。如果它是论文中的常数,那么它就是论文级别的主持人。如果版主可以在论文中变化,那么它是估计级别的版主。事实上,这也是 lme4
的运作方式。