如何画OR值及其95%CI 图。

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

如何画OR值及其95%CI 图。
类似如下图：如何用SAS实现？

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：如何用 如何

这种图，sas一般没有预置的procedure或者相关选项，一般要根据自己的需要利用sas的annotation facility画出来，可能有人写过相关的程序，如果能找到，根据自己的需要在人家的程序上修改一下就行了

第一种方法，在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;            

Imasasor 发表于 2012-11-2 22:38
这种图，sas一般没有预置的procedure或者相关选项，一般要根据自己的需要利用sas的annotation facility画出 ...

谢谢，我再看看。

ziyenano 发表于 2012-11-3 03:07
第一种方法，在logistic过程里实现data ex;
input
age sex opt;

是的，有原始数据，做的logistic model。我再研究研究，谢谢。

ziyenano 发表于 2012-11-3 03:07
第一种方法，在logistic过程里实现data ex;
input
age sex opt;

用第一种方法时，ERROR: Java 储存库中缺失文件。

webgu 发表于 2012-11-4 21:04
用第一种方法时，ERROR: Java 储存库中缺失文件。

新的作图模式是调用java完成的；
应该是装SAS的时候java文件没有配置好，
去support.sas.com找找解决办法，我这边能够作图。