我进行了两项荟萃分析,并希望证明两种荟萃分析的平均效应大小(在渔民z中)不同。
由于我在R中退出了新的而不是统计学方面的专家,你能否提供适当的测试以及如何在R中进行测试?
以下是我目前对两个元分析的结果:
> results1GN
Random-Effects Model (k = 4; tau^2 estimator: REML)
tau^2 (estimated amount of total heterogeneity): 0.0921 (SE = 0.0752)
tau (square root of estimated tau^2 value): 0.3034
I^2 (total heterogeneity / total variability): 99.98%
H^2 (total variability / sampling variability): 5569.05
Test for Heterogeneity:
Q(df = 3) = 22183.0526, p-val < .0001
Model Results:
estimate se zval pval ci.lb ci.ub
0.3663 0.1517 2.4139 0.0158 0.0689 0.6637 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> results1NN
Random-Effects Model (k = 72; tau^2 estimator: REML)
tau^2 (estimated amount of total heterogeneity): 0.0521 (SE = 0.0096)
tau (square root of estimated tau^2 value): 0.2282
I^2 (total heterogeneity / total variability): 95.98%
H^2 (total variability / sampling variability): 24.85
Test for Heterogeneity:
Q(df = 71) = 1418.1237, p-val < .0001
Model Results:
estimate se zval pval ci.lb ci.ub
0.2594 0.0282 9.2016 <.0001 0.2042 0.3147 ***
答案 0 :(得分:-1)
我详细说明了之前的回复,这有点粗鲁。
在第一步,您可以考虑结果的置信区间。 您可以查阅该主题的维基百科页面以获得快速摘要: http://en.wikipedia.org/wiki/Confidence_interval
无论如何,在您的示例中,您有两个具有重叠置信区间的效果大小:
结果1GN = 0.36(95%CI:0.069至0.66)
结果1NN = 0.26(95%CI:0.20至0.31)。
所以,这两个结果没有统计学差异。
包'metafor'还包含一个函数(anova)来比较“嵌套模型”,我引用: “对于”rma.uni“类的两个(嵌套)模型,该函数在模型拟合统计和似然比检验方面提供了完整与简化的模型比较。
注意“当指定两个模型进行比较时,该功能提供了比较两个模型的似然比检验。模型必须基于相同的数据集,并且应该嵌套以使似然比检验有意义”。
您可以在此处找到有关此功能使用的更多详细信息和示例: http://www.inside-r.org/packages/cran/metafor/docs/anova.rma.uni
另外,请参阅软件包的手册和软件包维护者的站点。
希望这可能有所帮助。