平行两组比较样本量
定义与适用范围
最常见的一类设计:受试者按 1:φ 随机分配到标准组 S 与试验组 T,比较单一结局。所有公式共享同一骨架——把「效应/标准差」塞进正态近似式,再乘分配比因子 (1+φ)²/φ:
N_total = [(1+φ)²/φ] × (z_{1−α/2} + z_{1−β})² / (δ_Plan/σ_Plan)² (等分配 φ=1 时因子=4)
第 1 章式 (1.7) 给出最简形态 N = 4(z_{1−α}+z_{1−β})²/Δ²(单侧 z、等分配、连续结局)1;第 5 章式 (5.4) 补上小样本的非中心 t 修正项 +z²_{1−α/2}/22。此后各章只是换「效应」与「方差」的表达:比例差、log OR、log HR、率差、位置移动、AUC 差。
证据与解释
二分类结局(第 3 章)
- χ² 框架(式 3.2):
N = ((1+φ)/φ)·[z√((1+φ)π̄(1−π̄)) + z√(φπ_S(1−π_S)+π_T(1−π_T))]²/δ_Plan²,其中π̄ = (π_S+φπ_T)/(1+φ)(式 3.3)3。 - OR 尺度(式 3.4):
N_OR = ((1+φ)²/φ)·(z+z)²/[(log OR_Plan)²·π̄(1−π̄)]。先由π_T = OR·π_S/(1−π_S+OR·π_S)(式 3.1)把 OR 换成比例3。 - 精确检验框架(式 3.5):先用式 (3.2)/(3.4) 得 n_S,再膨胀
n_S^Exact = n_S/4·[1+√(1+2(1+φ)/(φ·n_S·|δ_Plan|))]²3。 - 已知一个比例(式 3.6):
N = [z√(π_Known(1−π_Known)) + z√(π_T(1−π_T))]²/δ²,单组即可3。
交叉核对(本次实算):式 (3.5) 在等分配 φ=1 时退化为
n/4·[1+√(1+4/(n|δ|))]²,与诊断准确性样本量中 Akoglu 施加的 Yates 连续性校正公式完全相同;用 n=300.88、|δ|=0.05 代入两者都得 339.7。可见「精确检验膨胀」与「Yates 校正」在等分配下是同一件事,只是本书按分析检验(Fisher)命名、Akoglu 按校正来源(Beam)命名。
有序分类结局(第 4 章)
- 结局为 κ 个有序类别时,用比例优势(proportional odds)模型,效应量取单一
log OR,样本量为N = Γ · ((1+φ)²/φ) · (z+z)²/(log OR)²(式 4.3),其中Γ = 3/(1−Σπ̄ᵢ³)(式 4.4)4。 - Γ 是关键:κ=2 时
Γ = 1/(π̄(1−π̄)),式 (4.3) 立刻退回二分类式 (3.4)——即第 3 章是第 4 章的特例;Γ 的最大值为 4(π̄=0.5 时),类别越均匀、类别数越多,Γ 越小、越省样本4。 - 若不能假定比例优势(NPO),则各切点 OR 不同,需另用第 4 章的 NPO 方法(式 4.8 一带)。
- 另一条路线是 Mann-Whitney U 框架,用
P = Pr(Y_T > Y_S)作效应量。
连续型结局(第 5 章)
- 双样本 t(式 5.4):
N = [(1+φ)²/φ]·(z+z)²/(δ_Plan/σ_Plan)² + z²_{1−α/2}/22。 - 分配比的经济性(式 5.5):
N_Unequal = (1+φ)²/(4φ) × N_Equal。φ=1 时因子为 1;φ 偏离 1 越远,总样本量越大(1:2 时约 1.125 倍)2。 - 每组方差不等(式 5.6,Satterthwaite 修正):设 τ=σ_T²/σ_S²,
N_Satt = ((1+φ)/φ){[(τ+φ)σ²/δ²](z+z)² + (τ²+φ³)z²/(2(τ+φ)²)}——方差不等时必须用 τ 重算,不能直接套等方差式2。 - Wilcoxon-Mann-Whitney(式 5.7):在 t 公式上乘效率因子 f。正态假设下
f = π/3 = 1.0472(式 5.8)——即非参数检验要多 4.7% 的样本2。 - 无 ties 的简化式(式 5.9):
N = (1+φ)²(z+z)²/[12φ(P−0.5)²],直接用P = Pr(Y_T > Y_S)表达位置移动2。 - 与已知均值比较:单样本 t(式 5.10/5.11 一带)。
模型推断(本页):书中「最优分配」一节讨论的是 φ 的选取(式 5.5 给出偏离 1:1 的代价),而不是「按标准差比例分配」。后者(Neyman 分配)在本章未见明确表述,本页不当作本书结论。
率结局(第 6 章)
- 泊松(式 6.2):以「每人单位(PU)」计总暴露量
PU = ((1+φ)/φ)(z+z)²(λ_T+φλ_S)/(λ_T−λ_S)²,再除以随访时长换算成人数5。 - 负二项(过度离散)(式 6.4/6.5):
N = (1+φ)·V·(z+z)²/[log(λ_T/λ_S)]²,其中V = (1/τ)(1/λ_S + 1/(φλ_T)) + κ(1+φ)/φ——κ 是过度离散参数,必须预先指定5。 - 上市后监测:单组(无对照、已知背景发生率)、比较性研究、匹配病例对照、多种独立反应各有专门式子(式 6.6–6.12)5。
多组、剂量-反应、与参照组比较、析因(第 14 章,另见第 2 章)
- 无结构多组:先按「最强对比」做两组计算,再按组数放大。
- 剂量-反应:用线性趋势检验,
D = Σ(dᵢ−d̄)²(式 14.3)刻画剂量间距,等间距时最省样本6。 - 与参照组多次比较(Dunnett 型):先用最小有意义差异 δ_Min 算两组所需 N₂(式 14.7),总样本量
N_Total = N₂·√(g−1)(式 14.11)——不是简单相加,√(g−1) 就是多重比较的代价6。 - 析因设计:2×2 等设计可同时回答两个问题,样本量按「最保守的那个对比」定(式 14.x,第 2 章给出框架)67。
多重性与多结局(第 2 章)
- 多个结局或多个比较时,先排序(主/次/探索),或施加 Bonferroni 校正(α/m);重复显著性检验会显著抬高假阳性,必须在设计阶段就决定校正策略7。
- 协变量、组别结构(有序组、回归模型)、流行病学研究设计在第 2 章另有框架;第 2 章还给出一组「方案与发表中如何论证样本量」的实例,以及样本量与伦理的关系7。
例子与推断
模型推断(本页):Γ 与 f 这两个「惩罚因子」值得记住它们的来源:Γ 惩罚的是分类信息被压成二分类时的浪费(把有序结局二分化,等价于把 Γ 从 <4 抬回 4),f 惩罚的是放弃正态假设改用秩检验(+4.7%)。两者都是「分析选择」而非「数据」决定的,因此都可以在规划阶段主动规避。
模型推断(本页):式 (14.11) 的 √(g−1) 放大只在「所有剂量组都与同一参照组比较」时成立;若是「剂量组两两互比」,代价应按比较数而非 √(g−1) 计算。原书此处的推导细节本库未逐行核对,属待核实。
分歧、限制与待核实
| 争议点 | 来源 | 各自说法和范围 | 当前判断 | 下一步核实 |
|---|---|---|---|---|
| 二分类该用 χ² 还是精确检验 | 本书式 (3.2) vs 式 (3.5) | 同一设计两套公式,精确检验版样本量更大 | 不是冲突,是分析计划不同;须与分析方案一致 | 需确认书中对「何时必须用精确检验」的建议原文 |
| Yates/精确膨胀是否该用 | 本书式 (3.5) / 诊断两篇文献 | 本书按 Fisher 精确检验命名;Akoglu 按 Beam 的 Yates 建议施加;Hajian-Tilaki 未施加 | 等分配下公式相同;是否施加仍属惯例分歧(详见诊断页分歧表) | 核对 Beam 1992 原文 |
| Wilcoxon 的效率因子 | 本书式 (5.8) | 正态假设下 f=π/3=1.0472 | 采纳(有明确推导) | 非正态分布下 f 的取值需查原书表 |
限制:本书全部公式建立在正态近似上,小样本或极端比例(接近 0/1)时需用精确方法;第 3 章明确提示「某一比例接近 0 或 1 时需注意困难」3。
相关知识
- 精度导向样本量 — 同一骨架的「CI 宽度」版本(本书第 9 章)
- 非劣效性与等效性样本量 — 单侧假设下的同一骨架
- 配对与重复测量样本量 — 观测不独立时的修正
- 整群随机试验样本量 — 观测成簇时的设计效应
- 样本量手册导航 — 全书章节地图
参考文献
-
Machin D, Campbell MJ, Tan SB, Tan SH. Sample Sizes for Clinical, Laboratory and Epidemiology Studies. 4th ed. Wiley-Blackwell; 2018. DOI: 10.1002/9781118874905. ISBN 9781118874943 / 9781118874905. 依据章节:第 1 章「基本设计考量」(式 1.7、1.8–1.10)。 ↩︎