Stata做logistic回归89003.ppt
《Stata做logistic回归89003.ppt》由会员分享,可在线阅读,更多相关《Stata做logistic回归89003.ppt(36页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、应用应用Stata做做logistic回归回归何保昌何保昌 实际生活中经常会遇到因变量只有实际生活中经常会遇到因变量只有0和和1的二分类变量,的二分类变量,不能满足正态性和方差齐性,故不能直接使用线性模型不能满足正态性和方差齐性,故不能直接使用线性模型来拟合方程。来拟合方程。Logistic回归正是处理因变量是二分类或多分类变量的回归正是处理因变量是二分类或多分类变量的一种方法。现已广泛应用于队列研究,病例对照研究和一种方法。现已广泛应用于队列研究,病例对照研究和试验性研究,成为分类因变量的首选多变量分析模型。试验性研究,成为分类因变量的首选多变量分析模型。分类分类 按因变量性质,可分为二分类
2、、无序多按因变量性质,可分为二分类、无序多分类、有序多分类。分类、有序多分类。按是否匹配可分为非条件和条件按是否匹配可分为非条件和条件logistic回归。回归。Stata软件专门有一组命令用于做不同软件专门有一组命令用于做不同类型的类型的logistic回归,例如:回归,例如:logit、blogit、glogit、clogit、mlogit、ologit。1.Logistic 回归命令:命令:logit 因变量因变量 自变量自变量,选择项选择项在进行在进行logistic 回归时要注意资料的形式。通常,回归时要注意资料的形式。通常,用于用于logistic 回归的资料有三种形式:回归的资料
3、有三种形式:(1)分水平频数资料,一般自变量较少,且均为分类分水平频数资料,一般自变量较少,且均为分类变量,常以各变量变量,常以各变量(包括因变量、自变量包括因变量、自变量)各水平的组各水平的组合的频数表形式出现。如例合的频数表形式出现。如例1。拟合时仍用上述命令,。拟合时仍用上述命令,只是命令中增加只是命令中增加fw=频数变量频数变量选择项。选择项。(2)分组频数资料,一般自变量较少,且均为分组频数资料,一般自变量较少,且均为分类变量,常以各自变量分类变量,常以各自变量(不包括因变量不包括因变量)各各水平的组合的频数表形式出现,因变量常表水平的组合的频数表形式出现,因变量常表达为分子与分母。
4、如例达为分子与分母。如例2。用下列命令:。用下列命令:blogit 阳性数变量阳性数变量 总观察数变量总观察数变量,logit 命命令选择项令选择项或或 glogit 阳性数变量阳性数变量 总观察数变量总观察数变量,level(#)or(3)个体水平资料,即一个观察对象一条记录。个体水平资料,即一个观察对象一条记录。如例如例3,直接使用,直接使用logit 命令估计即可。命令估计即可。拟合模型后可以用指令拟合模型后可以用指令predict 得到预测概率,得到预测概率,然后进行模型诊断、应用等。然后进行模型诊断、应用等。例例1 本例是探讨妇女使用雌激素与患子本例是探讨妇女使用雌激素与患子宫内膜癌
5、之间关系的病例宫内膜癌之间关系的病例-对照研究对照研究资料,见表资料,见表1,请计算,请计算OR 及其及其95可信区间。再可信区间。再用用logistic 回归估计参数,写出回回归估计参数,写出回归方程,并说明回归系数与归方程,并说明回归系数与OR 的关的关系。系。也可以用也可以用logit命令命令可以利用例可以利用例3做逐步回归做逐步回归Stata 用于逐步回归分析的命令是在要执行的命用于逐步回归分析的命令是在要执行的命令前增加令前增加sw。osw 回归命令回归命令 因变量因变量自变量自变量,筛选变量的筛选变量的P 值值 选择项选择项o其中,筛选变量的其中,筛选变量的P 值有值有3 种组合种
6、组合opr(#)/*后退法后退法ope(#)/*向前法向前法opr(#)pe(#)/*逐步后退法逐步后退法opr(#)pe(#)forward/*逐步向前法逐步向前法 pr(#)是剔除变量的P 值,pe(#)是选入变量的P 值,如果只选pr(#),则表示用后退法,如果同时选用pr(#)和pe(#)表示逐步法。应用时,为防止计算进入死循环,pr(#)须略大于pe(#)。例如,pe(0.05),pr(0.051)。1.用逐步后退法,剔选变量的概率为:用逐步后退法,剔选变量的概率为:pe(0.05),pr(0.06),结果如,结果如下:下:sw logit y x1 x2 x3,pe(0.05)pr
7、(0.06)2.用逐步前进法,剔选变量的概率不用逐步前进法,剔选变量的概率不变变,结果如下:结果如下:sw logit y x1 x2 x3,pr(0.06)pe(0.05)forward条件条件logistic 回归回归非条件非条件logistic 回归适用于平行组设计的病例回归适用于平行组设计的病例-对对照研究,队列研究,而不适用于配比设计的病例照研究,队列研究,而不适用于配比设计的病例-对照研究。对于配比的病例对照研究。对于配比的病例-对照研究资料需要对照研究资料需要用条件用条件logistic 回归。其命令为:回归。其命令为:oclogit 因变量因变量 自变量自变量,group(配比
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Stata logistic 回归 89003
限制150内