第十三章

二項分類變量與多項有序分類變量關(guān)系的分析——二項分類Logistic回歸分析主講老師：陳青山數(shù)據(jù)變量類型：影響變量：多項有序分類分變量；結(jié)果變量：二項分類變量。分析目的：不同等級水平對結(jié)果是否發(fā)生的影響及其影響趨勢例子：不同新生兒出生體重是否影響新生兒黃疸的發(fā)生不同劑量組的療效是否不同基本原理二項分類的Logistic回歸模型：利用回歸系數(shù)或比值比OR值的大小，探討影響變量的作用及其作用強度。當(dāng)，比值比OR=1時，表示影響變量對疾病的發(fā)生不起作用；當(dāng)，比值比OR>1時，表示該影響因素為危險因素；當(dāng)，比值比OR<1時，表示該影響因素為保護因素。Logistic回歸分析需要對擬合的Logistic回歸模型進行檢驗，主要包括3部分：對回歸系數(shù)的檢驗；對Logistic回歸模型的擬合優(yōu)度檢驗；對Logistic回歸模型預(yù)測準(zhǔn)確度的檢驗。第一節(jié)數(shù)據(jù)庫數(shù)據(jù)的分析結(jié)果變量和影響變量用數(shù)據(jù)庫形式給出的二項分類變量與多項有序分類變量的分析。分析目的：分析不同等級水平對結(jié)果是否發(fā)生的影響及其影響趨勢實例數(shù)據(jù)特征：影響變量：多項有序分類變量，體重分組（group=1，2，3），有低體重兒、正常體重兒和巨大兒3結(jié)果變量：二項分類變量，新生兒黃疸是否發(fā)生（effect=1，0），有發(fā)生和未發(fā)生2種。分析目的：了解新生兒黃疸是否發(fā)生與新生兒出生體重的關(guān)系。統(tǒng)計方法：應(yīng)用二項分類Logistic回歸模型分析。實例分析SAS的計算程序SAS的計算結(jié)果反應(yīng)變量的響應(yīng)剖面（ResponseProfile）：反應(yīng)變量的水平排序為升序按照effect=1的概率擬合模型。SAS中對模型是否成立的檢驗使用了似然比檢驗、Score檢驗和Wald檢驗三種方法。似然比檢驗結(jié)果x2=39.65，P<0.05；Score檢驗結(jié)果x2=37.51，P<0.05；Wald檢驗結(jié)果x2=34.55，P<0.05。P值均小于0.05，故可以認(rèn)為模型成立。正常體重兒組（group2）與低體重兒組（group1）的OR值為0.177，OR95%CI為（0.091，0.347），差異有統(tǒng)計學(xué)意義；巨大兒組（group3）與低體重兒組（group1）的OR值為0.547，OR95%CI為（0.249，1.199），差異無統(tǒng)計學(xué)意義。以上結(jié)果表明，正常體重兒患新生兒黃疸的風(fēng)險是低體重兒的0.177倍，或者說，低體重兒患黃疸的風(fēng)險較高，是正常體重兒的5.650倍。SPSS的操作步驟：SPSS的計算結(jié)果：似然比x2=39.65，P=0.000<0.05，

故可認(rèn)為該模型成立。

（以低體重兒組為參照）：正常體重兒組（group2）與低體重兒組（group1）的OR值為0.177，OR95%CI為（0.091，0.347），差異有統(tǒng)計學(xué)意義；巨大兒組（group3）與低體重兒組（group1）的OR值為0.547，OR95%CI為（0.249，1.199），差異無統(tǒng)計學(xué)意義。以上結(jié)果表明正常體重兒患新生兒黃疸的風(fēng)險是低體重兒的0.177倍，即低體重兒患黃疸的風(fēng)險是正常體重兒的5.650倍。第二節(jié)頻數(shù)表數(shù)據(jù)的分析結(jié)果變量為二項分類變量，影響變量為多項有序分類變量，且應(yīng)用交叉頻數(shù)表形式給出的數(shù)據(jù)分析。分析目的：分析不同等級水平對結(jié)果是否發(fā)生的影響及其影響趨勢實例數(shù)據(jù)特征：影響變量：多項有序分類變量，劑量組（group=1，2，3），有高劑量組、中劑量組、低劑量組3組；結(jié)果變量：二項分類變量為療效（effect=1，0），分為有效、無效。分析目的：分析使用一段時間后療效與劑量高低的關(guān)系。統(tǒng)計方法：應(yīng)用二項分類Logistic回歸模型分析。實例分析SAS的計算程序SAS的計算結(jié)果反應(yīng)變量的響應(yīng)剖面（ResponseProfile）：反應(yīng)變量的水平排序為升序按照effect=1的概率擬合模型SAS中對模型是否成立的檢驗使用了似然比檢驗、Score檢驗和Wald檢驗三種方法。似然比檢驗結(jié)果x2=8.57，P=0.0138<0.05；Score檢驗結(jié)果x2=8.65，P=0.0135<0.05；Wald檢驗結(jié)果x2=8.28，P=0.0159<0.05。P值均小于0.05，故可以認(rèn)為模型成立。

二項分類Logistic回歸分析得（以低劑量組為參照）：高劑量組（group1）與低劑量組（group3）的OR值為3.136，OR95%CI為（1.389，7.082）；而中劑量組（group2），與低劑量組（group3）的OR值為1.970，OR95%CI為（0.945，4.104。以上結(jié)果表明，高劑量組的療效高于低劑量組，是低劑量組的3.136倍；盡管中劑量組的療效是低劑量組的1.970倍，但因其OR95%CI包含1，尚不能認(rèn)為中劑量組療效比低劑量組高。SPSS的操作步驟SPSS的計算結(jié)果似然比x2=8.57，P=0.014<0.05，故可認(rèn)為該模型成立。

參數(shù)檢驗結(jié)果顯示（以低劑量組為參照）：高劑量組（group1）與低劑量組（group3）的OR值為3.136，OR95%CI為（1.389，7.082）；中劑量組（group2），與低劑量組組（group3）的OR值為