meta分析 2. 固定效应和随机效应

固定效应模型

固定效应模型假设所有效应量都来自一个单一的相同总体。它指出所有研究都具有相同的真实效应量。这个真实的效果是我们想要在meta分析中计算的整体效果大小,表示为 \theta。

我们可以这样描述这种关系:

\hat\theta_k = \theta + \epsilon_k

虽然不知道研究的真实总体效应大小,但可以利用这种关系来对真实总体效应进行最佳估计 \hat\theta_k, 因此,在meta分析中汇集效应时,应该赋予具有更高精度(即更小的标准误差)的效应量更大的权重

w_k = \frac{1}{s^2_k}

想计算固定效应模型下的合并效应大小,只需使用所有研究的加权平均值。一旦我们知道权重,我们就可以计算加权平均值,我们对真实合并效应的估计 :

\hat\theta = \frac{\sum^{K}_{k=1} \hat\theta_kw_k}{\sum^{K}_{k=1} w_k}

这种方法是在meta分析中计算平均效应的最常用方法。因为我们使用方差的倒数,所以它通常被称为逆方差加权或简称为逆方差荟萃分析✍。

随机效应模型

固定效应模型假设我们所有的研究都是相同总体的一部分,并且观察到的效应差异的唯一原因是研究的抽样误差。如果我们要在没有抽样误差的情况下计算每项研究的效应量,那么所有真实的效应量将完全相同。

固定效应模型的假设在许多实际应用中可能过于简单💣。meta分析中的研究总是完全相同的,这是不现实的。研究往往会有所不同,即使只是微妙的方式。感兴趣的结果可能以不同的方式衡量。也许治疗的类型或治疗的强度和时间不完全相同。研究的目标人群可能并不完全相同,或者所使用的对照组可能存在差异。

meta分析中的研究很可能不仅在这些方面中的一个方面有所不同,而且在几个方面也有所不同。如果这是真的,我们可以预期真实效果的研究之间存在相当大的异质性

所有这些都使固定效应模型的有效性受到怀疑。例如,如果一些研究使用了不同类型的治疗,那么一种形式比另一种更有效似乎是完全正常的。假设这些差异只是由研究的抽样误差产生的噪音,那将是牵强附会。

恰恰相反,研究的真实效果大小存在真正差异的原因可能有无数种。随机效应模型解决了这个问题。它为我们提供了一个模型,通常可以更好地反映数据背后的现实🕶。

在随机效应模型中,我们想考虑这样一个事实,即效应量显示出比从单个相同总体中提取时更大的方差。因此,我们假设个别研究的效果不仅因抽样误差而产生偏差,而且还有另一个方差来源。这一额外的方差分量是由于研究并非源于一个单一群体这一事实而引入的。相反,每项研究都被视为来自人群“宇宙”的独立抽取。

与固定效应模型类似,随机效应模型首先假设观察到的效应大小 \hat\theta_k 是研究真实效应大小 \theta_k 的估计值, 受抽样误差的影响 \epsilon_k 。

\hat\theta_k =\theta_k+\epsilon_k

\hat\theta_k 代替 \theta 已经指出了一个重要的区别。随机效应模型仅假设\hat\theta_k是一项单一研究 k 的真实效应量. 它规定存在第二个误差源,表示为 \epsilon_k。第二个误差源是由以下事实引入的:即使是真实的效果大小 \theta_k 只是真实效应大小总体分布的一部分,均值 \mu 。

\theta_k = \mu + \zeta_k

通过将第二个公式插入第一个公式,随机效应模型:

\hat\theta_k =\mu + \zeta_k+ \epsilon_k

这个公式清楚地表明观察到的效应大小 \hat\theta_k 由于 \zeta_k 和 \epsilon_k 偏离了合并效应,这种关系如图所示:

随机效应模型.png

因此,在包括医学和社会科学在内的许多领域中,通常总是使用随机效应模型,因为几乎总是可以预料到一定程度的研究间异质性。只有当我们无法检测到任何研究间的异质性,并且我们有很好的理由假设真正的效应是固定的时,才可以使用固定效应模型。例如,当仅考虑一项研究的精确复制时,或者当对一项大型研究的子集进行meta分析时,可能就是这种情况。毋庸置疑,这种情况很少发生,固定效应模型“在野外”的应用也相当罕见。

研究间异质性的估计

与随机效应模型相关的挑战是我们必须接受错误 \zeta_k 。为此,我们必须估计真实效应量分布的方差。这个方差被称为 \tau^2 或tau 平方。一旦我们知道了 \tau^2 的值,我们可以在确定每个效应大小的反方差权重时包括研究间的异质性。

在随机效应模型中,调整后的随机效应权重 w^*_k对于每个观察。公式如下所示:

w^*_k = \frac{1}{s^2_k+\tau^2}

使用调整后的随机效应权重,我们然后使用逆方差方法计算合并效应大小,就像我们使用固定效应模型一样:

有几种估计 \tau^2 的方法:

  • DerSimonian-Laird ("DL")

  • Restricted Maximum Likelihood ("REML") or Maximum Likelihood ("ML")

  • Paule-Mandel ("PM")

  • Empirical Bayes ("EB")

  • Sidik-Jonkman ("SJ")

选择估算方法

确切地说应该使用哪个估计没有铁一般的规则。在很多情况下,各种估计产生的结果只会有细微的差别💬。

  1. 对于基于连续结果数据的效果大小,可以使用 ML 作为第一个开始。

  2. 对于二元效应大小数据,Paule-Mandel 估计是不错的首选,前提是样本大小没有极端变化。

  3. 当有充分的理由相信样本中效应的异质性非常大,并且如果避免误报具有非常高的优先级,使用 Sidik-Jonkman 估计。

  4. 如果希望其他人可以在R之外尽可能精确地重复的结果,选择DerSimonian-Laird 估计。

参考书籍:

  1. Doing Meta-Analysis in R.
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,684评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 87,143评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,214评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,788评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,796评论 5 368
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,665评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,027评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,679评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 41,346评论 1 299
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,664评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,766评论 1 331
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,412评论 4 321
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,015评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,974评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,203评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,073评论 2 350
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,501评论 2 343

推荐阅读更多精彩内容