样本量计算
分析目的
样本量计算在研究设计与性能评价方案阶段确定所需样本数或可获得的检验效能。ivdtools
提供三个 独立函数,分别对应三类常见问题:
sample_size_bland_altman():Bland–Altman
一致性评价的样本量/效能(Lu 等 2016);sample_size_proportion_ci():单比例置信区间达到指定半宽所需样本量;sample_size_proportion():单臂目标值检验(单样本比例对已知目标)的样本量/效能/显著性。
三个函数都由参数驱动、不需要原始数据文件。它们返回带属性的标量,print()
给出友好解释,as.numeric()
取出数值结果。示例中的参数(效能、置信水平、允许半宽等)均为教学设定,正式研究
必须使用方案规定的参数。
::: cell
library(ivdtools):::
函数概述
函数 主要用途 求解方向
sample_size_bland_altman()BA 一致性评价 由power求 n,或由n求 power
sample_size_proportion_ci()单比例置信区间 由d求 n,或由n
求d
sample_size_proportion()单臂目标值检验 任意两个给定
n/alpha/power
求第三个
示例一:Bland–Altman 一致性评价的样本量
由目标效能求样本量
::: cell
n_ba <- sample_size_bland_altman( power = 0.80, mu = 0.2, sd = 1, delta = 2.5 ) print(n_ba):::
::: cell
as.numeric(n_ba):::
在平均差异mu = 0.2、差异 SDsd = 1、临床一致性限delta = 2.5
的设定下,达到 80% 效能 约需 201 对样本。
参数说明:
mu为平均差异,sd为差异标准差(>0),delta
为临床可接受的一致性限(>0)。conf.level
是一致性限置信区间的置信水平,agree.level是 LoA
的覆盖度,两者是不同层次------ 前者影响估计的可靠性,后者决定 LoA
本身的宽度。
由样本量求效能
给定已有样本量,反查能达到的效能:
::: cell
sample_size_bland_altman(n = 201, mu = 0.2, sd = 1, delta = 2.5):::
与上面互为验证,n=201 恰好对应效能 0.80。
示例二:单比例置信区间的样本量
由目标半宽求样本量
::: cell
n_pci <- sample_size_proportion_ci(p = 0.3, d = 0.05) print(n_pci):::
对预期比例 p = 0.3,要使 95% 置信区间半宽不超过 0.05,约需 320 例。
由样本量求半宽
::: cell
sample_size_proportion_ci(p = 0.3, n = 320):::
n=320 时半宽约为 0.05,与上面互为验证。
参数说明:
p为预期比例(必须给定,0<p<1),d
为期望半宽(0<d≤0.5)。method可选 7 种置信区间方法(默认"wilson",另有"wald"、"wald-cc"、"agresti-coull"、"jeffreys"、"wilson-cc"、"clopper-pearson"),方法选择影响所需样本量,应在方案中说明。
示例三:单臂目标值检验的样本量
求样本量
::: cell
n_sp <- sample_size_proportion( p0 = 0.2, p1 = 0.35, alpha = 0.05, power = 0.8, alternative = "greater", method = "wilson" ) print(n_sp):::
在 p0 = 0.2(目标/零假设比例)与 p1 = 0.35(预期比例)、单侧 α =
0.05、效能 0.80 的设定下, 约需 47 例。
参数说明:
alternative可取"greater"(此时要求 p1 >
p0)、"less"或"two.sided"。 该方法基于置信区间反演:给定n/alpha/power中任意两个求第三个;三个都给出时作为一致性校验。
求效能
::: cell
sample_size_proportion( p0 = 0.2, p1 = 0.35, n = 47, alpha = 0.05, alternative = "greater" ):::
求显著性水平(alpha)
::: cell
sample_size_proportion( p0 = 0.2, p1 = 0.35, n = 47, power = 0.8, alternative = "greater" ):::
三种求解方向互相印证,结果一致(n=47、power=0.8、alpha≈0.05)。
结果判读要点
- 参数先行:样本量计算前确定效能、置信水平、允许误差/半宽、预期比例与单双侧方向。
- 结果是有条件的:样本量只在所设参数成立时有效,参数估计不准会导致样本量不准。
- 方法一致:比例类问题应提前确定 CI 方法(Wilson
等),不能事后挑选。 - 方向约束:单臂检验注意
alternative与 p0/p1 大小关系的一致性。 - 方案留档:把参数、函数调用与输出一并写入研究方案,便于复核与更新。