Pratique de Sas Windows ...9.3 Volume 2 - Ined

3.5.4 Amélioration du graphique de type statgraph par modification du template graphique
Devant le fait que les données de survie restent à un niveau élevé de probabilité de survie, il est nécessaire
d’aménager cette échelle de probabilité pour que visuellement, on puisse voir si les deux groupes/strates se
distinguent. Avec l’ ODS Graphic Editor (SGE vu en 3.3), il n’est pas possible de réajuster les valeurs des axes.
Il nous faut donc nous plonger dans la découverte du code du template de ce statgraph de nom
productlimitsurvival et il faudra essayer de comprendre quelque peu sa structure pour y découvrir l’élément à
modifier qui est le code tickvaluelist=(0 .2 .4 .6 .8 1.0))); apparaissant dans l’essai initial.
Voir remarque concernant le GTL, page 117 en fin du point13 : cela peut vous aider face à ce type de code Sas.
proc template;
define statgraph Stat.Lifetest.Graphics.ProductLimitSurvival;
dynamic NStrata xName plotAtRisk plotCensored plotCL plotHW plotEP labelCL labelHW labelEP
maxTime xtickVals xtickValFitPol method StratumID classAtRisk
plotBand plotTest GroupName yMin Transparency SecondTitle TestName pValue;
BeginGraph;
if (NSTRATA=1)
if (EXISTS(STRATUMID))
entrytitle "Kaplan Meier plot" " pour " STRATUMID;
else
entrytitle "Kaplan Meier plot";
endif;
if (PLOTATRISK)
entrytitle "avec le nombre de sujets à risque" / textattrs=GRAPHVALUETEXT;
endif;
layout overlay / xaxisopts=(shortlabel=XNAME offsetmin=.05
linearopts=(viewmax=MAXTIME tickvaluelist=XTICKVALS
tickvaluefitpolicy=XTICKVALFITPOL))
yaxisopts=(label="Probabilité de survie"
shortlabel="Survie" linearopts=(viewmin=0.8 viewmax=1
tickvaluelist=(.8 .87 .88 .89 .90 .91 .92 .93 .94 .95 .96 .97 .98 .99 1.00)));
if (PLOTHW=1 AND PLOTEP=0)
bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME /
modelname="Survie" fillattrs=GRAPHCONFIDENCE
name="HW" legendlabel=LABELHW;
endif;
if (PLOTHW=0 AND PLOTEP=1)
bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME /
modelname="Survie" fillattrs=GRAPHCONFIDENCE
name="EP" legendlabel=LABELEP;
endif;
if (PLOTHW=1 AND PLOTEP=1)
bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / modelname="Survie"
fillattrs=GRAPHDATA1 datatransparency=.55 name="HW" legendlabel=LABELHW;
bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / modelname="Survie"
fillattrs=GRAPHDATA2 datatransparency=.55 name="EP" legendlabel=LABELEP;
endif;
if (PLOTCL=1)
if (PLOTHW=1 OR PLOTEP=1)
bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / modelname="Survie"
display=(outline) outlineattrs=GRAPHPREDICTIONLIMITS name="CL" legendlabel=LABELCL;
else
bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / modelname="Survie"
fillattrs=GRAPHCONFIDENCE name="CL" legendlabel=LABELCL;
endif;
endif;
stepplot y=SURVIVAL x=TIME / name="Survival" rolename=(_tip1=ATRISK _tip2=EVENT)
tip=(y x Time _tip1 _tip2) legendlabel="Survie";
if (PLOTCENSORED=1)
scatterplot y=CENSORED x=TIME / markerattrs=(symbol=plus)
name="Censored" legendlabel="Censuré";
endif;
if (PLOTCL=1 OR PLOTHW=1 OR PLOTEP=1)
discretelegend "Censored" "CL" "HW" "EP" / location=outside halign=center;
else
if (PLOTCENSORED=1)
discretelegend "Censored" / location=inside autoalign=(topright bottomleft);
endif;
endif;
if (PLOTATRISK=1)
innermargin / align=bottom;
blockplot x=TATRISK block=ATRISK / repeatedvalues=true display=(values)
valuehalign=start valuefitpolicy=truncate
labelposition=left labelattrs=GRAPHVALUETEXT
valueattrs=GRAPHDATATEXT (size=7pt) includemissingclass=false;
129 INED Service Informatique Wielki