强烈不推荐重复测量方差分析?

临床研究,特别是随机对照研究,凡是碰到重复测量资料者,似乎大家都喜欢用重复测量方差分析。

它是一种被严重误用和滥用的方法!
 

第一,如果你是前后重复测量两次者,必定不能用重复测量方差分析。关于这一点,我在之前的推文已经介绍过了。

还没有王法?临床试验干预前后测量2次专家说要用重复测量方差分析?

第二,如果是测量多次呢?重复测量方差分析带来的是平均效应,平均效益很多时候不能代表者我们想要探讨的效应(特别是存在着交互效应时候)。

举个例子:

•一项临床试验,评价干预措施对新冠肺炎患者疗效。分别在治疗前、治疗后2周、4周、6周进行随访调查。主要结局指标EQ-5D生命质量。我计算治疗后2周、4周、6周与治疗后的差值。画成下面的图。

如果通过重复测量方差分析,得到的结果是2周、4周、6周效应的的平均值

这个3.94,不代表2周,也不代表4周和6周,不是好结果。

第三,重复测量方差分析不会因为重复测量资料增加而检验效能增加。

第四,重复测量方差分析无法克服缺失值的问题,而绝大多数临床试验,特别随访多次时,都存在缺失值的问题。

第五,重复测量方差分析无法克服结局偏态数据、等级数据的问题。

第六,一旦用了重复测量方差分析,意味着的临床试验的结局指标设定和评价定义是不清晰的,这不是好的临床试验!

第七、很多说应用重复测量方差分析因为关注时间效应!我想说的是,随机对照研究的时间效应,没有太大意义。

第八,临床试验更推荐协方差分析、广义估计方程或者混合线性模型、广义混合线性模型来进行分析。

早在2017年,Gosho and Sato 在J Clin Pharm Ther上发表了一篇题为<Trends in study design and the statistical methods employed in a leading general medicine journal>的文章,更加详细系统地将2015年全年发表在NEJM上的238项原始研究重新分类并汇总。包括研究设计类型(Table 1),根据干预性研究vs观察性研究分别汇总统计方法应

在“repeated-measures analysis(重复测量分析)”,最常用的是mixed-effects model混合效应模型(11%),广义估计方程(7%),广义线性混合模型和重复测量方差分析(各3%),重复测量资料已经很少用了!

### 如何在 Python 中进行重复测量方差分析 对于希望执行重复测量方差分析的研究者而言,在 Python 生态系统中有多种库可供选择。`pingouin` 是一个专注于提供易于使用的统计测试函数的库,非常适合处理此类问题[^3]。 #### 使用 `pingouin` 进行重复测量方差分析 安装所需库可以通过 pip 完成: ```bash pip install pingouin pandas numpy statsmodels seaborn matplotlib scipy ``` 下面是一个简单的例子展示如何利用 `pingouin` 来完成重复测量 ANOVA 的操作: ```python import numpy as np import pandas as pd from pingouin import rm_anova, read_dataset # 加载内置的数据集作为示例 data = read_dataset('rm_anova') print(data.head()) # 执行单因素重复测量方差分析 result_single_factor = rm_anova(dv='DesireToKill', within='Disgustingness', subject='Subject', data=data) print(result_single_factor[['Source', 'ddof1', 'ddof2', 'F', 'p-unc']]) # 对于多因素情况下的重复测量设计 multi_data = read_dataset('mixed_anova') result_multi_factors = rm_anova(dv='Scores', within=['Time', 'Metric'], subject='Student', data=multi_data, detailed=True) print(result_multi_factors[['Source', 'ddof1', 'ddof2', 'F', 'np2', 'eps', 'p-unc']]) ``` 上述代码片段展示了加载数据并调用 `pingouin.rm_anova()` 函数来进行一维或多维度上的重复测量方差分析的过程。该函数返回的结果包含了源(Source),分子自由度(ddof1),分母自由度(ddof2),F 值(F),未调整 p 值(p-unc)等要信息。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值