成组序贯设计与消耗函数(Group-Sequential Design)

定义

成组序贯设计:在试验进行中预先指定的若干次「中期分析」时间点上,对累积数据进行假设检验,必要时提前终止。关键问题是:在总体 α 不变的前提下,如何把 α 消耗分配到各次分析?答案是消耗函数 / 边界构造。

核心参数:

  • 总一侧 α(如 0.025)与总 β(如 0.1)。
  • 最大分析次数 k(如 3)。
  • 信息比率 frac = t / t_max(如 c(0.33, 0.67, 1.0))。
  • 消耗函数 sfu(如 O'Brien-Fleming, Pocock, Lan-DeMets with HSD/LDOF/Power)。
  • 非绑定/绑定无效边界(NonBindingFutility)。

资料陈述:O'Brien-Fleming、Pocock、Lan-DeMets 三类边界与消耗函数的定义见 Jennison & Turnbull 专著1;rpact 的实现说明见 Wassmer & Brannath2;PwrGSD 用 Lan-DeMets 消耗函数 + O'Brien-Fleming 型边界实现生存序贯监测3。

主流 R 实现

包 定位 典型入口 独特能力
gsDesign 经典成组序贯 / 两阶段 gsDesign(k=3, sfu=sfHSD, ...) 10+ 消耗函数 (sfHSD/sfLDOF/sfPower/sfTrimmed/sfXG1/2/3…);gsCP/gsPP/gsPosterior/gsSurvCalendar
rpact 适应性 / 确证性序贯 + 非劣等效 getDesignGroupSequential(), getSampleSizeSurvival() 参数互算 getHazardRatioByMedian/getPiByLambda;ASN getPowerAndAverageSampleNumber;仿真家族 getSimulation*
PwrGSD 生存成组序贯 + 加权 log-rank GrpSeqBnds(), wtdlogrank() Fleming-Harrington 权重处理延迟效应;CDFOR2LRR/CRRtoRR 率互算
binseqtest 二叉逐例序贯(逐例入组) designBinseq(), analyzeBinseq() 精确边界、逐例入组监视,与成组序贯互补

关键细节

  • 消耗函数库:gsDesign 自带 21 个 sf* 函数(sfHSD, sfLDOF, sfPower, sfTrimmed, sfXG1/2/3, sfTruncated 等),且可用户自定义。
  • 条件/预测功效:gsDesign 独有 gsCP(条件功效)、gsPP(预测功效)、gsPosterior(后验)。
  • 日历时间换算:gsDesign gsSurvCalendar 把入组/随访日历换算成信息时间。
  • 总样本量 vs. 事件数:生存设计中 gsDesign::nSurvival 返回总 N 而非每组 N;事件数由 nEvents() 给出。gsDesign 的 nNormal/nBinomial 亦返回总 N 而非每组 N。
  • 非 PH 场景:PwrGSD 的 wtdlogrank(formula, WtFun="FH", param=c(0,1)) 实现 Fleming-Harrington 加权 log-rank,适合免疫治疗延迟效应(已跑通:z=−3.87, p≈1e−4)。
  • rpact 独有:参数互算 getHazardRatioByMedian / getPiByLambda、ASN getPowerAndAverageSampleNumber、模拟系列 getSimulation*。
  • binseqtest:二叉逐例序贯边界(逐例入组,非成组),适合罕见病/极早期停止。

关键文献

方法 核心文献 教科书
O'Brien-Fleming O'Brien & Fleming (1979) Jennison & Turnbull (2000)
Pocock Pocock (1977) 同上
Lan-DeMets Lan & DeMets (1983) 同上;Wassmer & Brannath (2025) for rpact
Hwang-Shih-DeCani (HSD) Hwang, Shih & DeCani (1990) 亦见 gsDesign sfHSD
加权 log-rank (FH) Fleming & Harrington (1981) 亦见 PwrGSD 文档

⚠️ 注意:gsDesign 的 nSurvival/nEvents 文档引用 Schoenfeld (1981) Biometrika 68:316–319 与 Lachin & Foulkes (1986),不是 Schoenfeld (1983) Biometrics 39:499–503(后者是另一篇、也常被误引到此处)456。

相关

  • 《R 样本量包:功能重合、同例异构与文献引用比对》(本库私人分析稿,路径见 frontmatter sources) — §1/§2/§4
  • 《R 样本量包:各自的独门功能与实现实例》(本库私人分析稿,路径见 frontmatter sources) — §2/§3/§21
  • 精度导向样本量 — 另一大范式

参考文献

书目核对:上述条目经 Crossref 核对(Jennison & Turnbull 实为 1999 年出版,本页先前写作 2000;Wassmer & Brannath 为 2016 年初版,2025 年为第二版,本页先前的「2025」指第二版)。


  1. Jennison C, Turnbull BW. Group Sequential Methods with Applications to Clinical Trials. Boca Raton: Chapman & Hall/CRC; 1999. DOI: 10.1201/9780367805326. ↩︎

  2. Wassmer G, Brannath W. Group Sequential and Confirmatory Adaptive Designs in Clinical Trials. Cham: Springer; 2016(Springer Series in Pharmaceutical Statistics). DOI: 10.1007/978-3-319-32562-0. ↩︎

  3. PwrGSD 包文档(CRAN)。见 frontmatter source_urls。 ↩︎

  4. Schoenfeld DA. The asymptotic properties of nonparametric tests for comparing survival distributions. Biometrika. 1981;68(1):316–319. DOI: 10.1093/biomet/68.1.316. ↩︎

  5. Schoenfeld DA. Sample-size formula for the proportional-hazards regression model. Biometrics. 1983;39(2):499–503. DOI: 10.2307/2531021. ↩︎

  6. gsDesign 包文档(CRAN,nSurvival/nEvents 的 References 段)。见 frontmatter source_urls。 ↩︎