thank you, yes i think i have to do it like that... do you have also an idea how i can code 06.06.6666 and 07.07.7777 as special missing? ... View more

doesnt work 😞 ... View more

I found out by myself. juhuu 🙂 its because my variable is numeric not character and i replaced xc to x. Then everything worked. But still have a problem with my by-Variable. This is my code now: proc sql; create table try as select period, time,mean(Score) as mean format=6.2, mean(Score)-std(Score) as std_lower, mean(Score)+std(Score) as std_upper from final group by time; run; Proc sort data=try; by period; run; data test; set try; If std_lower LT 0 Then std_lower=0; xsys='2'; ysys='2'; when='b'; x=time; y=std_lower; function='move'; output; y=std_upper; function='draw'; output; xsys='2'; x=time; y=std_upper; function='move'; output; xsys='7'; x=-2; function='draw'; output; x=+4;function='draw'; output; If std_lower GT 0 Then do; xsys='2'; x=time; y=std_lower; function='move'; output; xsys='7'; x=-2; function='draw'; output; x=+4; function='draw'; output; end; run; axis1 LABEL=(angle=90 Score') Order=(0 1 2 3 4 5 6 7 8 9 10 11 12) minor=none; axis2 LABEL=('Time [min]') Order=(-60 to 360 by 60)offset=(,2); symbol1 value=none interpol=none color=black; symbol2 value=circle height=3pct interpol=none color=black; proc gplot data=try anno=test; plot std_lower*time=1 mean*time=2 std_upper*time=1 / overlay vaxis=axis1 haxis=axis2 hminor=2 vref=6 LVREF=20; run; quit; If i put a by-Statement to proc gplot i get for every period the same Graph. Do you have any idea how i can draw a Graph for each period with the corresponding data (without having an own dataset for each period)? Message was edited by: marieK ... View more

Hi Robert, thanks for the code. It is working on the example, but unfortunately not on my data. It draws the mean but not the std 😞 On my x-axis i have a time-variable (-60 to 360 min). On my y-axis i have my measured value (i.e. pulse (in my data i have a score from 0 to 12)). For every measure (at beginning every 20 min then every 40 min) i want mean and std of pulse/score about all my 100 patients for every period (i have 4 periods) (period is my by-variable). It worked well with: proc template; define statgraph errorplot; begingraph; layout overlay / xaxisopts=( Label="Time [min]" type=linear linearopts=(tickvaluelist=(-60 0 60 120 180 240 300 360) viewmin=-60 viewmax=360) ) yaxisopts=( Label="Score" type=linear linearopts=( tickvaluelist=( 0 1 2 3 4 5 6 7 8 9 10 11 12 ) viewmin=0 viewmax=12 ) ); scatterplot x=time y=MeanScore / yerrorlower=eval(MeanScore - STD) yerrorupper=eval(MeanScore + STD); seriesplot x=time_value y=MeanScore; ReferenceLine y=6 / clip=true Lineattrs=( Pattern=34); endlayout; endgraph; end; run; proc summary data=test nway; class time; var Score; by period ; output out=weights(drop=_type_ _freq_) mean=Mean std=STD; run; proc sgrender data=weights template=errorplot; By period; run; But unfortunately i cant customize my axis like i told above 😞 regards, marie ... View more

Hey, no i dont want a boxplot. i want a mean-and-error-plot (mean, mean+std, mean-std) 🙂 I also have the un-summarized data. Message was edited by: marieK ... View more

Hi Dan, thank you very much. It is working 🙂 ... View more

Hi Dan, thanks for your tip! I will try it the next days! ... View more

Thank you! The translate-function is good 🙂 ... View more

Hi, thank you for your tip. now its working 🙂 do you have any idea how i join the upper and lower std and mean with a vertical line? Like: http://web.anglia.ac.uk/numbers/biostatistics/mann_whitney/local_folder/C07fig02_ErrorPlot_Bones.jpg ... View more

thank you! But when i do a*b=by-variable then i get an overlay-Graph. And that i dont want 🙂 Any other idea? Message was edited by: marieK ... View more

Juhuuu 🙂 thanks! ... View more

mmhh i think the FRINGE plot is not a good solution for me. Do you maybe have an idea how i can do an errorplot (mean, mean+std, mean-std) in Gplot? I tried it with: proc sql; create table try as select time_value,mean(score) as mean format=6.2, mean(score)-std(score) as std from test group by time union select time_value,mean(score) as mean format=6.2, mean(score)+std(score) as std from test group by time; symbol2 color=black interpol=stdtj cv=black ci=black line=1 value=dot; axis1 label=('Time [min]')order=(-60 TO 360 BY 60); axis2 label=(angle=90 'score') order=(0 to 12 by 1) value=(color=black); proc gplot data=try; where std NE . and time NE . ; format mean std 3.0; plot (std mean)*time/ overlay haxis=axis1 vaxis=axis2; run; quit; but it doesnt work properly... Any other idea? Thanks! Marie ... View more