9logistic回归分析.ppt

《9logistic回归分析.ppt》由会员分享，可在线阅读，更多相关《9logistic回归分析.ppt（23页珍藏版）》请在三一办公上搜索。

1、STATA软件的应用,-Logistic 回归分析,非常简便、快速,终二恶左鸵贴溅沃赡香仑富剁蛹炼歉镶蚀暇蓝蔼绣怯僻罐酝氦郝锦亚肘眯9logistic回归分析9logistic回归分析,因变量（结局）：分类变量二分类变量(二分类logistic回归)有序多分类变量(有序多分类logistic回归)无序多分类变量（无序多分类logistic回归）自变量（各种影响因素）：可以是分类变量，也可以是连续型变量。,数据特征,欧求乐毗葫彤菇亿风棵弯沈而莉隔邦舱豫趴沉蛔起铣俩酗宦旅讽清涛多夏9logistic回归分析9logistic回归分析,二分类资料的分析,非条件logistic模型：成组病例对照研究资

2、料条件logistic模型：配比病例对照研究资料,仗遁空尉终堂衔僳讽碴棘租两阴垄壕染犀觉枢迂你紧球傅腐赦奉萌俯腆漾9logistic回归分析9logistic回归分析,非条件logistic回归模型,烟晋颖焰嗓炔主揽类讹攘着临惧聊市囤阮虽游辗肛献节住俞橙孪螟募蝎几9logistic回归分析9logistic回归分析,参数的估计：极大似然（MLE）假设检验：似然比检验:G=-2lnL-(-2lnL)Wald检验:z 检验,参数估计与假设检验,挣翅理柴剩朝奠产镁蹦容踞京马荣汾谐怎寸恕稳耸互赎书嘴善谓排骄竟骇9logistic回归分析9logistic回归分析,回归系数的解释,回归系数 表示当其它自

3、变量固定不变时，X每改变一个单位，优势对数的改变量（优势比的对数）。,开傲甲茹相阅灯神虹腿命犁东借坟薯喳丰掐磺勘尤抹鲸歼宝旷稳斯烦柜牟9logistic回归分析9logistic回归分析,回归系数的解释,回归系数的解释二分类变量：OR=exp()表示暴露组发病的危险是非暴露组的几倍。等级变量：一般以最小等级或最大等级作为参考组，并按等级顺序依次取为0，1，2，。此时，OR=exp()表示X增加一个等级时，发病危险变为原来的几倍。连续性变量：表示增加1(个计量单位)时，发病危险变为原来的几倍。,断呐鲤钙挚泼纸态徘肄状冒巳忻砚来莹摧巷月向呻选矛儒内叫霹袄烃别坪9logistic回归分析9logis

4、tic回归分析,回归系数的解释,多分类变量：哑变量(dummy variable)x1时：x11，x20，x30，x40 表示A型血x2时：x10，x21，x30，x40 表示B型血x3时：x10，x20，x31，x40 表示AB型血x4时：x10，x20，x30，x41 表示O型血产生哑变量:tab x,gen(x),月抿法姐醚噎轻歼蔗潍缝颐瓢伤跌靳悸游嫩沼宽殿善娇挨怖乘唱矛迹麦钦9logistic回归分析9logistic回归分析,非条件logistic回归,logit 因变量 自变量,选择项sw logit 因变量 自变量,选择项 选择项：or 指定结果中给出OR值，缺失时输出回归系数

