随机对照研究(连续型变量)的Meta分析——R软件实现(meta包计算ROM效应量)

发布于 2024年4月22日星期一 15:45:21 浏览：536

原创不易，转载请注明来源，感谢！

附件下载：

1716363066-随机对照研究(连续型数据).csv 请勿重复点击，如无响应请耐心等待或稍后再试。

在随机对照试验(Randomized controlled trail, RCT)中，某些测量结果是在某一特定范围内的任意取值，每个测量结果都有一个具体的数值，在统计学上称为连续型数据(也叫计量资料)，如身高、体重。对于连续型数据，假设Meta分析中的每项研究按组样本量、均数、标准差提取数据，则可形成六格表资料，其主要的效应量指标有均数差(Difference in means，MD)、标准化均数差(Standardised mean difference，SMD)、均数比(Ratio of means，RoM)、几何均数比(Ratio of geometric means，RoGM)等。本文将实例演示在R软件中使用meta包中的metacont命令计算ROM实现随机对照研究(连续型变量)的Meta分析相关操作。

关键词：R; Meta分析; 随机对照试验; 连续型数据; meta包; 均数比

一、案例介绍

本文以模拟数据为例，为研究脂肪摄入对小型家兔体重的影响，现收集了10篇研究报道，均以小型家兔体重增量(单位g)作为主要结局指标，以此来比较饮食中加入脂肪添加剂与不加入脂肪添加剂的家兔体重是否有统计学差异，数据见表1，案例数据可从“附件下载”处下载。

二、问题分析

对于连续型数据的Meta分析，要根据数据的类型、测量值类型、测量基线是否平衡、效应量等综合考虑。本案例使用的是聚合数据。在研究基线平衡的情况下，可以选择最终值和变化值(如本案例中的体重增量)进行合并，根据异质性检验结果，选择不同的效应模型和方法，如固定效应模型使用IV法，随机效应模型选用D-L法等；如果基线不平衡，可能需要限制使用Meta分析，也可以采用校正方法，但效能比较低。

三、软件操作及结果解读

(一) 数据整理及读取

mydata <- read.csv("随机对照研究(连续型数据).csv")  #导入CSV数据
View(mydata)  #查看数据

在数据栏目中可以查看全部数据情况，数据集中共有8个变量和10个观察数据，8个变量分别代表第一作者(Author)、发表年代(Year)、试验组例数(nt)、试验组均数(mt)、试验组标准差(st)、对照组例数(nc)、对照组均数(mc)、对照组标准差(sc)。

(二) 异质性检验

选择ROM为效应量，任意选用固定或随机效应模型进行异质性检验，可以直接得到Q统计量和I²统计量，再根据异质性检验的结果选择相应的效应模型。

代码如下：

library(meta)
result.mydata <- metacont(n.e = nt, mean.e = mt, sd.e = st, n.c = nc, mean.c = mc, sd.c = sc, data = mydata, sm = "ROM", method.tau = "REML", studlab = paste(Author, Year))
print(summary(result.mydata), digits = 2)

该命令参数介绍如下：