如何正确运用χ2检验
——生存资料的三种特殊检验与SAS实现
2021-11-04胡纯严胡良平
胡纯严 ,胡良平 ,2*
(1.军事科学院研究生院,北京 100850;2.世界中医药学会联合会临床科研统计学专业委员会,北京 100029*通信作者:胡良平,E-mail:lphu927@163.com)
在对生存资料进行假设检验时,除了在试验因素各水平组之间采取似然比检验和六种非参数检验之外,还有三种特殊的检验,即“多重比较”“趋势检验”和“协变量检验”[1]。这些假设检验方法在本质上都属于“χ2检验”。本文将详细介绍前述提及的三种特殊检验的计算公式,结合实例并借助SAS软件实现统计计算。
1 多组生存函数之间的多重比较
1.1 何为多重比较
设有K(K>2)组生存函数,若同时对它们进行比较,就叫做整体比较;若每次仅取出两组进行比较,并按以下两种情形之一完成全部比较,就称为多重比较:情形一,任何两组生存函数之间都要进行一次比较,若K=4,总共需要比较6次;情形二,选定一个对照组,其他组都与此对照组进行比较,若K=4,总共需要比较3次。
1.2 多重比较的校正
1.2.1 校正的必要性
当需要对K(K>2)组生存函数进行多重比较时,基于前述提及的两种情形,都至少需要比较2次。若每次比较都设定检验水准为α(通常设α=0.05),设总共需要比较m次,m的计算公式如下:
于是,m次比较所对应的犯I类错误的总概率PT就会大于α,其计算方法[1]如下:
PT=1-(1-α)m
例如,当m=6、α=0.05时,PT=0.264908>0.05。
为此,需要对每次比较所得的概率P值进行校正,以保证整个多重比较所对应的犯I类错误的总概率PT小于等于α。
1.2.2 校正方法的种类
在 SAS/STAT 的“LIFETEST”过程中[1],基于不同的统计学原理,对多重比较结果进行概率校正的方法有如下几种:Bonferroni校正法、Dunnett-Hsu校正法、Scheffé校正法、Šidák校正法、SMM校正法和Tukey校正法。……
