如何用SAS画森林图

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

此处是同事倾情奉献的程序，用于画森林图，做出的图如下（宝宝一向觉得template画图选项太多，没sgplot简单）

程序如下：

1. data _forest;
   
2. input group1 $1-34  group2 $35.OR LowerCL UpperCL ;
   
3. datalines;
   
4. Age                                                                  0.999 0.992 1.006
   
5. BMI                                                                  0.953 0.929 0.976
   
6. Gender                            (Female)                           0.640 0.526 0.777
   
7. Gender                            Male (control)                     .    .      .   
   
8. Diabetes                          Yes                                1.120 0.956 1.311
   
9. Diabetes                          No (control)                       .     .     .         
   
10. Hypertension                      Yes                                1.079 0.903 1.294
    
11. Hypertension                      No (control)                       .     .     .         
    
12. Smoking status                    (Past smoker)                      0.935 0.768 1.138
    
13. Smoking status                    (Current smoker)                   0.714 0.571 0.892
    
14. Smoking status                    Non smoker (control)               .     .     .         
    
15. Dyslipidemia                      Yes                                0.577 0.497 0.668
    
16. Dyslipidemia                      No (control)                       .     .     .
    
17. lipid lowering therapy            Yes                                3.540 2.879 4.388
    
18. lipid lowering therapy            No (control)                             .      .     .         
    
19. ;
    
20. run;
    
21. 
    
22. data _forest2;
    
23. length  name1 $70.;
    
24. set _forest;
    
25. if group2^="Yes" then name1=strip(group1)||" "||strip(group2);
    
26.         else name1=strip(group1);
    
27. if or=.   then delete;
    
28. or1='OR(95%CI)';
    
29. CI=right(strip(PUT(OR,best.))||"("||strip(PUT(LowerCL,best.))||", "||strip(PUT(UpperCL,best.))||") ");
    
30. run;
    
31. 
    
32. *用于加annotation，出来的效果要自己酌情调整;
    
33. data _f9_anno;
    
34.    length function $4 color style $20 label $40;
    
35.         retain width 40;
    
36.    function="text";
    
37.    color="black"; style="'Times New Roman'";
    
38.    xsys='3'; x1=58;
    
39.    ysys='3'; y1=98;
    
40.    label='Hazard Ratio and 95% CI';
    
41.    output;
    
42. run;
    
43. 
    
44. proc sgplot data=_forest2 sganno=_f9_anno noautolegend;
    
45.   scatter y=name1 x=or / xerrorupper=UpperCL xerrorlower=LowerCL
    
46.             markerattrs=graphdata1(symbol=squarefilled size=3);
    
47.   scatter y=name1 x=or1 / markerchar=CI x2axis ;
    
48. /*  scatter y=name1 x=or2 / markerchar=group1  ;*/
    
49.   refline 1 / axis=x lineattrs=(pattern=shortdash) transparency=0.5;
    
50.   inset 'Goal unachieved'  / position=bottomleft;
    
51.   inset 'Goal achieved'  / position=bottom;
    
52.   xaxis type=log offsetmin=0 offsetmax=0.35 min=0.1 max=5
    
53.           minor display=(nolabel) ;
    
54.   x2axis offsetmin=0.7 offsetmax=0.01 display=(noticks nolabel);
    
55.   yaxis display=(noticks nolabel) offsetmin=0.071 offsetmax=0.035;
    
56. run;

复制代码

xerrorupper& xerrorlower 功能同yerrorupper和yerrorlower，只是一个是Y轴的上下限，一个是X轴的，symbol决定图标样式，size决定图标大小。markerchar和x2axis用于添加次级X轴的值。Refline用于指定中心，此图以1为中心，画虚线分隔，设置虚线透明度。Inset语句就是添加类似label的语句，可以调节位置。

此处几个offset语句，是用于调节X轴和Y轴的位置的，offsetmin，离左边的距离，offsetmax，离右边的距离，0.35,0.1都是距离的百分比。

扫码加我 拉你入群

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

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

分享0 收藏7 回帖

关键词：如何用 森林图 Transparency Diabetes function 森林图

谢谢楼主分享

谢谢楼主分享

如果我想再加一列P值，程序该如何改动，多谢大神！

另外，左侧的分类条目怎么做成总类下面分条目的形式？

请问如何更改markerchar的字号

想问下 置信区间两边的小竖线 是什么选项 我做出来的就没有竖线

楼主没看到代码啊

你这不是森林图吧，你这只是简单的罗列了各个指标的OR(95%CI)，并没有包含合并的OR(95%CI)和p值呀