全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
14284 6
2012-11-02
如何画OR值及其95%CI 图。
类似如下图:如何用SAS实现?
附件列表
QQ截图20121102222302.png

原图尺寸 59.49 KB

QQ截图20121102222302.png

二维码

扫码加我 拉你入群

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

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

全部回复
2012-11-2 22:38:05
这种图,sas一般没有预置的procedure或者相关选项,一般要根据自己的需要利用sas的annotation facility画出来,可能有人写过相关的程序,如果能找到,根据自己的需要在人家的程序上修改一下就行了
二维码

扫码加我 拉你入群

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

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

2012-11-3 03:07:51
第一种方法,在logistic过程里实现data ex;
input
age sex opt;
cards;
51.00        1.00        .00
57.00        1.00        .00
46.00        1.00        .00
20.00        1.00        1.00
50.00        .00        .00
22.00        1.00        .00
40.00        1.00        .00
29.00        .00        1.00
68.00        1.00        .00
66.00        .00        .00
28.00        1.00        1.00
43.00        .00        1.00
43.00        .00        .00
53.00        .00        1.00
69.00        1.00        .00
63.00        .00        1.00
47.00        .00        1.00
67.00        .00        .00
65.00        .00        .00
66.00        1.00        .00
24.00        .00        .00
38.00        .00        1.00
24.00        1.00        .00
40.00        1.00        1.00
33.00        1.00        1.00
36.00        1.00        1.00
68.00        1.00        .00
28.00        .00        1.00
43.00        .00        1.00
58.00        1.00        1.00
28.00        .00        1.00
27.00        .00        1.00
38.00        .00        .00
;
run;
ods graphics on;

proc logistic data=ex plots(only)=(oddsratio(type=vertical));
model opt=age sex;
oddsratio 'age' age;
oddsratio 'sex' sex;
run;


第二种方法 ,使用ods output oddsratios =oddsratio;
将OR值保存到数据集种,使用gplot实现,
data anno;                                                                                                                              
   length function style color $8;                                                                                                      
   retain xsys ysys '2' when 'a';                                                                                                      
   set oddsratio;                                                                                                                                                                                             
   function='move'; xsys='2'; ysys='2'; yc=effect; x=lowercl; color='black'; output;                                                  
   function='draw'; x=uppercl; color='black'; size=1;output;


   function='move';xsys='2'; ysys='2';yc=effect; x=lowercl; color='black'; output;                                                   
   function='draw';x=lowercl; ysys='9'; y=1; size=1; output;                                                                       
   function='draw';x=lowercl; y=-2; size=1;output;  

   function='move';xsys='2'; ysys='2'; yc=effect; x=uppercl; color='black'; output;                                                   
   function='draw';x=uppercl; ysys='9'; y=1; size=1; output;                                                                       
   function='draw';x=uppercl; y=-2; size=1; output;   
run;                                                                                                                                    



axis1 label=none                                                                                                                        
      minor=none                                                                                                                        
      offset=(5,5);                                                                                                                     

axis2  order=(0 to 15 by 2.5)
label=('Odds Ratio')                                                                                                              
      minor=none;                                                                                                                       

title "plot of odds ratio";
symbol1 i=none color=black value=dot height=1;                                                                                 

proc gplot data=oddsratio;                                                                                                                  
   plot effect*OddsRatioEst / annotate=anno                                                                                                      
        nolegend                                                                                                                        
        vaxis=axis1                                                                                                                     
        haxis=axis2;                                                                                                                                                                                                                                       
run;                                                                                                                                    
quit;            



二维码

扫码加我 拉你入群

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

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

2012-11-4 11:08:55
Imasasor 发表于 2012-11-2 22:38
这种图,sas一般没有预置的procedure或者相关选项,一般要根据自己的需要利用sas的annotation facility画出 ...
谢谢,我再看看。
二维码

扫码加我 拉你入群

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

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

2012-11-4 11:09:50
ziyenano 发表于 2012-11-3 03:07
第一种方法,在logistic过程里实现data ex;
input
age sex opt;
是的,有原始数据,做的logistic model。我再研究研究,谢谢。
二维码

扫码加我 拉你入群

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

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

2012-11-4 21:04:22
ziyenano 发表于 2012-11-3 03:07
第一种方法,在logistic过程里实现data ex;
input
age sex opt;
用第一种方法时,ERROR: Java 储存库中缺失文件。
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群