You need to enable JavaScript to run this app.
优惠活动
大模型
产品
解决方案
定价
更多

脑脊液分流术限制性平均生存时间(RMST)Meta分析方法学正确性咨询

脑脊液分流术RMST单均数Meta分析方法学疑问

研究背景与操作步骤

针对脑积水的脑脊液分流术开展**限制性平均生存时间(RMST)**单均数Meta分析,操作流程如下:

  • 对已发表的生存曲线进行数字化处理
  • 借助R包IPDfromKM提取个体患者数据
  • 重建生存曲线并计算12个月RMST

已获取的研究数据

study <- c("study1", "study2", "study3", "study4", "study5")
n_patients <- c(535, 209, 111, 599, 434)
rmst_12 <- c(10.54759, 11.36175, 10.50244, 10.51183, 8.716552)
se_12 <- c(0.1463532, 0.1439506, 0.3246873, 0.1471398, 0.2374582)

# 数据框预览
> data
   study n_patients  rmst_12     se_12
1 study1        535 10.54759 0.1463532
2 study2        209 11.36175 0.1439506
3 study3        111 10.50244 0.3246873
4 study4        599 10.51183 0.1471398
5 study5        434 8.716552 0.2374582

Meta分析实施代码与结果

# 从标准误(SE)计算标准差(SD)
data$sd_12 <- data$se_12 * sqrt(data$n_patients)

# 加载meta包
require(meta)

# 用metamean函数开展Meta分析
mm_12 <- metamean(n = n_shunts,
                  mean = rmst_12,
                  sd = sd_12,
                  studlab = author_year, 
                  data = data, 
                  method.mean = "Luo",
                  method.sd = "Shi",
                  sm = 'MRAW',             
                  random = TRUE,          
                  warn = TRUE,
                  prediction = TRUE,
                  method.tau = "REML")

# 分析结果
> mm_12
Number of studies combined: k = 5
Number of observations: o = 1888

                        mean             95%-CI
Common effect model  10.5757 [10.4247; 10.7268]
Random effects model 10.3373 [ 9.4868; 11.1878]
Prediction interval          [ 7.0212; 13.6534]

Quantifying heterogeneity:
 tau^2 = 0.8975 [0.2937; 7.7528]; tau = 0.9474 [0.5419; 2.7844]
 I^2 = 95.6% [92.3%; 97.5%]; H = 4.78 [3.61; 6.34]

Test of heterogeneity:
     Q d.f.  p-value
 91.39    4 < 0.0001

Details on meta-analytical method:
- Inverse variance method
- Restricted maximum-likelihood estimator for tau^2
- Q-Profile method for confidence interval of tau^2 and tau
- Prediction interval based on t-distribution (df = 3)
- Untransformed (raw) means

疑问

分析过程无报错,但不确定上述操作在方法学层面是否正确,恳请专业人士给予指导。

内容的提问来源于stack exchange,提问作者CharlesLDN

相关产品推荐
方舟 Agent Plan

超全模态模型 × Harness 升级,最新支持 Deepseek-V4.1-Flash、GLM-5.3 系列、Doubao-Seedream-5.0-pro、Kimi-K3 (部分), 限时 9.9 元起

最近更新时间:2026.07.25 10:07:41