一般来说，JAMA、新英格兰医学杂志等顶级论文，其统计分析计划，都非常详细，动辄5、60页。

但我今天要说的这篇JAMA主刊的论文，统计分析计划，却只有4页纸，但却用浓缩的语言，给我们讲述了针对不同统计指标，如何进行统计分析计划。这是一份非常值得学习的统计分析方案。

我现在来全景式展现这篇文章的分析过程！

一、研究论文介绍

2021年10月份，JAMA主刊刊登了一篇文章，比较的是中度和轻度低温治疗对心脏骤停幸存者的影响研究。

这是一项单中心、双盲、随机对照研究研究。从2013 年8月到20203月，共纳入330例院外心脏骤停幸存者。患者1:1进行随机分组，分别接受两种低温治疗方法。干预组是中度低温（31℃），对照组是传统的轻度低温（34℃），研究希望探讨中度低温是否更有利于改善患者的临床结局。

该研究的主要结局是患者180日内容死亡或者神经系统恶化情况（复合结局），次要包括一系列的定量、定性的指标，共19个

结果发现，中度低温（31℃）和轻度低温（34℃）治疗24小时在180天内的死亡率没有统计学差异（阳性结局发生率48.4% vs 45.4，RR=1.07 [95% CI, 0.86-1.33]; P = 0.56)。在总共19个次要结局中，18个没有统计学差异，只有ICU停留时间有统计学差异，中度低温（31℃）治疗ICU时间较长（10 vs 7天）。

因此可以总体上说，中度低温（31℃）和轻度低温（34℃）两种治疗手段在临床疗效上并没有差异。

二、它的统计分析方法

这篇文章在附件提供了统计分析计划，只有4页纸，作者也够懒的，幸好他们是加拿大的研究团队。

不过，这4页纸的统计分析计划，高度浓缩，针对研究的总共20个结局指标（1个主要，19个次要），分门别类做了介绍。而这些统计学方法也是我们日常所接触的经典方法。

为什么我说它是值得学习呢？因为我们国人写统计分析方案往往也很精要，因此这篇文章的方案就可以拿过来学习、模仿和改进了！

数据分析集

首先，临床试验首先要鉴定数据分析集。作者在统计分析计划，详细介绍了三种分析集：Modified intent-to-treat (mITT)、Population:Intent-to-treat (ITT) Population、Per-protocol Population (PP)。而该研究主要采用的是mITT，其他两种分析集作为敏感性分析所用。

The primary analysis was performed including all patientswho were randomized to a study therapy group and receivedthe therapy。

基线描述与差异性比较

连续型资料用t检验或者秩和检验，分类数据用卡方或者Fisher检验进行比较。Continuous variables will be tested for baseline comparability between the therapy groups moderate hypothermia and mild hypothermia using the Student’s t-test or the Wilcoxon rank- sum test. Categorical variables will be tested for baseline comparability with the chi-square test or Fisher’s exact test.

主要结局指标的分析方法

主要结局指标采用2套方法来进行，一套为主方法，另外一套是敏感性分析。由于是二分类结局，简单的方法采用卡方检验，复杂的方法采用logistic回归来分析。

第一套方法，是卡方检验

For the primary outcome analysis, the χ2 test was usedto compare the 31 °C and 34°C targettemperature groups, and relative risks (RRs) and 95% CIs werecalculated aswell as risk differences and their 95% CIs.

第二套方法是logistic回归

A logistic regression analysis was conducted to assess the consistency of the target temperature on the primary outcome, taking the randomization stratification factor (ventricular fibrillation or pulseless ventricular tachycardia initialrhythm) into consideration

次要结局指标

次要结局指标共19种，类型非常多。作者也一一做了解释

对于生存时间资料，采用生存分析的两种：logRank和Cox回归。

In addition, the cumulative incidence of mortality at 180 days was estimated using theKaplan-Meier method and hazard ratios (HRs) and 95% CIswere calculated.

该研究，有离散型的定量变量，离散型定量变量可以考虑采用类poisson回归的对数回归分析，计算RR值（generalized linear model using the binomial distribution and a log link）

对于连续型定量数据，则采用t检验、秩和检验和线性回归分析。

第一套方法是t检验和秩和检验

For these analysis, the mean, standard deviation, median and IQR for each secondary endpoint will be determined and 95% CIs calculated; the Student’s t-test or the Wilcoxon rank-sum test will be considered.

第二套方法是多重线性回归

In addition, a multiple regression model will be used to assess the consistency of the therapy group effect on the primary outcome taking the randomization stratification factor (VT/VF initial rhythm) into consideration;

安全性分析

对于研究可能存在着不良反应，作者采用卡方或者Fisher法进行两组间的比较。

If number of events is sufficient, moderate TH and mild TH therapy groups will be compare on these adverse events using the chi-square test or Fisher’s exact test

缺失数据处理

临床研究，缺失数据不可避免，作者的开展了缺失数据的填补，采用多重填补法

Missingness is considered to be missing at random (MAR) and multiple imputation (MI) techniques will be used for handling missing data. Multiple imputations were used for any missing outcomes usingthe Markov chain Monte Carlo method with Jeffreys noninformative prior distribution.

多重比较的一类错误控制

既然本研究同时有多达20个指标要进行比较，就不得不思考一定的方法矫正一类错误（即假阳性率）了。不过作者没有矫正次要结局指标的假阳性错误，所以他说结果是探索性的。

Because of the potential for type I error due to multiple comparisons, findings for analyses of secondary end pointsshould be interpreted as exploratory。

三、小结

所有的高质量临床试验论文，对统计学方法要“如数家珍”，什么结局指标用什么统计学方法。其方法不见得复杂，但是得有针对性。这篇文章涉及了连续型定量、离散型定量、分类变量的数据分析，所以囊括了几乎所有传统的方法。

我之所以推崇它，是因为它的统计学分析计划，麻雀虽小五脏俱全，是我们很多临床研究设计者学习和模仿的材料。

现在，咱们的大夫越来越懂得统计了，但是我们得知道如何将它合理地嵌入到科研设计理念中去，要把一个统计故事讲完整了！

本文转自：医学论文与统计分析

往期精品(点击图片直达文字对应教程)

机器学习

后台回复“生信宝典福利第一波”或点击阅读原文获取教程合集