微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

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

如何解决使用 metafor 的多级元分析 - 代码

我是 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 举报,一经查实,本站将立刻删除。