是这样的。我的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值?