（SAS作点线图加上标准差）怎样在SAS中做点线图并加上标准差

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

如上图所示，
A组数据是一列观察值，用方块表示；

B组数据是每一行是3次重复的平均值，用三角表示，该组数据同时给出了标准差。
问题是：怎样在SAS中做如上图所示的点线图并给B组数据加上标准差？

数据用例如下：

A	B	B的标准差
33.99	33.95	1.5
29.54	29.73	2.5
25.80	25.91	3.6
29.45	29.50	5
21.92	22.11	4
25.10	25.38	0.99
23.53	23.61	1
26.92	27.35	2
25.53	25.81	4
28.04	28.09	7

扫码加我 拉你入群

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

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

分享0 收藏5 回帖

关键词：标准差 ABB 平均值 标准差 平均值

1. PROC DATASETS LIBRARY=WORK KILL NOLIST;QUIT;
   
2. 
   
3. DATA DATAIN;
   
4. X=_N_;
   
5.         INPUT A B BSTD;
   
6. DATALINES;
   
7. 33.99 33.95 1.5
   
8. 29.54 29.73 2.5
   
9. 25.80 25.91 3.6
   
10. 29.45 29.50 5
    
11. 21.92 22.11 4
    
12. 25.10 25.38 0.99
    
13. 23.53 23.61 1
    
14. 26.92 27.35 2
    
15. 25.53 25.81 4
    
16. 28.04 28.09 7
    
17. ;
    
18. RUN;
    
19. 
    
20. DATA DATAUSE;
    
21.         SET DATAIN(RENAME=(A=VALUE) DROP=B BSTD);
    
22.         STD=0;GROUP="A";OUTPUT;
    
23.         SET DATAIN(RENAME=(B=VALUE BSTD=STD) DROP=A);
    
24.         GROUP="B";OUTPUT;
    
25. RUN;
    
26. DATA DATAUSE;
    
27.         SET DATAUSE;
    
28.         YLOW=VALUE-STD;
    
29.         YUP=VALUE+STD;
    
30. RUN;
    
31. 
    
32. 
    
33. 
    
34. *Template Codes;
    
35. proc template;
    
36. define statgraph barChart;
    
37. begingraph/backgroundcolor=white border=False;
    
38.         EntryTitle "LineChart with ErrorBars showing (*ESC*){unicode '00B1' x} SE"/textattrs=(size=10pt) pad=(bottom=10);
    
39.         layout LATTICE /columns=1;
    
40.                    layout OVERLAY/
    
41.                 Xaxisopts=(gridDisplay=ON LABEL="Observation" LABELATTRS=(weight=bold color=black)
    
42.                                         LINEAROPTS=(TICKVALUEFITPOLICY=rotate TICKVALUEPRIORITY=TRUE))
    
43.                 Yaxisopts=(gridDisplay=ON LABEL="Mean" LABELATTRS=(weight=bold color=black));
    
44.                            ReferenceLine y=30;
    
45.                         *ErrorBar;
    
46.                         Scatterplot X=X  Y=VALUE/ primary=true index=trtnum yerrorlower=ylow yerrorupper=yup
    
47.                 name="scatter"  
    
48.                                 GROUP=group /***Group A&B here in a new colume**/
    
49.                          errorbarattrs=(thickness=2) markerattrs=(size=8 weight=bold);
    
50.                         *Series;
    
51.                            SERIESPLOT         X=X  Y=VALUE/primary=true name="Series1" group=group;
    
52.                   endlayout;
    
53.         *DiscreteLegend;
    
54.         sidebar / align=bottom spacefill=false;
    
55.                        discretelegend "Series1" / across=2 valueattrs=(size=10pt) pad=10px border=False;
    
56.            endsidebar;
    
57.         endlayout;
    
58. Endgraph; end; run;
    
59. 
    
60. *Call template;
    
61. proc sgrender data=DATAUSE template=barChart; run;

复制代码

Damn *****! No mind on work!

Hi, xiaodanzi，多谢你，太厉害了！

yukiooy 发表于 2013-1-29 15:12
Hi, xiaodanzi，多谢你，太厉害了！

{:2_25:}

真是厉害啊

用gplot试试：
DATA ex;
x=_n_;
INPUT A B B_STD;
y=b-b_std;
output;
y=b+b_std;
output;
DATALINES;
33.99 33.95 1.5
29.54 29.73 2.5
25.80 25.91 3.6
29.45 29.50 5
21.92 22.11 4
25.10 25.38 0.99
23.53 23.61 1
26.92 27.35 2
25.53 25.81 4
28.04 28.09 7
;
RUN;

symbol1 i=join v=square c=blue;
symbol2 i=join v=triangle c=red;
symbol3 i=hiloctj c=green;
proc gplot data=ex;
plot a*x=1 b*x=2 y*x=3/overlay;
run;
quit;

Please don't make SAS even more complicated! Jingju

1. DATA ex;
   
2. INPUT A B B_STD;
   
3. t ++1;
   
4. b_l =b-b_std;
   
5. B_u =b+b_std;
   
6. DATALINES;
   
7. 33.99 33.95 1.5
   
8. ...
   
9. ;
   
10. proc sgplot noautolegend;
    
11. scatter x =t y =b/yErrorLower =b_l yErrorUpper =b_u markerattrs =(size =0);
    
12. series x =t y =b/markers markerattrs =(symbol =squarefilled) name ="sb";
    
13. series x =t y =a/ lineAttrs =(color =blue) markers markerattrs =(symbol =TriangleFilled) name ="sa";
    
14. keyLegend "sa" "sb"/location =inside;
    
15. xaxis values =(1 to 10 by 1) ;
    
16. yaxis grid display =(noLabel);
    
17. run;

复制代码

看看··借用··

有道理！

非常棒的答案