NM条件logistic回归分析在统计软件上的实现.pdf

·93· N：M条件logistic回归分析在统计软件上的实现堆 广州南方医科大学公共卫生与热带医学学院生物统计学系(510515)郑一男 曹佩华 欧春泉△ 本文在STATA、R、SAS常用统计软件上实现N lj ：M条件logistic回归分析，并指出Cox风险比例模型 ．t exp(∑y“工∥) Pr(Yi 进行条件logistic回归分析时结点调整方法的适用情 “1 ∑exp(∑以邶) 况。 d-ESi f21 t 问题的提出 其中以等于0或1且∑d。=露“，置为七。。个1与 在有关疾病危险因素的流行病学调查研究中，常 也，个0所有可能组合的集合，显然这样的组合有 将病例与对照按年龄、性别等因素进行匹配，以消除某 I 7’．、 I。。1种。为精简计算，采用递推法计算(1)式的分母， 些混杂因素对研究结果的干扰作用¨。J。大家熟悉的 ＼庀“／ 是1：1或1：M匹配，即给每个病例选择1个或M个对 Ti 照与之匹配。匹配后往往再按匹配因素进行分层分 令Z(I州k)=∑exp(∑ditx，∥) diE5i t=l 析，这样可使每一个匹配层中都有一定数目的病例与 对照。或者直接采用N：M匹配，每一匹配层中有N 个病例和M个对照，而且各层中病例数和对照数的比 由此得到条件对数似然函数lnL【4】，进而采用最大似 例可以不固定，研究实施更为灵活。配对的病例一对照 然法估计参数。 研究数据均采用条件logistic回归分析。关于l：1和l t t In (2) ：M条件logistic回归分析的软件应用有不少讨论，但 L=∑{∑YitXit卢一l。甑(t，¨) 对N：M条件logistic回归分析却鲜有介绍，笔者在实2．Cox比例风险模型 际工作中遇到不少非统计专业的研究人员错误地运用 分析【5】。由于同一层内的生存时间均相同，故生存时 SPSS软件的Cox模块进行N：M条件logistic回归分 析，导致结果出现较大偏差。本文将全面介绍STA． 间的结点(ties)问题非常突出，如何处理结点是分析 TA、R和SAS统计软件如何对N：M配对病例一对照研 的一个关键。目前主流统计软件中有四种处理方法： Breslow、Efton、Exact 究资料进行条件logistic回归分析，并详细阐述在使用 marginal Cox模块进行条件logistic回归分析时的注意事项。然法)与Exactpartial 读者可根据自身的使用偏好选择任一软件进行分析。 但计算速度略慢，两种方法在处理大量结点数据时均 原理及方法 会使参数估计结果出现严重偏差。精确边际似然法实 1．条件l