Logistic 里面如何实现分层?急！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

是这样的。我的logistic模型为

disease（0 =没病 ，1=有病）=1.5sex（0=女，1=男）+1.6somke（0=从不吸烟 1=戒烟  2=现在吸烟）+...+1.5area;

proc logistic data=pht2  descending;
  class    smoke   /param=ref ref=first ;
model disease= sex  smoke  。。。 area ;

run;

其中area的取值为0=城市，1=农村，上述模型的用意是探讨危险因素，可以解释为：以area为例：农村与城市相比，发生疾病的风险提高了50%。

那么现在问题是 我发现城市和农村的发病风险是不一样的，所以我要分地区探讨这些危险因素对疾病的影响，如果直接把两个地区的人群分开再运行一遍程序，会不会忽略了一些综合的影响？有没有什么分层语句可以应用，像是  strata 之类的，使得在方程中综合考虑了所有的危险因素后探讨危险因素在城市和农村分布的异同？也就是是 加一个分层的语句后使得原本出现1个OR值变为出现两个OR值？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：logistic ogistic logisti logist logis logistic 分层

logistic回归不能做分层，理论上行不通。人为按危险因素分层，会破坏原回归方程的贡献率，还会影响其他因因素危险度的估计。logistic回归可以配对，但不能分层。分层法只能直接用一些分层技术如Mantel-Haenszel。也可以建立在一些新方法上，如倾向指数(propensity score).

你可以针对城市和农村做两个模型

2# andy162639
也许我没有真正理解strata 的含义，不过proc logistic 的确有strata statement。即使你所提到的cmh test，proc logistic也有对应的code：

1. 
   
2. proc freq ;
   
3.    tables A*B*Y/ cmh;
   
4. run;
   
5. 
   
6. proc logistic;
   
7.      class B / param=ref;
   
8.       model Y=B;
   
9.       strata A;
   
10.   run;

复制代码