topic GPLOT Plot mean and confidence limits with error bars in Graphics Programming

GPLOT Plot mean and confidence limits with error bars

pakiprinsesse — Sun, 13 Jan 2019 16:27:28 GMT

I want to plot mean and confidence intervals in gplot. I have tried modifying a code found in the community (see below) for my data but it only produces the Means with a joined line and no error bars and I get the error message:

NOTE: ERROR DETECTED IN ANNOTATE= DATASET WORK.ANNO.

VARIABLE SHOWN IS NOT OF THE PROPER DATATYPE YC

NOTE: ERROR LIMIT REACHED IN ANNOTATE PROCESS. PROCESSING IS TERMINATED.

NOTE: PROCESSING TERMINATED BY INDIVIDUAL ERROR COUNT.

NOTE: 1 TOTAL ERRORS.

Could anyone please tell me how to solve the issue?

Thanks,

Anna

proc summary data=PAwPDAPFO;

by dage;

var CO;

output out=means (drop=_:) mean=mean lclm=lclm uclm=uclm;

run;

data anno;

length function color $8;

retain when 'a';

set means;

/*draw horisontal line from limit to limit*/

function='move'; yc=var; xsys='2'; ysys='2'; x=lclm; color='black'; output;

function='draw'; x=uclm; size=1; output;

/*draw tick line for lower*/

function='move'; xsys='2'; ysys='2'; yc=var; y=.; x=lclm; color='black'; output;

function='draw'; x=lclm; y=+1; ysys='9'; size=1; output;

function='draw'; x=lclm; y=-2; size=1; output;

/*draw tick line for upper*/

function='move'; xsys='2'; ysys='2'; yc=var; y=.; x=uclm; color='black'; output;

function='draw'; x=uclm; y=+1; ysys='9'; size=1; output;

function='draw'; x=uclm; y=-2; size=1; output;

run;

symbol1 value=dot height=1 color=black;

proc gplot data=means;

plot mean*dage/annotate=anno nolegend;

run; quit;

Re: GPLOT Plot mean and confidence limits with error bars

PaigeMiller — Sun, 13 Jan 2019 17:36:12 GMT

According to the error message:

VARIABLE SHOWN IS NOT OF THE PROPER DATATYPE YC

You have the wrong variable type for YC. It is expecting a character variable, and you have provided a numeric variable. But if you are trying to draw confidence intervals, you probably want to use variable Y and not YC.

Re: GPLOT Plot mean and confidence limits with error bars

DanH_sas — Sun, 13 Jan 2019 21:34:39 GMT

Try this instead. It's a lot less code 🙂

proc sgplot data=PAwPDAPFO noautolegend;
vline dage / response=CO stat=mean limitstat=clm markers;
run;

Hope this helps!

Dan

Re: GPLOT Plot mean and confidence limits with error bars

pakiprinsesse — Mon, 14 Jan 2019 07:57:33 GMT

Thank you!.

That worked as well.

In the mean time I succeeded with plot and the following code:

procsortdata=PAwPDAPFO;

bydage;

run;

procsummarydata=PAwPDAPFO;

bydage;

varCO;

outputout=means (drop=_:) mean=meanlclm=lclmuclm=uclm;

run;

/*CONVERT DAGE FROM NUMERIC TO CHARACTER*/

datameansCO; setmeans;

char_dage=put(dage,$2.);

dropdage;

renamechar_dage=dage;

run;

/*GRAPH*/

dataanno;

lengthfunction color $8;

retain/*xsys ysys '2' */when 'a';

setmeansco;

/*draw VERTICAL line from limit to limit (=ERROR BARS)*/

function='move'; Xc=dage; xsys='2'; ysys='2'; Y=lclm; color='black'; output;

function='draw'; Y=uclm; size=1; output;

/*draw tick line for lower CONFIDENCE LIMIT*/

function='move'; xsys='2'; ysys='2'; Xc=dage; X=.; Y=lclm; color='black';output;

function='draw'; Y=lclm; X=+1; Xsys='9'; size=1;output;

function='draw'; Y=lclm; X=-2; size=1; output;

/*draw tick line for upper CONFIDENCE LIMIT*/

function='move'; xsys='2'; ysys='2'; Xc=dage; X=.; Y=uclm; color='black';output;

function='draw'; Y=uclm; X=+1; Xsys='9'; size=1;output;

function='draw'; Y=uclm; X=-2; size=1; output;

run;

/*SET AXIS AND SYMBOLS*/

symbol1value=dot height=1color=black;

axis1order=0to28by1;

axis2order=0.4to1.75by0.1;

/*PLOT*/

procgplotdata=meansco;

plotMEAN*DAGE/annotate=anno nolegendhaxis=axis1 vaxis=axis2;

run;quit;

/******************************************************************************************************/

Re: GPLOT Plot mean and confidence limits with error bars

pakiprinsesse — Mon, 14 Jan 2019 07:58:26 GMT

Thank you. Changing the variable to character helped me a lot along the way!

Re: GPLOT Plot mean and confidence limits with error bars

PaigeMiller — Mon, 14 Jan 2019 12:58:14 GMT

@pakiprinsesse wrote:

Thank you. Changing the variable to character helped me a lot along the way!

Just in case there's any interest in the PROC GPLOT solution ... I didn't advise to change the variable to character. I advised that you use Y instead of YC, with a numeric variable.