临床研究常用统计方法概述金雪娟周俊时智英葛均波(复旦大学附属中山医院,上海市心血管病研究所,上海200032)经过周密设计和科学实施的临床研究还需要规范的数据管理和统计分析,才能得到可靠的结论。随着计算机技术和统计分析软件发展,近年来,统计理论和方法发展非常迅速。临床医师日常繁忙的工作使得他们很少有时间系统学习医学统计理论,及时了解一些实用、有效的新方法。在此,我们介绍目前临床研究最常用的一些统计分析方法,以实用、易懂为原则,重点综述各种方法的适用条件。1几个基本概念和统计量1.1数据的类型数据(Data)是统计分析的基础。统计分析方法的选择取决于不同的数据类型。最常见的数据类型有两种,分类数据(CategoricalData)或称定性数据(QualitativeData)和定量数据(QulantitativeData)或称计量数据(NumericalData)。分类数据类型:分类数据的分层大于2时,又称为多分类数据(PolytomousData)。分类数据类型有无序(NominalCategorial)和有序(OrdinalCategorieal)。无序数据如性别(男、女)、血型(A、B、O、AB型)等。有序数据如肿瘤的分级(I级、II级、III级)、疼痛的程度(轻、中、重)等,以及在临床研究设计中,经常看到的“非常好、好、一般、差”这样的数据类型。不同类型的分类数据在统计分析方法上也不同,并不是大家所熟悉的x2检验所能全部涵盖的。定量数据类型:包括连续性数据(ContinuousData),如身高、体重以及不连续性数据(DiscreteData),如妇女的产次,疾病的复发次数等。1.2常用的描述性统计量最常用的描述集中趋势的统计量为算术均数(ArithmeticMean),但其值易受极端值影响。可以采用中位数(Median)、修整均数(TrimmedMean,去除最大和最小值后的算术均数)或Winsorized均数(WinsorizedMean,极端值用最接近的非极端值替代后的算术均数)来代替。对于数值呈几何分布的资料,则可采用几何均数(GeometricMean)。临床研究论文中常采用均数±标准差或均数±标准误来表示定量数据的分布特征。标准差(StandardDeviation)为方差(Variance)的平方根,表示个体数值与样本均数间的离散程度;标准误(StandardError)为均数的标准差,表示样本统计量与总体参数间的离散程度,标准误越小,总体均数的95%可信区间(confidentinterval,CI)越窄,也就是说样本均数对总体均数的代表性越好。虽然不同的统计学家对论文中应该引用哪种表达方式有争议,但两种方式均用于描述正态分布的计量数据。在医学论文中,采用标准差或标准误应该说明。对于非对称数据只用均数±标准差或标准误表达是不恰当的,可以采用中位数结合四分位数间距(Inter-quartileRange)表示。1.3显著性水平(a)和P值为了说明统计学显著性差异,通常把P值与参考值a比较,如果Pa,则拒绝H0,(按此a水平,存在统计学显著性差异)。a水平通常选用0.05,其实0.05并非是一个具有魔力的数字,其它常用的a值还有0.01和0.10等,甚至在某些特殊的情况下,a值可以取0.20等。a值说明I类错误的概率,如a等于0.05,则I类错误的概率为5%(拒绝正确H0的概率)。研究论文报告结果时,应该把a和P值,同时加以说明。参阅已经发表的文献,常发现这样的问题,研究者设立了A、B和C组3个组,在统计检验时,先把A与B比较,然后B与C比较,再A与C比较,这是需要极力避免的。多重比较会增加I类错误的概率。假设有K个中心,同时参与某临床实验,如果按各个中心分别统计,则a=1-(1-0.05)K。当K等于1、2、3、4时,a值依次为0.050、0.098、0.243、0.185,随着独立检验次数的增加,a值也增加。1.4统计分析软件的选择统计软件非常多,目前临床研究常用并且比较权威的统计软件有SAS(StatisticsAnalysisSystem)、SPSS(StatisticalPackagefortheSocialScience)、STATA(Statisties/DataAnalysis)和Splus等。STATA灵巧方便,价格也能为个人用户所承受。SSPS的菜单式操作,使用简便,而且介绍SPSS的书籍比较多,目前已经成为国内非统计专业人员统计的首选软件。SAS是主要针对专业统计用户设计的软件,在数据处理和统计分析领域,被誉为国际上的标准软件系统。国际上大部分著名高校和生物统计机构均使用SAS作为统计分析工具,一些最新的统计方法在SPSS和STAT中没有包括,需要选用SAS处理。2临床研究常用统计方法选择的原则2.1单组资料的分析如果数据呈正态分布,则采用单样本t检验(One-SampleTest);如果数据呈非正态分布,可以采用相应的非参数统计方法Wilcoxon符号等级检验(WilcoxonSigned-RankTest)。2.2两组资料的分析首先,需要判别数据为定量数据还是分类数据。如果是定量数据,数据的分布特征呈正态,则选择两样本t检验(Two-SampleTest);如果定量数据呈非正态分布,则选择Wilcoxon等级和检验(WileoxonRankSumTest)。对于分类数据,卡方检验(Chi-SquareTest)被广泛运用。但值得注意的是,如果行列表中有1/5以上的格子理论频数小于5,或有飞个格子理论频数小于1,卡方检验将导致分析的偏性。此时,可以采用Fishe's精确概率法(Fishe'sExactTest)计算P值。2.3三组或以上资料的分析与两组资料分析类似,如果数据为定量资料呈正态分布,则采用单因素方差分析(OnewayANOVA);如果定量数据,呈非正态分布,则选择Kruskal一Wallis检验(Kruskal一WallisTest)。对于分类数据,多分类无序数据采用卡方检验(Chi-SquareTest)或Fish's精确概率法;多分类有序数据可采用Cochran-Mantel-Haenszel检验(Cochran-Mantel-HaenszelTest)。2.4生存分析(AnalyzingResearehQuestionsaboutSurvival)分析一段时间后生存、死亡或其它事件发生情况需要采用生存分析,例如,研究者想了解心脏移植后病人生存天数是否与不同的手术方式有关。生存分析的目的通常是为了描述研究人群的事件发生时间(生存时间、suvivaltime)的分布特征,比较不同组的生存时间或研究生存时间是否与研究变量有关。单因素生存分析可以采用Log-Rank检验(Log-RankTest);多因素可以考虑选择比例风险模型(CoxProportionalHazardsModel)。需要注意的是,在临床研究中经常包含重复测量数据,如病人从心脏移植至死亡发生期间,重复测量了多次心功能值,对于这种资料,可以采用SASPHREG中,重复测量资料的cox模型的运用。2.5相关性分析(ResearehQuestionSaboutRelationshipsamongVariables)许多临床研究涉及对一组研究对象2个连续性变量的相互关系的研究。如用2个不同的指标测定心功能,拟评价这2个指标是否一致,这就需要涉及相关(Correlation)和一致性(Agreement)的评价。当数据为正态分布时,Pearson相关系数(Pearson'srelationCoefficient)可以评价2个指标的相关性。当数据分布非正态,相应的非参数统计量为Spearman's等级相关系数(Spearman'sRankCorrelationCoeffielent)以及Kendall'sTau-b等级相关系数(Kendall'sTau-bRankCorrelationCoeffielent),两者类似,但更多的统计学家推荐使用后者。对一致性的评价,定量数据可以采用Concordance相关系数(ConcordanceCorrelationCoeffieient),分类数据采用Kappa分析(Cohen'sKappastatistic)。值得注意的是,相关性与一致性的区别。如在临床研究中,希望评价一个新的方法是否等同于原来的方法,需要使用一致性分析。2.6多因素分析(StatlstiealMethodSfo:MultipleVariables)临床研究的对象常常为病人,与有严格实验条件控制的动物实验不同,除了研究的因素外,常需要控制许多混杂因素或协变量,统计分析需要采用多因素模型对协变量进行校正。统计分析软件和程序的使用为多因素分析提供了可能。根据反应变量的类型,可以采用多元线形回归(MuiripleLinearRegression)、协方差分析(ANOVA,AnalysisofCovariance)以及Logistic回归(LogisticRegression)等。2.7重复测量数据的分析(MethodsforAnalyzingRepeatedMeasuresData)定量数据可以采用重复测量方差分析(Repeated-MeasuresANOVA)以及混合效应模型(MixeD-effectslinearModel)。对于分类数据,可以广义估算方程(GeneralizedEstimatedEvluation,GEE)拟合Logistic模型。表1总结了常用的多因素分析方法以及选择原则以及对应的SAS步。表1多因素分析方法与SAS自变量反应变量----------------------------------------连续性数据分类数据重复测量计量数据(正态分布)多元线性回归(PROCREG)多因素方差分析(PROCGLE)重复测量方差分析(PROCMIXED)分类数据(无序)Logistic回归(PRCOLOGISTIC)Logistic回归(PRCOLOGISTIC)重复测量Logistic回归(PRCOLOGISTIC)分类数据(有序)比例比Logistic回归(PRCOLOGISTIC)比例比Logistic回归(PRCOLOGISTIC)重复测量比例比(PRCOGENMOD)生存分析COX回归(PROCPHREG)Logrank/WilcoxonCox回归(PROCLIFETEST或PROCPHREG)重复测量Cox回归(SUIDAAN)3.小结临床工作者花费了大量心血收集到的宝贵资料,需要有适当的统计分析和结果呈现。目前国内统计书籍中介绍的方法常常不能满足充分利用研究信息的需要。我们对各种常用统计方法的选择作了概述,实际运用时,还需要参阅相应的SAS程序。