研究生医学统计中logistic回归课件.ppt

1、Logistic Logistic regressionregression研究生医学统计中logistic回归1 第一节第一节.非条件非条件logisticlogistic回归回归第二节第二节.条件条件logisticlogistic回归回归第三节第三节.应用及其注意事项应用及其注意事项研究生医学统计中logistic回归2 医学研究中常碰到的可能取值仅有两个（即），如发病与未发病、阳性与阴性、死亡与生存、治愈与未治愈、暴露与未暴露等，显然这类资料不满足多元（重）回归的条件 什么情况下采用什么情况下采用LogisticLogistic回归回归研究生医学统计中logistic回归3 Brown

2、(1980)在术前检查了53例前列腺癌患者，拟用年龄(AGE)、酸性磷酸酯酶(ACID)两个连续型的变量，X射线(X_RAY)、术前探针活检病理分级(GRADE)、直肠指检肿瘤的大小与位置(STAGE)三个分类变量与手术探查结果变量NODES（1、0分别表示癌症淋巴结转移与未转移）建立淋巴结转移的预报模型。实例研究生医学统计中logistic回归4（一）53例接受手术的前列腺癌患者情况 研究生医学统计中logistic回归5（二）26例冠心病病人和28例对照进行病例对照研究 研究生医学统计中logistic回归626例冠心病病人和28例对照者进行病例对照研究 研究生医学统计中logistic回

3、归7一、logistic回归模型 研究生医学统计中logistic回归8概率预报模型概率预报模型)(11011011011011)(exp11)exp(1)exp(ppXXppppppeXXXXXX研究生医学统计中logistic回归9二、模型的 Logistic回归参数的估计通常采用最大似然法最大似然法(maximum likelihood，ML)。最大似然法的基本思想是先建立似然函数与对数似然函数，再通过使对数似然函数最大求解相应的参数值，所得到的估计值称为参数的最大似然估计值。研究生医学统计中logistic回归10 研究生医学统计中logistic回归11三、回归三、回归 研究生医学统

4、计中logistic回归12优势比及其可信区间优势比及其可信区间 研究生医学统计中logistic回归13标准化回归标准化回归用于评价各自变量对模型的贡献大小用于评价各自变量对模型的贡献大小研究生医学统计中logistic回归14SAS程序程序研究生医学统计中logistic回归15 The LOGISTIC Procedure The LOGISTIC ProcedureAnalysis of Maximum Likelihood Analysis of Maximum Likelihood EstimatesEstimates研究生医学统计中logistic回归16 预报模型预报模型)91

