如何用spss17.0進行二元和多元logistic回歸分析

#歸分析歸分析如何用spss17?0進行二元和多元logistic―、二元logistic歸分析二元logistic回歸分析的前提為因變量是可以轉(zhuǎn)化為0、1的二分變量，如:死亡或者生存，男性或者女性，有或無，Yes或No,是或否的情況。下面以醫(yī)學中不同類型腦梗塞與年齡和性別之間的相互關系來進行二元logistic回歸分析。（―）數(shù)據(jù)準備和SPSS選項設置第一步原始數(shù)據(jù)的轉(zhuǎn)化:如圖1-1所示具中腦梗塞可以分為ICAS、ECAS和NCAS三種，但現(xiàn)在我們僅考慮性別和年齡與ICAS的關系，因此將分組數(shù)據(jù)ICAS、ECAS和NCAS轉(zhuǎn)化為1、0分類，是ICAS賦值為1,否賦值為0。年齡為數(shù)值變量，可直接輸入到spss中，而性別需要轉(zhuǎn)化為（1、0）分類變量輸入到spss當中，假設男性為1,女性為0,但在后續(xù)分析中系統(tǒng)會將1,0置換（下面還會介紹）,因此為方便期間我們這里先將男女賦值置換，即男性為"0”,女性為“1”I"e?a曲□並噸零母?歩I"H書.iCAS[-QMIrtftE(KM116301.1541130TH[111.3700*116100.?1n1切901陶11115501.023n180its1E.190□I11KQ1.154aHD14570011GO021130?DD和C>00*11720OKi3501601W1eiu11114&1a752闌n17D100[iS.S71011旳11.1?dWi1?01n.w1011E3Q1.1744BaI6Usa04?0n1173QIE：43ai：3n14葩012B2061QQ4?i扣iHQ15JH1i|20詢1i：i?1闌nUQ10015.J611205701.D03121no也&7.?1f70羽□03912Di^a弐］15iSB0□20陌1L5I44B1120蟲］0叫112057Q1血34B114Dsa15330020?800J4212nI切100[IB.921120肝0a73261価9010J?0o|20隔11.?3網(wǎng)n戚別11ll205302113]1E9D00305711VEB115Uioa1迪1120E?□t&j4品0iSl兇0iJijOB1120330Q702nUQ901B.'C0120旳1餛41110tfl15.i71120Iow2□14Dso01341II-4j*幻電n□垃it于一if:獻律腿可筋制P"…4■耳岀1|T44H-5W5S..|隔POoifHi七&I壘=Ark&口:i比或Pir44Sii.e圖1-1第二步：打開“二值Logistic回歸分析”對話框：沿著主菜單的“分析（Analyze）—回歸（Regression）—二元logistic（BinaryLogistic）”的路徑（圖1-2）打開二值Logistic回歸分析選項框（圖1-3）。如圖1-3左側(cè)對話框中有許多變量，但在單因素方差分析中與ICAS顯著相關的為性別、年齡、有無高血壓，有無糖尿病等（PvO.05），因此我們這里選擇以性別和年齡為例進行分析。

J-llZ.電ffw|二已砂iamr詡尸..]?豈tpmgi卜E^Jim-aamiauttah山：.二紀4石41^血七..[當:MffiiwiT.Ga打二…S.爐畀幺吞矽程牌圖1-2圖1-3在圖1-3中，因為我們要分析性別和年齡與ICAS的相關程度，因此將ICAS選入因變量(Dependent)中，而將性別和年齡選入?yún)f(xié)變量(Covariates)框中，在協(xié)變量下方的“方法(Method)”一欄中，共有七個選項。采用第一種方法，即系統(tǒng)默認的強迫回歸方法(進入“Enter”)。接下來我們將對分類(Categorical)，保存(Save)，選項(Options)按照如圖1-4、1-5、1-6中所示進行設置。在“分類”對話框中，因為性別為二分類變量，因此將其選入分類協(xié)變量中，參考類別為在分析中是以最小數(shù)值“0(第一個)”作為參考，還是將最大數(shù)值“1(最后一個)”作為參考，這里我們選擇第一個'0”作為參考。在“存放”選項框中是指將不將數(shù)據(jù)輸出到編輯顯示區(qū)中。在“先項”對話框中要勾選如圖幾項，其中exp(B)的CI(X)”一定要勾選，這個就是輸出的0R和CI值，后面的95%為系統(tǒng)默認，不需要更改。