5、pr(#)是剔除变量的P值 pe(#)是选入变量的P值,节呸拼分衬翁锯耀殆谦恋肩炽疫拟颈匝厢问托律例恤下沼藕娩肢缄起疽蹦9logistic回归分析9logistic回归分析,例1(成组病例对照研究)某单位研究胸膜间皮瘤与接触石棉的关系，资料见下表。试对其进行分析。,方法1：卡方检验方法2：拟合logistic回归模型，即,白硼硬具釜漱峡盏障烃定顺漆魁郁跌骏址进妈夫财产键牙则笑根恿坏心暇9logistic回归分析9logistic回归分析,数据结构,寂翁酥杭的渝群涝妄斯溪树碟生俗注握历但遁著肮纤傅斑句坟冻田芥约钦9logistic回归分析9logistic回归分析,方法1：2 检验,cc cas

6、e exposure fw=f Proportion|Exposed Unexposed|Total Exposed-+-+-Cases|40 36|76 0.5263 Controls|9 67|76 0.1184-+-+-Total|49 103|152 0.3224|Point estimate|95%Conf.Interval|-+-Odds ratio|8.271605|3.4193 21.33091(exact)Attr.frac.ex.|.8791045|.7075425.9531197(exact)Attr.frac.pop|.4626866|+-chi2(1)=28.94 P

7、rchi2=0.0000,吏抓宜磊啸或橙脆蜕眨腹寻来咕仙锑绞姻冒局虹迎患臣踢跋咀疆岔赊俯雅9logistic回归分析9logistic回归分析,方法2：logistic回归输出回归系数,logit case exposure fw=fLogit estimates Number of obs=152 LR chi2(1)=30.67 Prob chi2=0.0000Log likelihood=-90.024994 Pseudo R2=0.1455-case|Coef.Std.Err.z P|z|95%Conf.Interval-+-exposure|2.112829.4228578 5.00

8、 0.000 1.284043 2.941615 _cons|-.6211737.2066474-3.01 0.003-1.026195-.2161522-似然比2=30.67，P0.0000,因此可以认为模型有意义。,咒如踞摇丙屑盎第酒翟只蠢癣浸筛猪崇亮秸券脉疡滚遮膏哑苟洁然粤襟娟9logistic回归分析9logistic回归分析,方法2：logistic回归输出OR,logit case exposure fw=f,orLogit estimates Number of obs=152 LR chi2(1)=30.67 Prob chi2=0.0000Log likelihood=-90

9、.024994 Pseudo R2=0.1455-case|Odds Ratio Std.Err.z P|z|95%Conf.Interval-+-exposure|8.271605 3.497713 5.00 0.000 3.611209 18.94641-胸膜间皮瘤与接触石棉的关系有统计学意义，OR=8.27（95%CI:3.61-18.95）,可认为接触石棉者发生间皮瘤的危险性为未接触者的8.27倍。,住父验红宵肠梁橱虫嘻勒乞篱挺晒胸痞腋提嗡坤裙导估服嘿磁渴晾学贝兑9logistic回归分析9logistic回归分析,分类资料的分析,非条件logistic模型：成组病例对照研究资料 队列

10、研究资料条件logistic模型：配比病例对照研究资料,癌涟根缓驴苯猖恭魏掠瀑粒戈谆些兑抒廊沙颐捐去姨结矮咙例惹嗜秩政侣9logistic回归分析9logistic回归分析,回归的常数项因同一层病例和对照的基线患病（发病）概率相同被抵消掉了，因此不能作预测，只能作因素分析。回归系数：表示病例与对照变量值之差与患病优势的关系，即exp()表示病例与对照暴露水平相差一个单位时患病的优势比。,条件logistic回归模型,缎擅穷琳咯俄酌舵两烩纱瘁楚蔚顾滥桓勇旋吵缩凭许拦旱宗涕铁窥毁驳践9logistic回归分析9logistic回归分析,例3 在子宫内膜癌与使用雌激素关系的研究中，运用了1:4的病例

11、-对照研究，配比因素为年龄，共调查了63对，315例。各变量定义如下：,辐蟹酗毒蔚搪圣童薪俩引羊宏芭减逞怪摄谦谋隅企玩潦蜡腐料舜锄郑谚膀9logistic回归分析9logistic回归分析,数据结构,党椽愧倾瞎斧傈狄菊桩锁茬由溶迪浆外陀瓤澳窍坎恢犀惺撕像烦得羹窥苹9logistic回归分析9logistic回归分析,条件logistic回归,clogit 因变量 自变量,group(配比变量)or group()是必选项，它是用来区分各配比组的。or的意义同logit。,渊句肛裕刨撇企戮上锈笑紧脑栏矗棕诈捍抚佣呸捅箕筛增锦凋埋臀膏愤孕9logistic回归分析9logistic回归分析,clo

12、git y ht est drug,group(match)or,Conditional(fixed-effects)logistic regression Number of obs=315 LR chi2(3)=39.45 Prob chi2=0.0000Log likelihood=-81.671759 Pseudo R2=0.1945-y|Odds Ratio Std.Err.z P|z|95%Conf.Interval-+-ht|.8954587.3026653-0.33 0.744.4616781 1.736808 est|6.903735 2.997048 4.45 0.000

13、2.948225 16.16619 drug|2.142038 1.05669 1.54 0.123.8145551 5.632922-,腮禾雀醋茸栋奖纫盈匠栓忙汕娩业瘪投爬统赊虚恢私耪余膝戒延鲁汛闻棠9logistic回归分析9logistic回归分析,sw clogit y ht est drug,group(match)or pe(0.10)pr(0.11),sw clogit y ht est drug,group(match)or pe(0.10)pr(0.11)begin with full modelp=0.7439=0.1100 removing htp=0.1314=0.11

14、00 removing drugConditional(fixed-effects)logistic regression Number of obs=315 LR chi2(1)=36.87 Prob chi2=0.0000Log likelihood=-82.959065 Pseudo R2=0.1818-y|Odds Ratio Std.Err.z P|z|95%Conf.Interval-+-est|8.175451 3.429501 5.01 0.000 3.592884 18.60289-,班伞甭烂孙谐芬计句拖瞻卜漏歌状忙矾加蜘崖小惧遂瀑淄厕憾勿扰弃遁咎9logistic回归分析9

15、logistic回归分析,结论,子宫内膜癌与使用雌激素关系有统计学意义，OR=8.18（95%CI:3.59-18.60）,可认为使用雌激素者发生子宫内膜癌的危险性为未使用者的8.18倍。,啪联珊香尼况削究吟巡笼泞俊右旧忧谗芭疑遮莹泡槛斡研金埔照锨疼抗僵9logistic回归分析9logistic回归分析,课后练习,利用帮助功能学习STATA软件中的Logistic回归分析help logit 二分类（非条件）help clogit 二分类（条件）help ologit 有序多分类help mlogit 无序多分类,淆羊育率椎济卖送估笋长矣滞忍探港退牙掺神蝉铱缸臭潭睛莆卒稻尤峻题9logistic回归分析9logistic回归分析,