5、69.16443.08896.5(8181818111)9169.16443.08896.5(exp11)9169.16443.08896.5exp(1)9169.16443.08896.5exp(XXeXXXXXX研究生医学统计中logistic回归17 The LOGISTIC Procedure The LOGISTIC ProcedureAnalysis of Maximum Likelihood Analysis of Maximum Likelihood EstimatesEstimates研究生医学统计中logistic回归18 预报模型预报模型123451234512345ex

6、p(0.06182.04530.76141.56410.06930.0243)1exp(0.06182.04530.76141.56410.06930.0243)11exp(0.06182.04530.76141.56410.06930.0243)iXXXXXXXXXXXXXXX研究生医学统计中logistic回归19四、回归四、回归 X 当只有一个自变量时，以相应的预报概率 为纵轴，自变量 为横轴，可绘制出一条S形曲线。回归参数的正负符号与绝对值大小，分别决定了S形曲线的方向与形状研究生医学统计中logistic回归20优势比改变优势比改变exp(exp(j j)个单位个单位研究生医学统计中

7、logistic回归21研究生医学统计中logistic回归22五、整个回归模型五、整个回归模型 研究生医学统计中logistic回归23似然比检验（似然比检验（likelihood ratio test）研究生医学统计中logistic回归24ROC曲线模型评价曲线模型评价研究生医学统计中logistic回归25ROC曲线模型评价曲线模型评价图图16-2 Logistic16-2 Logistic回归预报能力的回归预报能力的ROCROC曲线曲线研究生医学统计中logistic回归26六、六、logistic逐步回归（变量筛选）逐步回归（变量筛选）MODEL语句加入选项“SELECTION=S

8、TEPWISE SLE=0.100.10 SLS=0.100.10;”常采用似然比检验：决定自变量是否引入或剔除。)ln2()ln2()ln(ln210102LLLLLR研究生医学统计中logistic回归27模型中有模型中有X5、X6、X8，看是否引入看是否引入X1模型含X5、X6、X8的模型的负二倍对数似然为：50.402模型含X1、X5、X6、X8的模型的负二倍对数似然为：46.224)ln2(0L)ln2(1LX1,10.00410.0,1,178.4224.40402.502引入变量PLR研究生医学统计中logistic回归28 研究生医学统计中logistic回归29 111010

9、111exp()()nmiijipijpi pjLXXXX研究生医学统计中logistic回归30 研究生医学统计中logistic回归31 研究生医学统计中logistic回归32 研究生医学统计中logistic回归33研究生医学统计中logistic回归34 第三节第三节 应用及其注意事项应用及其注意事项应变量应变量为（二项）分类的资料为（二项）分类的资料（预测、判别、危险因素分析等等）（预测、判别、危险因素分析等等）研究生医学统计中logistic回归35 注意事项注意事项1.分类自变量的哑变量编码 2.为了便于解释，对二项分类变量一般按0、1编码，一般以0表示阴性或较轻情况，而1表示

10、阳性或较严重情况。如果对二项分类变量按+1与-1编码，那么所得的 ，容易造成错误的解释。)2exp(OR研究生医学统计中logistic回归36西、中西、中三种疗法哑变量化 其它中西其它西012011XX原资料原资料姓名姓名性别性别年龄年龄 疗法疗法张山张山150中西中西李四李四120西西王五王五018中中刘六刘六070中中赵七赵七135中西中西孙八孙八029西西哑变量化哑变量化姓名姓名性别性别年龄年龄X1X2张山张山15001李四李四12010王五王五01800刘六刘六07000赵七赵七13501孙八孙八0291000100121中中西西疗法XX研究生医学统计中logistic回归37 注意

11、事项注意事项2.自变量的筛选 不同的筛选方法有时会产生不同的模型。实际工作中可同时采用这些方法，然后根据专业的可解释性、模型的节约性和资料采集的方便性等，决定采用何种方法的计算结果。研究生医学统计中logistic回归38 注意事项注意事项3.交互作用 交互作用的分析十分复杂，应根据临床意义与实际情况酌情使用。研究生医学统计中logistic回归39 注意事项注意事项4.多分类logistic回归 心理疾病分为精神分裂症、抑郁症、神经官能症等（名义变量名义变量nominal variables)；疗效评价分为无效、好转、显效、痊愈(有序变量有序变量ordinal variables)。参见第1

12、7章应变量研究生医学统计中logistic回归40SPSS软件计算Analyze Regression Binary LogisticDependent:yCovariates:x1 x8 Method:Forward WardSave Predicted Values Probabilities Group membershipOption CI for exp 95%Probability for Stepwise Entry:0.1 Removal 0.15研究生医学统计中logistic回归41DATA samp16_1;INPUT x_ray grade stage age acid

13、 nodes;CARDS;.;PROC LOGISTIC DESCENDING;MODEL nodes=x_ray grade stage age acid/RISKLIMITS;OUTPUT OUT=pred PROB=pred;PROC PRINT DATA=pred;RUN;研究生医学统计中logistic回归42 The SAS System 22:07 Monday,November 29,2005 1 The LOGISTIC Procedure Model Information Data Set WORK.SAMP16_1 Response Variable nodes Num

14、ber of Response Levels 2 Number of Observations 53 Model binary logit Optimization Technique Fishers scoring Response Profile Ordered Total Value nodes Frequency 1 1 20 2 0 33 Probability modeled is nodes=1.研究生医学统计中logistic回归43 Model Convergence Status Convergence criterion(GCONV=1E-8)satisfied.Mode

15、l Fit Statistics Intercept Intercept and Criterion Only Covariates AIC 72.252 60.126 SC 74.222 71.948 -2 Log L 70.252 48.126 Testing Global Null Hypothesis:BETA=0 Test Chi-Square DF Pr ChiSq Likelihood Ratio 22.1264 5 0.0005 Score 19.4514 5 0.0016 Wald 13.1406 5 0.0221研究生医学统计中logistic回归44 The SAS Sy

16、stem 22:07 Monday,November 29,2005 2 The LOGISTIC Procedure Analysis of Maximum Likelihood Estimates Standard Wald Parameter DF Estimate Error Chi-Square Pr ChiSq Intercept 1 0.0618 3.4599 0.0003 0.9857 x_ray 1 2.0453 0.8072 6.4208 0.0113 grade 1 0.7614 0.7708 0.9759 0.3232 stage 1 1.5641 0.7740 4.0

17、835 0.0433 age 1 -0.0693 0.0579 1.4320 0.2314 acid 1 0.0243 0.0132 3.4230 0.0643研究生医学统计中logistic回归45 The SAS System 22:07 Monday,November 29,2005 2 The LOGISTIC Procedure Analysis of Maximum Likelihood Estimates Standard Wald Parameter DF Estimate Error Chi-Square Pr ChiSq Intercept 1 0.0618 3.4599

18、0.0003 0.9857 x_ray 1 2.0453 0.8072 6.4208 0.0113 grade 1 0.7614 0.7708 0.9759 0.3232 stage 1 1.5641 0.7740 4.0835 0.0433 age 1 -0.0693 0.0579 1.4320 0.2314 acid 1 0.0243 0.0132 3.4230 0.0643 Odds Ratio Estimates Point 95%Wald Effect Estimate Confidence Limits x_ray 7.732 1.589 37.614 grade 2.141 0.

19、473 9.700 stage 4.778 1.048 21.783 age 0.933 0.833 1.045 acid 1.025 0.999 1.051研究生医学统计中logistic回归46 Association of Predicted Probabilities and Observed Responses Percent Concordant 84.5 Somers D 0.694 Percent Discordant 15.2 Gamma 0.696 Percent Tied 0.3 Tau-a 0.332 Pairs 660 c 0.847 Wald Confidence

20、Interval for Adjusted Odds Ratios Effect Unit Estimate 95%Confidence Limits x_ray 1.0000 7.732 1.589 37.614 grade 1.0000 2.141 0.473 9.700 stage 1.0000 4.778 1.048 21.783 age 1.0000 0.933 0.833 1.045 acid 1.0000 1.025 0.999 1.051研究生医学统计中logistic回归47 Obs no x_ray grade stage age acid nodes _LEVEL_ pred 1 1 0 1 1 64 40 0 1 0.25511 2 2 0 0 1 63 40 0 1 0.14633 3 3 1 0 0 65 46 0 1 0.21842 4 4 0 1 0 67 47 0 1 0.06459 。50 50 1 0 1 64 89 1 1 0.80302 51 51 0 1 0 59 99 1 1 0.29880 52 52 1 1 1 68 126 1 1 0.94215 53 53 1 0 0 61 136 1 1 0.76730研究生医学统计中logistic回归48