logistic回歸分析課件

STATA軟件的應用---Logistic回歸分析非常簡便、快速因變量（結局）：分類變量二分類變量(二分類logistic回歸)有序多分類變量(有序多分類logistic回歸)無序多分類變量（無序多分類logistic回歸）自變量（各種影響因素）：可以是分類變量，也可以是連續(xù)型變量。數據特征二分類資料的分析非條件logistic模型：成組病例對照研究資料

條件logistic模型：配比病例對照研究資料參數的估計：極大似然（MLE）假設檢驗：似然比檢驗:G=-2lnL-(-2lnL’)

Wald檢驗:

z檢驗參數估計與假設檢驗回歸系數的解釋回歸系數表示當其它自變量固定不變時，X每改變一個單位，優(yōu)勢對數的改變量（優(yōu)勢比的對數）。回歸系數的解釋回歸系數的解釋二分類變量：OR=exp()表示暴露組發(fā)病的危險是非暴露組的幾倍。等級變量：一般以最小等級或最大等級作為參考組，并按等級順序依次取為0，1，2，…。此時，OR=exp()表示X增加一個等級時，發(fā)病危險變?yōu)樵瓉淼膸妆?。連續(xù)性變量：表示增加1(個計量單位)時，發(fā)病危險變?yōu)樵瓉淼膸妆??；貧w系數的解釋多分類變量：啞變量(dummyvariable)x＝1時：x1＝1，x2＝0，x3＝0，x4＝0表示A型血x＝2時：x1＝0，x2＝1，x3＝0，x4＝0表示B型血x＝3時：x1＝0，x2＝0，x3＝1，x4＝0表示AB型血x＝4時：x1＝0，x2＝0，x3＝0，x4＝1表示O型血產生啞變量:

tabx,gen(x)例1(成組病例對照研究)某單位研究胸膜間皮瘤與接觸石棉的關系，資料見下表。試對其進行分析。

方法1：卡方檢驗方法2：擬合logistic回歸模型，即數據結構方法2：logistic回歸—輸出回歸系數logitcaseexposure[fw=f]LogitestimatesNumberofobs=152

LRchi2(1)=30.67Prob>chi2=0.0000Loglikelihood=-90.024994PseudoR2=0.1455------------------------------------------------------------------------------case|Coef.Std.Err.zP>|z|[95%Conf.Interval]-------------+----------------------------------------------------------------exposure|2.112829.42285785.000.0001.2840432.941615_cons|-.6211737.2066474-3.010.003-1.026195-.2161522------------------------------------------------------------------------------似然比2=30.67，P＝0.0000,因此可以認為模型有意義。方法2：logistic回歸—輸出OR

logitcaseexposure[fw=f],orLogitestimatesNumberofobs=152LRchi2(1)=30.67Prob>chi2=0.0000Loglikelihood=-90.024994PseudoR2=0.1455------------------------------------------------------------------------------case|OddsRatioStd.Err.zP>|z|[95%Conf.Interval]-------------+----------------------------------------------------------------exposure|8.2716053.4977135.000.0003.61120918.94641------------------------------------------------------------------------------胸膜間皮瘤與接觸石棉的關系有統計學意義，

OR=8.27（95%CI:3.61-18.95）,可認為接觸石棉者發(fā)生間皮瘤的危險性為未接觸者的8.27倍。分類資料的分析非條件logistic模型：成組病例對照研究資料隊列研究資料條件logistic模型：配比病例對照研究資料例3在子宮內膜癌與使用雌激素關系的研究中，運用了1:4的病例-對照研究，配比因素為年齡，共調查了63對，315例。各變量定義如下：數據結構條件logistic回歸clogit因變量自變量,group(配比變量)[or]

group()是必選項，它是用來區(qū)分各配比組的。or的意義同logit。swclogityhtestdrug,group(match)orpe(0.10)pr(0.11)

swclogityhtestdrug,group(match)orpe(0.10)pr(0.11)beginwithfullmodelp=0.7439>=0.1100removinghtp=0.1314>=0.1100removingdrugConditional(fixed-effects)logisticregressionNumberofobs=315LRchi2(1)=36.87Prob>chi2=0.0000Loglikelihood=-82.959065PseudoR2=0.1818------------------------------------------------------------------------------y|OddsRatioStd.Err.zP>|z|[95%Conf.Interval]-------------+----------------------------------------------------------------est|8.1754513.4295015.010.0003.59288418.60289------------------------------------------------------------------------------結論子宮內膜癌