成组序贯设计与消耗函数(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、ASNgetPowerAndAverageSampleNumber、模拟系列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」指第二版)。
-
Jennison C, Turnbull BW. Group Sequential Methods with Applications to Clinical Trials. Boca Raton: Chapman & Hall/CRC; 1999. DOI: 10.1201/9780367805326. ↩︎
-
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. ↩︎
-
PwrGSD 包文档(CRAN)。见 frontmatter
source_urls。 ↩︎ -
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. ↩︎
-
Schoenfeld DA. Sample-size formula for the proportional-hazards regression model. Biometrics. 1983;39(2):499–503. DOI: 10.2307/2531021. ↩︎
-
gsDesign 包文档(CRAN,
nSurvival/nEvents的 References 段)。见 frontmattersource_urls。 ↩︎