1、1124021321313004130010310200103020140201261220271201281111291201260110270021280301290400SPSS 配對調(diào)查資料的條件 Logistic 回歸分析（1:1 或 1:n）1 .1:1 病例對照研究的基本概念在管理工作中，我們也經(jīng)常要開展對照調(diào)查。例如為什么有的人患了胃癌，有的人卻不會患胃癌？如果在同一居住地選取同性別、年齡相差僅2 歲的健康人作對照調(diào)查，調(diào)查他們與患胃癌有關(guān)的各種影響因素，這就是醫(yī)學(xué)上很常用的所謂“1:1 病例對照研究”。病例對照研究資料常用條件 Logistic 回歸分析。條件 Logisti

2、c 回歸模型（conditionallogisticregressionmodel,CLRM）,下稱 CLRM 模型。2 .條件 Logistic 回歸模型的一個實(shí)例某地在腫瘤防治健康教育、社區(qū)干預(yù)工作中做了一項(xiàng)調(diào)查，內(nèi)容是三種生活因素與胃癌發(fā)病的關(guān)系。調(diào)查的三種生活因素取值見表 11-6。請擬合條件 Logistic 回歸模型，說明胃癌發(fā)病的主要危險因素。表 11-6 三種生活因素與胃癌發(fā)病關(guān)系的取值取值范圍X1（不良生活習(xí)慣）0,1,2,3,4 表示程度（0 表示無，4 表示很多）X2（喜吃鹵食和鹽腌食物）0,1,2,3,4 表示程度（0 表示不吃，4 表示喜歡吃、吃很多）表刁，差，1 表

3、不好表 11-750 對胃癌病例（S=1）與對照（S=0）三種生活習(xí)慣調(diào)查結(jié)果病例對照病例對照NoSX1X2X3NoSX1X2X3NoSX1X2X3NoSX1X2X3變量名3 .條件 Logistic 回歸模型的擬合原理與方法本例以 SPSS 軟件包來擬合 CLRM 模型。SPSS 雖有一個 LOGISTICREGRESSIO 麗令,卻無法擬合條件 Logistic 模型(CLRM),但有擬合 Cox 回歸模型的固有命令(COXREG。回顧一下 Cox 回歸模型的公式：h(t,X)=h0(t)exp(X131+X232+.+Xp3p)將兩側(cè)同時取對數(shù)可得：LN(h(t,X)=LN(h0(t)+

4、(X131+X232+.+Xp3p)如果是含有分層變量的 Cox 回歸模型，則公式變?yōu)椋篖N(h(t,X)=LN(h0i(t)+(X11+X22+.+Xpp)各層的基線風(fēng)險函數(shù) h0i(t)可以完全無關(guān)，而協(xié)變量的系數(shù)則在所有層中保持不變。作為半?yún)?shù)的方法，Cox 模型在擬合時并不估計基線風(fēng)險函數(shù) h0i(t),只估計各協(xié)變量的系51301612207131081300912201011001113001213401311111412211512301612411711101811311913412010202113212211002313002411112511205 00006 00107

5、 02108 02009010110020011001112032013020014002115020016000117001118000119020020000021031022020123022024001125020030131031110132142133140134120135112036120037101138100139130140120141120042130143121144120145111146101147121048120149112150120130002131000032010133020134000135020136020137011038040039001040

6、0301410101420001430000440100450001460000470000480110490001500031數(shù)值 3。這則與條件 Logistic 回歸模型不關(guān)心“i 的大小，只解出系數(shù)值 3 的思路一致。而且兩者都以最大似然法(LR)來擬合。把配對因素作為分層因素，即可消除配對因素的作用，因此可用分層變量控制法的原理，來擬合條件 Logistic 回歸模型。在使用 SPSS 的 COXRE 殺令擬合條件 Logistic 模型日(AnalyzefSurvivalfCoxRegression),要回答下列 5 個對話框：Time:給每行記錄一個虛才的生存時間（Time）,一

7、般默認(rèn)“病例”生存時間短（例如為 T）,“對照”生存時間長（例如為“2”）。只要“對照”的生存時間（例如“2”“100”）長于“病例”即可。因?yàn)?SPSS 會把在最短的完全數(shù)據(jù)的生存時間之前的截尾數(shù)據(jù)（Censored）全部去掉而不參加分析。輸出分層狀態(tài)時，Event=1;Censored=2。Status 與 DefineEvent:選入虛擬生存狀態(tài)變量。定義事件時病例全為1（Singlevalue）,為完全數(shù)據(jù)；“對照”全“0”，為截尾數(shù)據(jù)。此值“1”表示事件已發(fā)生，本例即“病例”。Covariates:選入要進(jìn)行分析的協(xié)變量（即自變量 X1,X2,X3）。Method:用 FSTEP（L

8、R）,即似然比法。Strata:錄入配對的對子號，“1”指第 1 個對子（含“病例”為“1”一行與“對照”者為“0”一行）；“2”指第 2 個對子。本例指變量“No”。4 .擬合條件 Logistic 回歸模型所用程序文件所用程序文件名為 CondLogiRegre.sps。ConditionalLogisticRegression;filename:CondLogiRegre.sps.*43112114302000441120144021004511111450200146110114602000471121047020004811201480211049111214902001501120

9、15002031ENDDATA.LISTVARIABLES=ALL/FORMAT=NUMBERED/CASES=FROM1TO50.* Whenprograming:.* Status:Case=1;Control=0.* Time:Case=1;Control=2,thatis,thesurvivaltimeofControlshouldlongerthanCase.* DefineEvent:Singlevalue:1meansCase.* Covariate:Inputetheindependencevariables.* Method:ENTER,FSTEP,BSTEP,etc.* S

10、trata:TheNo.ofperpairs,thatis,No.COXREGTime/STATUS=Status(1)/STRATA=No/METHOD=FSTEP(LR)X1X2X3*FangJi-Qian:Med.Stat.andComputerExperiment,P.492:.*UsingCOX-RegressiontofitConditionalLogisticRegression:DATALISTFREE/NoStatusTimex1x2x3.BEGINDATA.11124031130051130171131091122011113001311111151123017111101

11、911341211132123113002511120271120129112013111101331140135111203711011391130110231030220150200070221090210111020111302200150220017020111902200210231023022202502200270202129024003102000330220135022013702110390201021132141130061122081130010111001211340141122116112411811131201102022111002411111261122028

12、111113011310321142134112013611200381100140112012020104022016020108022001002200120232014020211602001180200120020002202201240201126021102802301300202132021013402001360220138024004002301/PRINT=DEFAULTCI(95)/CRITERIA=PIN(0.05)POUT(0.1)ITERATE(20).*5.擬合條件 Logistic 回歸模型的主要運(yùn)算結(jié)果運(yùn)行程序文件 CondLogiRegre.sps 可得到下

13、列主要結(jié)果:VariablesintheEquationBSEWalddfSig.Exp(B)95.0%CIforExp(B)LowerUpperStep1X10.6590.2208.97310.0031.9321.2562.973Step2X10.7850.2579.35110.0022.1931.3263.629X20.8140.3077.04210.0082.2571.2374.118回歸系數(shù) B 表示病例與對照變量值之差與患病優(yōu)勢的關(guān)系，即 Exp（B）表示病例與對照暴露水平相差1 個單位時患病的比數(shù)比（OddRatio,OR）。從表中可見，X3（精神狀況）未被引入回歸模型，只引入影響因素 X1 與 X2,說明精神狀況的好壞與胃癌發(fā)病無顯著相關(guān)，而 X1（不良生活習(xí)慣）和 X2（喜吃鹵食和鹽腌食物