为什么我在 "meta" 包中的荟萃分析中得到不同的 95% 置信区间?
Why I am getting different 95% confidence interval for some studies in meta-analysis in "meta" package in r?
我正在尝试使用风险比、下限和上限 95% 置信区间进行荟萃分析,但例如 CARDIa 研究,获得上限和下限 95%CI ([2.1560; 9.9858]) 是与原始值 (1.33-6.16) 不同,我不知道如何获得确切的数字。
如有任何建议,我们将不胜感激。
使用代码:
library(meta);library(metafor)
data<-read.table(text="studlab HR LCI UCI
Blazek 1.78 0.84 3.76
PRECOMBAT 1.20 0.37 3.93
LE.MANS 1.14 0.3 4.25
NOBLE 2.99 1.66 5.39
MASS-II 2.90 1.39 6.01
CARDIa 4.64 1.33 6.16
BEST 2.75 1.16 6.54
", header=T, sep="")
metagen(log(HR), lower = log(LCI), upper = log(UCI),
studlab = studlab,data=data, sm = "HR")
得到的结果
HR 95%-CI %W(fixed) %W(random)
Blazek 1.7800 [0.8413; 3.7659] 16.4 16.5
PRECOMBAT 1.2000 [0.3682; 3.9109] 6.6 7.1
LE.MANS 1.1400 [0.3029; 4.2908] 5.2 5.7
NOBLE 2.9900 [1.6593; 5.3878] 26.6 25.0
MASS-II 2.9000 [1.3947; 6.0301] 17.2 17.2
CARDIa 4.6400 [2.1560; 9.9858] 15.7 15.8
BEST 2.7500 [1.1582; 6.5297] 12.3 12.7
Number of studies combined: k = 7
HR 95%-CI z p-value
Fixed effect model 2.5928 [1.9141; 3.5122] 6.15 < 0.0001
Random effects model 2.5695 [1.8611; 3.5477] 5.73 < 0.0001
Quantifying heterogeneity:
tau^2 = 0.0181 [0.0000; 0.9384]; tau = 0.1347 [0.0000; 0.9687];
I^2 = 9.4% [0.0%; 73.6%]; H = 1.05 [1.00; 1.94]
Test of heterogeneity:
Q d.f. p-value
6.63 6 0.3569
Details on meta-analytical method:
- Inverse variance method
- DerSimonian-Laird estimator for tau^2
- Jackson method for confidence interval of tau^2 and tau```
CI 输出在所有研究中与原始 CI 匹配到小数点后两位,除了 CARDIa,我认为输入错误(如果我错了请原谅我,但我可以'看不到任何其他解释)。
您可以通过手动计算标准误差然后重新计算置信区间来了解这一点,就像 metagen
函数所做的那样。
library(meta)
se <- meta:::TE.seTE.ci(log(data$LCI), log(data$UCI))$seTE; se
#[1] 0.3823469 0.6027896 0.6762603 0.3004463 0.3735071 0.3910526 0.4412115
data$lower <- round(exp(ci(TE=log(data$HR), seTE=se)$lower), 3)
data$upper <- round(exp(ci(TE=log(data$HR), seTE=se)$upper), 3)
data
studlab HR LCI UCI lower upper
1 Blazek 1.78 0.84 3.76 0.841 3.766 #
2 PRECOMBAT 1.20 0.37 3.93 0.368 3.911 #
3 LE.MANS 1.14 0.30 4.25 0.303 4.291 #
4 NOBLE 2.99 1.66 5.39 1.659 5.388 #
5 MASS-II 2.90 1.39 6.01 1.395 6.030 #
6 CARDIa 4.64 1.33 6.16 2.156 9.986 # <- this one is incorrect.
7 BEST 2.75 1.16 6.54 1.158 6.530 #
CARDIa 的正确 95% CI 应该在 (2.16 - 9.99) 左右。我会验证您输入的值是否正确。
我正在尝试使用风险比、下限和上限 95% 置信区间进行荟萃分析,但例如 CARDIa 研究,获得上限和下限 95%CI ([2.1560; 9.9858]) 是与原始值 (1.33-6.16) 不同,我不知道如何获得确切的数字。
如有任何建议,我们将不胜感激。
使用代码:
library(meta);library(metafor)
data<-read.table(text="studlab HR LCI UCI
Blazek 1.78 0.84 3.76
PRECOMBAT 1.20 0.37 3.93
LE.MANS 1.14 0.3 4.25
NOBLE 2.99 1.66 5.39
MASS-II 2.90 1.39 6.01
CARDIa 4.64 1.33 6.16
BEST 2.75 1.16 6.54
", header=T, sep="")
metagen(log(HR), lower = log(LCI), upper = log(UCI),
studlab = studlab,data=data, sm = "HR")
得到的结果
HR 95%-CI %W(fixed) %W(random)
Blazek 1.7800 [0.8413; 3.7659] 16.4 16.5
PRECOMBAT 1.2000 [0.3682; 3.9109] 6.6 7.1
LE.MANS 1.1400 [0.3029; 4.2908] 5.2 5.7
NOBLE 2.9900 [1.6593; 5.3878] 26.6 25.0
MASS-II 2.9000 [1.3947; 6.0301] 17.2 17.2
CARDIa 4.6400 [2.1560; 9.9858] 15.7 15.8
BEST 2.7500 [1.1582; 6.5297] 12.3 12.7
Number of studies combined: k = 7
HR 95%-CI z p-value
Fixed effect model 2.5928 [1.9141; 3.5122] 6.15 < 0.0001
Random effects model 2.5695 [1.8611; 3.5477] 5.73 < 0.0001
Quantifying heterogeneity:
tau^2 = 0.0181 [0.0000; 0.9384]; tau = 0.1347 [0.0000; 0.9687];
I^2 = 9.4% [0.0%; 73.6%]; H = 1.05 [1.00; 1.94]
Test of heterogeneity:
Q d.f. p-value
6.63 6 0.3569
Details on meta-analytical method:
- Inverse variance method
- DerSimonian-Laird estimator for tau^2
- Jackson method for confidence interval of tau^2 and tau```
CI 输出在所有研究中与原始 CI 匹配到小数点后两位,除了 CARDIa,我认为输入错误(如果我错了请原谅我,但我可以'看不到任何其他解释)。
您可以通过手动计算标准误差然后重新计算置信区间来了解这一点,就像 metagen
函数所做的那样。
library(meta)
se <- meta:::TE.seTE.ci(log(data$LCI), log(data$UCI))$seTE; se
#[1] 0.3823469 0.6027896 0.6762603 0.3004463 0.3735071 0.3910526 0.4412115
data$lower <- round(exp(ci(TE=log(data$HR), seTE=se)$lower), 3)
data$upper <- round(exp(ci(TE=log(data$HR), seTE=se)$upper), 3)
data
studlab HR LCI UCI lower upper
1 Blazek 1.78 0.84 3.76 0.841 3.766 #
2 PRECOMBAT 1.20 0.37 3.93 0.368 3.911 #
3 LE.MANS 1.14 0.30 4.25 0.303 4.291 #
4 NOBLE 2.99 1.66 5.39 1.659 5.388 #
5 MASS-II 2.90 1.39 6.01 1.395 6.030 #
6 CARDIa 4.64 1.33 6.16 2.156 9.986 # <- this one is incorrect.
7 BEST 2.75 1.16 6.54 1.158 6.530 #
CARDIa 的正确 95% CI 应该在 (2.16 - 9.99) 左右。我会验证您输入的值是否正确。