使用 metafor 的多级元分析 - 代码

编程问答 2022-05-23

问题描述

我是 R 的新手。我使用 excel 和综合元分析软件进行了相同的元分析，但我想使用多级元分析来运行模型。我可以遵循网站上的论文中提供的代码，例如 Konstantopoulos,S. (2011)、Assink 和 Wibbelink 的教程论文 (2016) 或 https://bookdown.org/MathiasHarrer/Doing_Meta_Analysis_in_R/fitting-a-three-level-model.html

然而，当我尝试复杂化模型时，我无法正确编码。测试以下内容的正确代码如何：

数据集中的变量： 国号 学习号 效果大小 y 影响大小（相关性） v 方差 出版状态（0 和 1） - 研究等级 年 - 学习等级 设计（0 和 1） - 学习 lvl IGC 国家得分 1 - 国家等级 LPT 国家得分 2 - 国家等级 联合国失业 - 国家等级 RL法治-国家lvl

我想测试控制研究等级因素和国家等级失业率和法治的国家/地区水平是否得分中等/预测。

我可以使用单个预测器运行：

model1 <-rma.mv(y,v,mods = ~ IGC,random = list(~ 1 | CountryID,~ 1 | StudyID),tdist= TRUE,data = data1)

Multivariate Meta-Analysis Model (k = 43; method: REML)

Variance Components:

            estim    sqrt  nlvls  fixed     factor 
sigma^2.1  0.0283  0.1682     14     no  CountryID 
sigma^2.2  0.0102  0.1008     34     no    StudyID 

Test for Residual Heterogeneity:
QE(df = 41) = 365.4720,p-val < .0001

Test of Moderators (coefficient 2):
F(df1 = 1,df2 = 41) = 0.0527,p-val = 0.8195

Model Results:

         estimate      se     tval    pval    ci.lb    ci.ub 
intrcpt   -0.3023  0.0885  -3.4160  0.0014  -0.4810  -0.1236  ** 
IGC       -0.0039  0.0169  -0.2297  0.8195  -0.0379   0.0302

谁能帮我插入其他预测变量并运行完整模型？ 非常感谢，

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

metaformultilevel-analysisrr