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; /******************************************************************************************************/ /******************************************************************************************************/
... View more