我希望在逻辑回归后将所有置信区间图放在一个图中,用于所有的层变量。例如,我的SAS代码是:
proc logistic data=data1;
model y = x;
strata cv1;
output out=out1 unknown1=x_beta1 unknown2=lowerbound unknown3=upperbound unknown4=strata_variable;
run;
我不知道我可以在输出语句中使用哪些变量名称(unknown1 unknown2 unknown3)。与sas支持页面一样,它表示“如果指定了STRATA语句,则只有PREDICTED =,DFBETAS =和H =选项可用”,这里是link。
我的情节陈述将是:
proc sgplot data=out1;
scatter y=strata_variable x=x_beta1 / xerrorlower=lowerbound xerrorupper=upperbound
markerattrs=(symbol=circlefilled size=9);
run;
此page中的第一个图表显示了我想要的内容。对不起,我不能插入任何情节,因为我的声誉不够高。
答案 0 :(得分:0)
我找到另一种方法来完成这个。我写了一个宏do循环来获取每个分层数据。然后添加
ods output OddsRatios=odds_temp;
获得估计和置信区间并将所有阶层合并在一起以形成我需要的情节。