跳到主内容

四种饮食方式对体重的影响

研究背景

营养师想比较四种流行饮食方式(低碳水、地中海、生酮、均衡饮食)对3个月体重变化的影响。招募了60名志愿者,随机分配到四组(每组15人)。数据严重右偏(少数人减重很多,多数人变化不大)。

研究问题

四种饮食方式的体重变化是否存在差异?哪种方式效果最好?

数据说明

四个独立组别的体重变化数据(kg)。数据不满足正态性,改用Kruskal-Wallis检验。

示例数据

  • 低碳水:-2, -3, -1, -4, -2, -5, -3, -2, -4, -1, -3, -6, -2, -4, -3
  • 地中海:-1, -2, 0, -1, -2, -1, -3, 0, -2, -1, -2, -1, -3, 0, -2
  • 生酮:-4, -5, -3, -6, -4, -7, -5, -4, -6, -3, -5, -8, -4, -6, -5
  • 均衡:0, -1, 1, 0, -1, 0, -2, 1, -1, 0, -1, 0, -2, 1, -1

分析步骤

  1. 步骤 1:为什么用Kruskal-Wallis而不是ANOVA

    ANOVA要求每组数据近似正态且方差齐性。本案例数据严重右偏,Shapiro-Wilk检验各组p<0.05。此时应使用Kruskal-Wallis检验(非参数方法,相当于单因素ANOVA的非参数版本)。

    提示:Kruskal-Wallis检验比较的是四组的秩次分布,而非原始数值。

  2. 步骤 2:解读H统计量

    H=38.72, df=3, p<0.001。说明至少有两组的体重变化分布存在显著差异,但不知道是哪两组。

  3. 步骤 3:事后两两比较(Dunn检验)

    需要进行事后检验(通常用Dunn检验+Bonferroni校正)。结果显示:生酮 vs 均衡(p<0.001)、低碳水 vs 均衡(p=0.002)、生酮 vs 地中海(p=0.01),其他比较不显著。

结果解读

结论

Kruskal-Wallis检验显示,四种饮食方式的体重变化差异具有统计学意义,H(3)=38.72, p<0.001。事后Dunn检验表明,生酮饮食(Mdn=-5.0, IQR=-6 to -4)减重效果最好,显著优于均衡饮食(Mdn=0, IQR=-1 to 1)和地中海饮食(Mdn=-1.0, IQR=-2 to 0)。

为什么这样解读?

Kruskal-Wallis的H检验显著只说明"不全相同",必须做事后检验才能知道具体差异模式。报告时应关注中位数和IQR,而非均值。效应量可用η²=H/(N-1),这里η²=0.66,说明饮食方式解释了66%的体重变化变异。

如何在论文中报告

采用Kruskal-Wallis检验比较四种饮食方式的体重变化效果(α=0.05)。结果显示,饮食方式对体重变化存在统计学显著关联,H(3)=38.72, p<0.001, η²=0.66。事后Dunn检验(Bonferroni校正)表明,生酮饮食(Mdn=-5.0 kg, IQR=-6 to -4)减重效果显著优于均衡饮食(Mdn=0 kg, IQR=-1 to 1, p<0.001)和地中海饮食(Mdn=-1.0 kg, IQR=-2 to 0, p=0.01)。这表明短期内生酮饮食的减重效果最佳。

常见错误

Kruskal-Wallis显著后直接说"四组都不同"

H检验显著只说明"至少有两组不同",必须做事后两两比较(Dunn检验)才能确定具体差异模式。

报告均值和标准差

非参数检验应报告中位数和四分位距(IQR)。均值受极端值影响大,不能代表偏态数据的"典型值"。

忽略事后检验的多重比较校正

做6次两两比较(4组C(4,2)=6),每次α=0.05,总体犯第一类错误的概率会膨胀。必须用Bonferroni校正(α/6=0.0083)或其他方法控制整体α。

动手跑一遍

用本案例数据打开 Kruskal-Wallis H 计算器,亲手复现统计量、p 值与置信区间,数据只在浏览器内处理。

打开 Kruskal-Wallis H →