第04讲logistic回归.pptVIP专享

第4讲logistic回归2南京医科大学流行病与卫生统计学系logistic的主要内容logistic回归条件logistic回归多类结果的logistic回归有序结果的累计比数logistic回归logistic回归的正确应用3南京医科大学流行病与卫生统计学系logistic回归模型简介传统方法与logistic回归的关系回归系数的解释模型的估计与假设检验建模策略4南京医科大学流行病与卫生统计学系多元线性回归的适用条件计量资料(均数)＝计量、等级、分类二分类资料(率)＝计量、等级、分类？5南京医科大学流行病与卫生统计学系1模型简介P发病的概率,0≤P≤1。1122mmPxxx1122()11mmxxxPeCox(1970)logit变换6南京医科大学流行病与卫生统计学系logistic函数11zPe1122mmzxxxz-5-4-3-2-10123450.2.4.6.81P7南京医科大学流行病与卫生统计学系logistic回归模型P发病的概率。为“机会”或“优势”(odds)logitP为优势之对数(logodds)mmxxxP2211logitPP1PPP1lnˆlogit8南京医科大学流行病与卫生统计学系logistic回归模型的几种形式mmxxxP2211logit1122()11mmxxxPe112211221mmmmxxxxxxePe9南京医科大学流行病与卫生统计学系2传统方法与logistic回归四格表资料(病例对照)与logistic的关系分层四格表资料与logistic的关系10南京医科大学流行病与卫生统计学系例四格表资料探讨妇女使用雌激素与患子宫内膜癌的病例-对照研究雌激素使用过未使用过病例55(a)128(b)对照19(c)164(d)7089.31281916455)1()1(0011bcadPPPPOR5616.6~0964.2164119112815511.96exp3.7089:%95CI11南京医科大学流行病与卫生统计学系四格表资料的logistic回归X=1表示使用过雌激素X=0表示未使用过雌激素xP1.3107+0.2478-=logitxPP1.3107+0.2478-=-1ln12南京医科大学流行病与卫生统计学系四格表资料与logistic回归X=1时X=0时3107.11ln1ln)1/()1/(ln)ln(00110011PPPPPPPPOR11ln=-0.24781.31071-PP00ln=-0.24781-PP1100(1)(1)PPORPP1.31073.7089ORe13南京医科大学流行病与卫生统计学系例分层四格表资料与logistic回归年龄组X=1,使用过X=0,未使用过Y=1Y=0Y=1Y=01:25~1946222242:30~34933123903:35~39426333304:40~4469653625:45~496593301按年龄分层的心肌梗死与近期口服避孕药的关系14南京医科大学流行病与卫生统计学系Mantel-Haenszel法结果.mhoddsyx[fw=f],by(age)Comparingx==1vs.x==0byage---------------------------------------------------------------age|OddsRatiochi2(1)P>chi2[95%Conf.Interval]------+--------------------------------------------------------1|7.2258066.780.00921.26334041.328772|8.86363628.640.00003.36912823.318813|1.5384620.580.44500.5049524.6873034|3.7128216.580.01031.26658710.883615|3.8838715.530.01871.14749013.14561---------------------------------------------------------------Mantel-Haenszelestimatecontrollingforage----------------------------------------------------------------OddsRatiochi2(1)P>chi2[95%Conf.Interval]----------------------------------------------------------------3.96989534.720.00002.4180416.517702----------------------------------------------------------------TestofhomogeneityofORs(approx):chi2(4)=6.27Pr>chi2=0.179715南京医科大学流行病与卫生统计学系logistic回归结果.xi:logityxi.age[fw=f]i.age_Iage_1-5(naturallycoded;_Iage_1omitted)LogitestimatesNumberofobs=1976LRchi2(5)=151.47Prob>chi2=0.0000Loglikelihood=-643.06749PseudoR2=0.1054----------------------------------------------------------------y|Coef.Std.Err.zP>|z|95%CI--------+-...