圖1-4分類圖1-5保存圖1-4分類圖1-5保存圖1-6選項另外在“選項”對話框中，“輸出”一欄中，系統(tǒng)默認為“在每個步驟中”，這里更改為“在最后一個步驟中”，即：輸出結(jié)果將僅僅給出最終結(jié)果，而省略每一步的計算過程。由于我們采用強迫回歸，逐步回歸概率選項可以不管此外還有一個選項需要說明。一是分類臨界值(Classificationcutoff)，默認值為0.5,即按四舍五入的原則將概率預測值化為0或者1。如果將數(shù)值改為0.6，則大于等于0.6的概率值才表示為1,否則為0。其情況余依此類推。二是最大迭代值(MaximumIterations),規(guī)定系統(tǒng)運算的迭代次數(shù)，默認值為20次，為安全起見，我們將迭代次數(shù)增加到50。原因是，有時迭代次數(shù)太少，計算結(jié)果不能真正收斂。三是模型中包括常數(shù)項(Includeconstantinmodel)，即模型中保留截距。除了迭代次數(shù)之外，其余兩個選項均采用系統(tǒng)默認值。完成后，點擊各項中“繼續(xù)(Continue)”按鈕。返回圖1-3,單擊“確定”按鈕。(二)結(jié)果解讀其他結(jié)果參照文章《利用SPSS進行Logistic回歸分析》中解讀，這里重點將兩點：第一，分類變量編碼(圖1-7)，由于這里包括性別分類變量，而我們對性別賦值為1和0,但在spss中系統(tǒng)會默認把我們的數(shù)值進行置換，即1—參數(shù)編碼0,0—參數(shù)編碼1,而最終輸出結(jié)果是以1來計算的，而0為參考數(shù)據(jù)。所以這也就是為什么我么之前要對研究組男性的賦值進行置換了。如果男性為1那么spss中最終輸出的將是女性的分析結(jié)果。頻率參數(shù)編碼(1)性別a521.000129.000圖1-7第二，最終輸出數(shù)據(jù)(圖1-8)在該結(jié)果中，Exp(B)即為文獻中提及的OR值，而EXP(B)的95%C.I.即為文獻中提及的CI值。其中Exp(B)表示某因素(自變量)內(nèi)該類別是其相應參考類別具有某種傾向性的倍數(shù)。而有的文獻中提到的CrodeOR和AdjustOR則分別為單因素優(yōu)勢率(Crodeoddsratio)和多因素優(yōu)勢率(Adjustoddsratio)，即僅對性別單個變量的單因素分析或者對性別和年齡等多個變量進行多因素分析后所得到的不同結(jié)果oCI則為可信區(qū)間(Confidenceinterval)。Sig.即我們常說的P值,P<0.05為顯著(無效假說不成立，具有統(tǒng)計學意義)，P>0.05為不顯著(無效假說成立，不具有統(tǒng)計學意義)。

B8.E.WaisdfSig.Exp⑶EXF1向的95%C.1.下限二匪卿卩年齡0^3.0222.2771.1311.0339901.C7S性別⑴3：0.503.4521.5011.^103.5233.7C3-2.7941.39540141.045.061在歩賤1中錨鼻的吏呈:干賂哇別.圖1-8歸分析二、多項（多元、多分類、Multinomial）logistic歸分析前面講的二元logistic回歸分析僅適合因變量Y只有兩種取值（二分類）的情況，當Y具有兩種以上的取值時，就要用多項logistic回歸（MutinomialLogisticRegression）分析了。這種分析不僅可以用于醫(yī)療領域，也可以用于社會學、經(jīng)濟學、農(nóng)業(yè)研究等多個領域。如不同階段（初一、初二、初三）學生視力下降程度，不同齲齒情況（輕度、中度、重度）下與刷牙、飲食、年齡的關系等。下面我們以圖1-2中，對apoba1（ApoB/AI）項中數(shù)值做四分位數(shù)后，將病人的ApoB/AI的比值劃分為低、較低、中、高四個分位后利用多項logistic回歸分析其與ICAS之間的相互關系。首先來做四分位數(shù)，很多人在做四分位數(shù)的時候都是自己算出來的，其實在SPSS里面給出了做四分位數(shù)的程度即分析（Aanlyze）-描述統(tǒng)計（DescriptiveStatistics）f頻率（Frequencies）。打如圖2-1開頻率對話框。將我們要分析的數(shù)值變量Apobal選入到變量對話框中。圖2-1選擇統(tǒng)計量扌安照圖2-2中勾選四分位數(shù)選項其他選項按照自己需要勾選，

然后點擊圖2-1中的確定按鈕，開始運算。在圖2-3中可以讀取我們的四分位

