第3课 统计推断基础
课程主题: 理解 P 值之外的信息
课程定位: 从"样本为什么不能直接等同于总体"出发,建立抽样误差、参数估计、置信区间、假设检验、Ⅰ/Ⅱ类错误、效应量和多重比较之间的完整逻辑。\
一、学习目标
完成本节后,应能够:
- 区分总体、样本、总体参数和样本统计量;
- 理解为什么同一总体反复抽样会得到不同结果,以及什么是抽样误差;
- 区分 SD(标准差)与 SE(标准误);
- 理解点估计与置信区间分别提供什么信息;
- 正确解释 95% CI,避免"95%概率包含真值"等常见误解;
- 理解假设检验中 H₀、H₁、检验统计量、P 值和 α 的关系;
- 正确解释 P<0.05 与 P>0.05,尤其避免把 P>0.05 解释成"证明两组相同";
- 理解Ⅰ类错误、Ⅱ类错误和统计功效(power);
- 理解为什么不能只报告 P 值,而应结合效应量和 95% CI;
- 理解多重比较为什么会增加假阳性风险,并识别论文中的常见问题。
3.1 从描述统计到统计推断:我们真正想知道的通常是"总体"
上一课解决的是:
手中的数据是什么?应该怎样描述?
例如某临床试验纳入 100 名患者,平均收缩压为:
138.6 ± 15.2 mmHg
这首先是在描述这 100 名患者。
但医学研究通常不会停在这里。研究者真正关心的往往是:
这 100 名患者所代表的更大患者总体,其平均收缩压大约是多少?
或者:
在研究样本中观察到的治疗组与对照组差异,能否推断到目标总体?
这就从描述性统计(descriptive statistics)进入了统计推断(statistical inference)。
统计推断的核心问题可以概括为:
总体通常无法全部观察
↓
从总体中抽取样本
↓
观察和分析样本
↓
利用样本信息推断总体
↓
但样本存在随机波动
↓
因此推断必须表达“不确定性”
3.2 四个最基础的概念:总体、样本、参数、统计量
3.2.1 总体 Population
总体是研究希望推断的全部研究对象。
例如:
希望研究某种降压药对目标高血压患者群体的降压效果。
这里真正关心的是符合研究定义的目标患者总体,而不是仅仅已经入组的 200 名患者。
3.2.2 样本 Sample
现实中通常无法测量总体中的每一个对象,因此研究者从总体中抽取一部分研究对象进行观察。
例如:
目标总体:符合条件的高血压患者
↓
抽取 200 人
↓
研究样本:200 名患者
3.2.3 总体参数 Parameter
描述总体特征的数值称为参数,例如:
- 总体均值 μ;
- 总体标准差 σ;
- 总体比例 p;
- 总体中的真实均值差;
- 总体中的真实 OR、RR、HR 等。
这些才是研究者通常真正希望了解的量。
3.2.4 样本统计量 Statistic
从样本数据计算出来的数值称为样本统计量,例如:
- 样本均值 x̄;
- 样本标准差 s;
- 样本比例;
- 样本中的均值差;
- 样本估计得到的 OR、RR、HR 等。
最重要的关系是:
样本统计量是我们观察到的;总体参数通常是我们希望估计的。
3.3 为什么样本结果每次都可能不同:抽样误差
假设某个总体的真实平均收缩压为:
μ = 140 mmHg
现在每次随机抽取 100 名患者。
第一次:
样本均值 = 138.2 mmHg
第二次:
样本均值 = 141.0 mmHg
第三次:
样本均值 = 139.3 mmHg
三个样本都来自同一个总体,但样本均值并不完全相同。
原因并不是"统计软件算错了",而是:
不同随机样本包含的个体不同,因此样本统计量会发生随机波动。
这种由随机抽样导致的样本统计量与总体参数之间、或不同样本统计量之间的随机差异,就是理解**抽样误差(sampling error)**的起点。
需要区分:
抽样误差 ≠ 数据录入错误
抽样误差 ≠ 测量仪器错误
抽样误差 ≠ 研究设计错误
它是抽样研究天然存在的不确定性。
3.4 从抽样误差到标准误:我们需要描述"估计有多不稳定"
假设我们不是只抽一次样,而是从同一个总体反复抽取很多个同样大小的样本。
每个样本都会产生一个样本均值:
样本1 → 138.2
样本2 → 141.0
样本3 → 139.3
样本4 → 140.5
……
这些"样本均值"本身也形成一个分布。
这就是理解**抽样分布(sampling distribution)**的关键。
而标准误(standard error, SE)描述的正是:
样本统计量在重复抽样中的变异程度。
对于样本均值,一个常见形式为:
其中:
- s:样本标准差;
- n:样本量。
因此,在其他条件相同时:
n 增大
↓
SE 减小
↓
估计通常更精确
3.5 SD 和 SE 必须分清
这是医学论文中非常重要、也非常容易混淆的一组概念。
| 指标 | SD | SE |
|---|---|---|
| 中文 | 标准差 | 标准误 |
| 主要描述 | 个体数据的离散程度 | 统计量估计的不确定性/精度 |
| 回答的问题 | 患者之间差异有多大? | 这个样本估计有多稳定? |
| 与样本量关系 | 不会简单因 n 增大而趋近 0 | 通常随 n 增大而减小 |
| 常见用途 | 描述原始数据 | 构造 CI、统计推断 |
例如:
Age = 56.2 ± 12.0 years
如果这里的 12.0 是 SD,说明患者年龄之间存在一定离散程度。
如果作者把结果写成:
Age = 56.2 ± 1.2 years
而 1.2 实际上是 SE,则它表达的是均值估计的精度,不是患者年龄本身的离散程度。
因此:
SD 和 SE 不能因为都写成"Mean ± 某个数"就认为含义相同。
在描述研究对象的原始连续型数据时,尤其是 Table 1,应仔细核实作者报告的是 SD 还是 SE。
3.6 参数估计:点估计只能给一个数,区间估计才能表达不确定性
3.6.1 点估计 Point estimate
假设治疗组平均收缩压下降:
10.2 mmHg
对照组平均下降:
5.2 mmHg
两组样本均值差为:
10.2 − 5.2 = 5.0 mmHg
这个 5.0 mmHg 是对总体真实治疗差异的一个点估计。
但只报告:
Mean difference = 5.0 mmHg
还不能告诉我们这个估计有多精确。
例如下面两项研究的点估计完全相同:
研究A:Mean difference = 5.0 mmHg
研究B:Mean difference = 5.0 mmHg
但如果:
研究A:95% CI 4.2–5.8
研究B:95% CI -1.0–11.0
两项研究提供的信息显然完全不同。
3.7 置信区间:不仅看"有没有效",还要看"可能有多大"
置信区间(confidence interval, CI)给出与数据和模型相容的一组参数值,并反映估计的精度。
例如:
Mean difference = 5.0 mmHg
95% CI: 1.0–9.0 mmHg
这里至少提供了两类信息:
点估计:5.0 mmHg
区间估计:1.0–9.0 mmHg
区间较窄通常意味着估计较精确;区间较宽意味着不确定性较大。
3.7.1 95% CI 应该怎样理解
在经典频率学派框架下,更严谨的理解是:
如果在相同条件下反复抽样,并每次使用同样的方法构造 95% CI,那么长期来看,约 95% 的这些区间会覆盖真实总体参数。
不建议把已经计算出来的某一个固定区间直接解释为:
"真实值有 95% 的概率落在这个区间里。"
因为在频率学派框架下,总体参数被视为固定值,随机的是抽样及由此构造出的区间。
培训和稿件审核中,更实用的理解是:
95% CI 同时帮助我们判断效应方向、估计精度以及与哪些真实效应大小相容。
3.8 先认识"无效值":不同效应指标的判断基准不同
置信区间能否跨过"无效值",对结果解释非常重要。
差值类指标
例如:
Mean difference
Risk difference
无效值通常是:
0
因为"差值 = 0"意味着两组没有差异。
例如:
Mean difference = 5.0
95% CI: 1.0–9.0
CI 不跨 0。
而:
Mean difference = 5.0
95% CI: -1.0–11.0
CI 跨 0。
比值类指标
例如:
OR
RR
HR
无效值通常是:
1
因为比值 = 1 表示两组相等。
例如原课件在 Logistic 回归质量评价中强调:
OR = 2.31
95% CI: 1.42–3.76
OR 位于其 CI 内,且 CI 不跨 1。
如果:
OR = 2.31
95% CI: 0.82–5.64
虽然点估计 OR > 1,但 CI 跨 1,不能仅根据 OR 的点估计声称存在明确的风险增加证据。
3.9 为什么还需要假设检验
参数估计回答的是:
效应大约是多少?不确定性有多大?
假设检验(hypothesis testing)则从另一个角度提出问题:
如果总体中实际上不存在研究者所关注的差异或关联,我们现在观察到这样的数据是否罕见?
这就需要建立两个假设。
3.10 原假设 H₀ 与备择假设 H₁
假设研究一种新降压药。
原假设 Null hypothesis, H₀
例如:
H₀:治疗组与对照组总体平均降压幅度无差异
可以写成:
备择假设 Alternative hypothesis, H₁
双侧检验时:
H₁:两组总体平均降压幅度存在差异
即:
接下来根据研究设计和数据类型选择相应的统计检验,计算检验统计量和 P 值。
3.11 P 值究竟是什么
P 值最容易被误解。
在给定统计模型和原假设成立的前提下,P 值表示:
观察到当前检验统计量这样极端或更加极端结果的概率。
可以把逻辑理解成:
先假设 H₀ 成立
↓
问:如果真的没有差异
↓
现在这样的数据有多“不寻常”?
↓
用 P 值量化
例如:
P = 0.03
不能解释成:
“原假设为真的概率只有 3%”
也不能解释成:
“有 97% 的概率治疗有效”
更不能解释成:
“结果由偶然造成的概率是 3%”
这些都不是经典 P 值的定义。
3.12 α = 0.05 与 P<0.05 到底意味着什么
研究者通常会事先设定显著性水平 α,例如:
α = 0.05
然后比较:
P 与 α
在传统假设检验框架中:
P < 0.05
→ 在 α=0.05 的规则下拒绝 H₀
→ 通常表述为“差异有统计学意义”
如果:
P ≥ 0.05
→ 在该显著性水平下不拒绝 H₀
注意这里应说:
不拒绝 H₀
而不是:
接受 H₀ / 证明 H₀ 为真。
3.13 P>0.05 不等于"两组相同"
假设某项小样本研究得到:
Mean difference = 5.0 mmHg
95% CI: -2.0–12.0 mmHg
P = 0.15
正确的信息是:
当前数据在设定的检验框架下没有提供足够证据拒绝"无差异"的原假设;同时 CI 很宽,提示估计存在较大不确定性。
不能直接写成:
"两种治疗效果相同。"
为什么?
因为 CI 同时容许:
轻微反方向差异
无差异
具有一定临床意义的正向差异
所以:
P>0.05
≠ 证明没有差异
≠ 证明两组等效
真正要证明"等效"或"非劣",需要相应的研究设计、预设界值和统计方法。
3.14 P<0.05 也不等于"差异很大"
反过来同样成立。
假设一项非常大的研究:
n = 20,000
Mean difference = 0.5 mmHg
95% CI: 0.2–0.8 mmHg
P < 0.001
P 值很小。
但临床问题应该继续问:
0.5 mmHg 的差异是否具有实际临床意义?
因此:
统计学显著
≠ 效应一定很大
≠ 临床意义一定重要
这就是为什么不能只看 P 值。
3.15 为什么统计决策可能犯错
真实世界中我们并不知道 H₀ 究竟是真是假,只能根据样本作出统计决策。
因此存在四种情况:
| 真实情况 | 不拒绝 H₀ | 拒绝 H₀ |
|---|---|---|
| H₀ 为真 | 正确决策 | Ⅰ类错误 |
| H₀ 为假 | Ⅱ类错误 | 正确发现 |
3.16 Ⅰ类错误:假阳性
Ⅰ类错误(Type I error)指:
H₀ 实际为真,却被错误拒绝。
直观理解:
实际上没有差异
↓
研究结果却判断“有差异”
↓
假阳性
Ⅰ类错误概率通常记为:
例如设:
α = 0.05
就是在相应检验框架和假设下,将Ⅰ类错误率控制在预设水平。
原课件在多组数据反复使用 t 检验的问题中明确提醒:多组数据进行不恰当的重复两两 t 检验,会增加Ⅰ类错误,即假阳性风险。
3.17 Ⅱ类错误:假阴性
Ⅱ类错误(Type II error)指:
H₀ 实际为假,但研究未能拒绝 H₀。
直观理解:
实际上存在差异
↓
研究却没有检测出来
↓
假阴性
Ⅱ类错误概率记为:
例如真实存在有意义的治疗差异,但样本量太小、数据变异过大,研究可能得到:
P > 0.05
这并不自动说明真实差异不存在。
3.18 Power:发现真实效应的能力
统计功效(power)通常定义为:
它描述:
当预设的真实效应确实存在时,统计检验正确拒绝 H₀ 的概率。
例如:
Power = 80%
意味着在设计所设定的效应大小、变异程度、显著性水平和其他条件下,研究有相应的概率检测到该效应。
Power 与多种因素有关,包括:
- 样本量;
- 效应大小;
- 数据变异;
- α 水平;
- 研究设计和统计方法。
因此样本量设计并不是简单地追求"人数越多越好",而是要结合预期效应、变异程度、α 和目标 power 等因素进行规划。
3.19 P 值之外为什么还要看效应量
P 值主要告诉我们:
数据与 H₀ 的不相容程度如何?
但临床研究还需要回答:
差异究竟有多大?
这就是效应量(effect size)的价值。
原课件在比较常用检验方法时列出了相应效应量:
| 分析类型 | 效应量/效应指标 |
|---|---|
| t 检验 | Cohen's d;原始均值差 |
| 单因素 ANOVA | η²;组均值/两两均值差 |
| 秩和检验 | r;中位数差/Hodges–Lehmann 估计 |
| 卡方检验 | Cramer's V / φ;率差、OR、RR |
因此完整结果不应该只剩:
P = 0.003
而应该尽可能回答三个层次的问题:
1. 效应方向是什么?
2. 效应有多大?
3. 这个估计有多精确?
即:
效应估计 + 95% CI + P值
3.20 原始效应量与标准化效应量
为了理解课件中的效应量,可以进一步区分两类。
3.20.1 原始尺度效应
例如:
Mean difference = -5.0 mmHg
Risk difference = -8%
OR = 2.3
RR = 1.8
它们通常具有较直接的医学解释。
例如:
平均收缩压降低 5 mmHg
比单独报告 Cohen's d 更容易直接联系临床意义。
3.20.2 标准化效应量
例如:
Cohen's d
η²
r
Cramer's V
φ
这些指标将效应转换到某种标准化尺度,便于描述效应强弱或跨尺度比较。
但不能把所有研究机械地只用一个标准化效应量。
应根据:
研究问题
+
数据类型
+
统计方法
+
医学解释需求
选择合适指标。
3.21 一个案例:为什么必须同时看 P、效应量和 CI
比较两个研究。
研究 A
n = 40
Mean difference = 5.0 mmHg
95% CI: -0.6–10.6
P = 0.08
研究 B
n = 10,000
Mean difference = 0.5 mmHg
95% CI: 0.2–0.8
P < 0.001
如果只看 P:
A → “不显著”
B → “显著”
但进一步看:
研究 A 的点估计较大,但 CI 较宽,说明数据仍然存在较大的不确定性。
研究 B 的估计非常精确,统计学证据较强,但实际差异只有 0.5 mmHg,是否具有临床重要性需要结合专业背景判断。
因此:
P 值回答的不是"效应有多大",CI 和效应量提供了 P 值无法替代的信息。
3.22 P 值与 95% CI 为什么经常可以相互核对
对于常见的双侧检验,在模型、显著性水平和 CI 构造方法相互对应时,P 值和 95% CI 通常具有一致关系。
差值类指标
95% CI 不跨 0
↔ 通常 P < 0.05
比值类指标
OR / RR / HR 的 95% CI 不跨 1
↔ 通常 P < 0.05
原课件在 Logistic 回归审核中就使用了这一原则:
OR = 2.565
95% CI: 1.382–4.761
P = 0.003
CI 不跨 1,P<0.05,方向一致。
而如果稿件出现:
OR = 2.31
95% CI: 1.42–3.76
P = 0.28
则应重点核查。
但编辑审核时不要把这个规律机械绝对化,因为在不同检验方法、CI 构造方法、单/双侧检验、舍入等情况下可能存在差异。
3.23 为什么多重比较会增加假阳性
假设研究者只做一次检验,并设:
α = 0.05
现在假设研究者进行很多次独立检验,每次都仍使用 0.05 作为阈值。
即使所有 H₀ 都是真的,只要检验次数不断增加:
至少出现一次假阳性的机会就会增加。
在"20 次相互独立且每次 α=0.05"的理想化例子中:
约为:
64%
这个例子不是说所有现实多重比较都恰好是 64%,而是用于说明:
重复进行大量未经校正的检验,会导致整体Ⅰ类错误风险膨胀。
3.24 一个医学研究中的典型多重比较错误
假设研究有 4 个治疗组:
A
B
C
D
作者分别进行:
A vs B
A vs C
A vs D
B vs C
B vs D
C vs D
一共 6 次两两比较。
如果每一次都直接按照:
P < 0.05
判断,而完全不考虑整体多重检验问题,假阳性风险可能被放大。
这正对应原课件强调的问题:
多组数据不能简单通过反复两两 t 检验解决;应先采用适合多组比较的总体方法,并在需要进一步定位差异时采用适当的事后多重比较策略。
3.25 多重比较的基本处理思路
本节只建立概念,不深入具体统计方法。
常见思路包括:
Bonferroni
如果计划进行 m 次比较,可将单次比较阈值调整为:
例如:
α = 0.05
m = 5
则:
α* = 0.01
Bonferroni 简单、保守,主要用于控制 family-wise error rate(FWER)。
Holm
Holm 方法也是控制 FWER 的常用方法,通常比简单 Bonferroni 更有统计效率。
FDR
在基因组学、转录组学等一次进行大量假设检验的场景中,经常更关注 false discovery rate(FDR)。
本节不要求记忆各种算法,而要记住:
只要同一研究问题中进行了大量检验,就应该主动想到 multiplicity。
3.26 "总体检验显著"以后为什么还要事后比较
假设有 4 个治疗组。
ANOVA 总体检验得到:
P < 0.001
这首先说明:
至少存在某些组之间的总体均值差异证据。
但它没有直接告诉我们:
A vs B?
A vs C?
A vs D?
B vs C?
B vs D?
C vs D?
因此需要根据研究设计和分析方案进行适当的事后比较(post hoc comparisons)。
原课件在方差分析部分明确采用了:
主效应
↓
若需要明确具体哪些水平不同
↓
事后多重比较
所以:
总体检验解决"有没有总体差异",事后比较进一步解决"差异具体在哪里"。
3.27 把整套统计推断串起来
现在可以把本课所有概念连成一条逻辑链:
研究关注总体
↓
但只能观察样本
↓
随机抽样产生抽样误差
↓
样本统计量存在不确定性
↓
SE 描述统计量的抽样变异
↓
点估计给出“最核心的估计值”
↓
CI 表达估计的不确定性
↓
假设检验评估数据与 H₀ 的相容程度
↓
P 值提供相应统计证据
↓
统计决策存在Ⅰ类/Ⅱ类错误
↓
Power 描述检测预设真实效应的能力
↓
效应量回答“效应有多大”
↓
CI回答“估计有多精确”
↓
多重检验时必须控制错误率膨胀
3.28 医学论文中一个结果应该怎样完整阅读
假设论文报告:
Treatment group: -10.2 ± 8.1 mmHg
Control group: -5.2 ± 7.9 mmHg
Mean difference = -5.0 mmHg
95% CI: -7.8 to -2.2 mmHg
P = 0.001
不要只看最后一个:
P = 0.001
应该依次阅读:
① 效应方向
Mean difference = -5.0 mmHg
首先明确差异方向及组别定义。
② 效应大小
两组平均变化相差约:
5 mmHg
然后判断这个大小是否具有医学意义。
③ 估计精度
95% CI: -7.8 to -2.2
CI 给出了与数据相容的效应范围,并反映估计精度。
④ 统计学证据
P = 0.001
在相应 H₀ 和模型假设下,数据对"无差异"假设提供较强的反对证据。
因此完整解释的顺序更接近:
效应方向
→ 效应大小
→ 95% CI
→ P值
→ 临床意义
而不是:
先找 P<0.05
→ 找到了
→ 分析结束
3.29 医学论文编辑审核:P 值之外应该检查什么
第一层:描述值和推断值有没有混淆
例如:
Mean ± SD
主要描述个体数据分布。
而:
Mean difference + 95% CI
用于表达组间效应估计及其不确定性。
检查作者是否把:
SD
SE
CI
混为一谈。
第二层:P 值解释是否正确
重点警惕:
P=0.03 → “H₀为真的概率是3%” ×
P=0.03 → “治疗有效概率为97%” ×
P>0.05 → “两组完全相同” ×
P<0.05 → “具有重要临床意义” ×
第三层:P 值和 CI 是否基本一致
对于常见双侧检验:
差值类 CI 是否跨 0?
OR/RR/HR 的 CI 是否跨 1?
P 值是否与之基本一致?
如果明显冲突,应进一步核查。
第四层:点估计是否落在 CI 内
例如:
OR = 2.31
95% CI: 1.42–3.76
合理。
如果:
OR = 2.31
95% CI: 0.82–1.95
点估计本身甚至不在所报告的 CI 内,属于明显需要核查的问题。原课件也将这一点列为 Logistic 回归结果审核重点。
第五层:是否只报告 P 值而没有效应信息
例如只写:
P < 0.001
却不提供:
均值差
OR / RR / HR
相关系数
效应量
95% CI
读者就很难判断效应大小和估计精度。
第六层:多组比较有没有反复做未经校正的两两检验
尤其看到:
3组、4组、5组……
却大量使用:
两两 t 检验
应立即想到:
Ⅰ类错误膨胀
+
多重比较问题
3.30 本课最常见的八个错误认识
| 错误认识 | 正确理解 |
|---|---|
| P<0.05 就说明效应很大 | P 值不直接表示效应大小 |
| P>0.05 就说明两组相同 | 只能说当前分析未提供足够证据拒绝 H₀ |
| P=0.03 表示 H₀ 为真的概率是 3% | P 值不是 H₀ 的后验概率 |
| 95% CI 表示真值有 95% 概率在这个已计算区间里 | 频率学派 CI 应从重复抽样覆盖率理解 |
| SD 和 SE 都是“离散程度”,可以互换 | SD 描述个体变异,SE 描述统计量估计的不确定性 |
| 样本越大,效应一定越大 | 样本量影响精度和检验能力,不决定真实效应大小 |
| 做很多次 P<0.05 检验没有问题 | 多重检验会增加整体假阳性风险 |
| 总体检验显著就说明每两组之间都显著 | 总体显著不等于所有两两比较均显著 |
3.31 最终判断框架
看到论文中的统计结果时,可以按下面的顺序思考:
这项研究想推断什么总体参数?
↓
样本给出的点估计是什么?
↓
效应方向和大小是什么?
↓
95% CI 有多宽?
是否跨无效值?
↓
P 值是多少?
与 CI 是否基本一致?
↓
是否可能存在Ⅰ类/Ⅱ类错误问题?
↓
样本量和 power 是否需要关注?
↓
是否报告了合适的效应量?
↓
是否进行了多次检验?
是否考虑 multiplicity?
↓
最后再综合判断统计意义和医学意义
3.32 本节小结
统计推断的核心并不是记住:
P < 0.05
而是理解:
我们只观察了样本,却希望认识总体,因此所有统计推断都伴随着不确定性。
这条主线可以浓缩为:
总体
↓
样本
↓
抽样误差
↓
SE
↓
点估计 + 置信区间
↓
假设检验 + P值
↓
Ⅰ类错误 / Ⅱ类错误 / Power
↓
效应量
↓
多重比较控制
因此,阅读医学论文结果时,最终应形成一个稳定习惯:
不要只问"P 是否小于 0.05",还要问:效应有多大?方向是什么?95% CI 有多宽?结果有多精确?是否具有临床意义?是否存在多重检验和错误率问题?
下一课进入具体统计方法以后,t 检验、ANOVA、秩和检验、卡方检验、相关分析和回归分析中的 P 值、CI 和效应量,就可以放回这套统一框架中理解。