R语言中固定与随机效应Meta分析 - 效率和置信区间覆盖
原文链接:http://tecdat.cn/?p=6330
Meta分析是综合现有证据的重要工具。它通常用于医学和临床环境中,以评估关于治疗或暴露对感兴趣结果的影响的现有证据。基本思想是将先前研究中感兴趣的影响的估计汇总在一起。在进行Meta分析 时必须做出的选择是固定效应和随机效应之间的选择。在这篇文章中,我们将看看这种选择的一些后果,实际上研究正在测量不同的效果。
设置
我们假设我们对n个研究对我们感兴趣的影响进行了估计。
固定效应Meta分析
在固定效应荟萃分析中,我们假设所包含的每项研究都在估计相同的基础参数。
随机效应Meta分析
相比之下,在随机效应荟萃分析中,我们假设每项研究都在估计研究特定的真实效果。观察到的估计异质性归因于两个来源:1)研究之间的真实效应异质性,以及2)研究内抽样误差。
固定与随机效应Meta分析
我们使用哪种方法会影响我们获得的估计总体效果及其相应的95%置信区间,因此决定哪种方法适合在任何给定情况下使用是很重要的。我个人认为,这个决定应该基于对成分研究的了解,而不是基于实际查看点估计。
两种方法之间的另一个明显差异是,固定效应荟萃分析的计算标准误差小于随机效应荟萃分析的计算标准误差。这可能会使研究人员相信固定效应估计更精确。
R中的模拟研究
为了进行模拟研究,我们将模拟30项研究的重复荟萃分析。
然后,我们执行固定效应和随机效应荟萃分析。代码如下所示:
fixedEffectEsts < - array(0,dim = c(nSims))
fixedCI < - array(0,dim = c(nSims,2))
randomEffectEsts < - array(0,dim = c(nSims))
randomCI < - array(0,dim = c(nSims,2))
for(i in 1:nSims){
trueEffects < - rnorm(nStudies,mean = 1,sd = 0.1)
#sample在研究标准偏差范围内
withinStudySD < - rgamma(nStudies,shape = 2.5,scale = 0.04)
studyEstimate < - rnorm(nStudies,mean = trueEffects,sd = withinStudySD)
#fixed effects meta-analysis
maFixed < - (studyEstimate~1,S = withinStudySD ^ 2,method =“fixed”)
fixedEffectEsts [i] < - coef(maFixed)
fixedCI [i,] < - c(coef(maFixed)-1.96 * maFixed $ vcov ^ 0.5,coef(maFixed)+ 1.96 * maFixed $ vcov ^ 0.5)
#random-effects元分析
maRandom < - (studyEstimate~1,S = withinStudySD ^ 2,method =“mm”)
randomEffectEsts [i] < - coef(maRandom)
randomCI [i,] < - c(coef(maRandom)-1.96 * maRandom $ vcov ^ 0.5,coef(maRandom)+ 1.96 * maRandom $ vcov ^ 0.5)
}
ci覆盖范围
mean((fixedCI [,1] <1)&(fixedCI [,2]> 1))
mean((randomCI [,1] <1)&(randomCI [,2]> 1))
R代码的最后几行计算1000次模拟中固定和随机效应估计的平均值和SD,然后计算95%置信区间的覆盖范围。当我运行我获得的代码时:
> mean(fixedEffectEsts)
[1] 0.9990649
> sd(fixedEffectEsts)
[1] 0.04939415
>
> mean(randomEffectEsts)
[1] 1.000967
> sd(randomEffectEsts)
[1] 0.0242558
>
> #ci
> mean((fixedCI [,1] <1)&(fixedCI [,2]> 1))
[1] 0.322
> mean((randomCI [,1] <1)&(randomCI [,2]> 1))
[1] 0.926
首先要注意的是,固定效应方法仍然是无偏的,即使数据是基于随机效应模型进行模拟的。然而,我们发现固定效应方法的SD要大得多(随机效应为0.049对0.024)。或者换句话说,随机效应估计器是一种更精确的估计。
对这种明显矛盾的解释是,当实际存在研究异质性之间时,由固定效应方法计算的标准误差是无效的。
结论
我从这个小型模拟研究中得出的结论是,应该警惕使用固定效应分析,除非人们相信meta分析中的研究正在估计相同的共同效应。
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 如何优化Nginx的处理性能
- 一键生成前后端代码,Mybatis-Plus代码生成器让我舒服了
- 各类好玩免费API推荐,强烈建议收藏
- 回敬Python蹭女神热度助发朋友圈,Java实现头像分成九宫图,再生成拼图头像
- java作业:实现数组翻转
- windows下搭建spark测试环境
- 微信小程序实时语音识别实践
- 轻量安全的部署方案
- 阿里2020.7.31笔试
- 聊聊dubbo-go的PrometheusReporter
- 聊聊dubbo-go的apolloConfiguration
- 看了这篇,关于浏览器缓存你还有哪些疑问?
- 深入分析Volatile的实现原理
- Java Collection Framework : List
- String类的equals方法