數(shù)值。圖中百分數(shù)表示的是對該變量做的四分位數(shù)的百分比，25表示前25%的，50表示前50%的，75表示前75%的。每一項對應的后面數(shù)值即為相應的四分位數(shù)，如0.5904，即為前25%的個體與后75%個體的分位數(shù)。N有效81N有效81缺失0百分位數(shù)25.590450.8600751.0006Am詢圖2-2圖2-3按照如上方法得出ApoB/AI的比率后我們可以把該比值劃分為四個區(qū)間，即當ApoB/AI的比率＜0.5904為低、當0.5904＜ApoB/AI的比率＜0.88時為較低、當0.89＜ApoB/AI的比率＜1.0886時為中，當ApoB/AI的比率＞1.0886時為高。然后將這一劃分如圖1-1中“四分位數(shù)”一項用分類數(shù)值表示即1代表低，2代表較低，3代表中，4代表高。這里還要強調(diào)的是我們要研究其與ICAS之間的相互關系，那么我們需要將其設為二分類變量，即是ICAS的情況為1,否則為0,但多項logistic回歸分析也會將1,0置換，所以我們需要在這里將我們需要研究的情況置換為0,然后將其他置換為1。下面就可以進行多項logistic回歸分析了。如圖2-4打開多項logistic回歸分析對話框（圖2-5）。如圖2-5所示，在”因變量”中選入剛才我們輸入的四分位數(shù)分類變量，在因子中輸入分類變量ICAS（這里一定是分類變量，可以是一個也可以是多個）,在協(xié)變量”中輸入數(shù)值變量如年齡（這里一定是數(shù)值變量，可以是一個也可以是多個），但因本次沒有對年齡進行分析，僅對ICAS進行了單因素分析，所以我們把年齡移出協(xié)變量選項。在SPSS中對因變量的定義是，如果因變量Y有J個值（即Y有J類），以其中一個類別作為參考類別，其他類別都同他相比較生成J-1個冗余的Logit變換模型，而作為參考類別的其模型中所有系數(shù)均為0。在SPSS中可以對所選因變量的參考類別進行設置，如圖2-5在因變量對話框下有一“參考類別”選項。點擊后會彈出圖2-6對話框。在該對話框中我們選中設定，輸入數(shù)值1,這代表

'-rmHv我們以分類數(shù)值1所代表的類別作為參考類別，即最低數(shù)值作為參考類別。單擊繼續(xù)。當然也可以選擇“第一類別”和“最后類別”，入選中分別表示以最低數(shù)值或最高數(shù)值作為參考類別。其他設置與二元Logistic分析相似，將我們要輸出的項勾選即可，點擊圖2-5中確定，輸出數(shù)據(jù)。'-rmHvMf制'ei:悶W1jm：ri-jMHUe)rxittHstO響厲畫直EE2匝tfliati110*D90I50EODIQGBII亠21口￥45D豈1I20EC09SO口031D￡B0i陽1w他0恒宙11410T90I3：i0180K1価10S1012411106001BJ6001610nn>r>p盤如472□IBCI艮1D|II11r-l14565?0I-E.itlHK聲?017Dal1001a&112DEd15.771a|B1016SCim閒n1D.7211^*iOID10240*140鮒n10.710o|衛(wèi)*=ii10rSBiSj□1%snnlEfifl1i|FAW*-Jl臼曲121口K□IGDsDID€E1i|13|0RMLl2□22D130n108(111I4102i01001100ma10I510bJR：fi理二JSiiMk或0140備u10.6311IE10園峽鼻*WtSl￡l6Cl1<i110i6.000-0|171口ttUQLi4S0I4C<EC012^1口□I101口舊D144□I4C>ECDins(]1llIB1Q.ViImW'2.1￥B30113DHl0TISP□a|為10610ite01S0ffl15.491121101.023場01001B01B.190o|雄I01.15J65PQl?si1收0□|2310Q7B2eaD24II3D7DDZJ.ID0□|1口Cl913?2D豈0IGOl?]11i|10a両21960170ira05B?1D*SrfBiD削0悶t>Hi[u**iTj*1他j田駐遜耳Ji珂戯色Eo^V4扌程最阿[■■*!]?wss職岳看命￡切獨疏潘■卜■1ZtH*叫「廠_廠命耐|番」宀障相ms辭J創(chuàng)頤azT.a..].丄.*曲4葉..|口口卿曲s*...|対矽血y...|二圖2-4圖2-5圖2-6輸出數(shù)據(jù)基本與二元Logistic分析相似，我們重點講下最后一項“參考估計”如圖2-7所示，其中參考類別為ICAS=1的分類情況，而其中的ICAS=0分為2、3、4三種，分別給出了ICAS=0時的數(shù)值。而其中Exp(B)(即OR值)表示某因素（自變量）內(nèi)該類別是其相應參考類別具有某種傾向性的倍數(shù)。如Exp（B）=2.235時，即表示在較輕這一類別下ICAS患者數(shù)為其他類別（ECAS和NCAS）的2.235倍。這里面的顯著水平即為P值。這里要強調(diào)的是，一些文獻中在輸出數(shù)據(jù)的時候經(jīng)常會給出“Referent（參考）"項，這里的Referent，即為我們這里所選的參考類別1,因為1作為參考類別，所以其所有數(shù)值為0,即無數(shù)據(jù)輸出。因此