1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44
| ods graphics off;
PROC LIFETEST DATA=&Tab METHOD=KM OUTSURV=&TabSortie &affichage;
TIME &Delai*&VarCensure(0);
STRATA &Strate/ TEST=LOGRANK;
FORMAT AAAA BBB.;
RUN;
**-----------correction du graphe, cette partie du code permet de prolonger les graphes -- *;
*-specifier juste la table de sortie lifeftest- pas besoin de modifier le code--**;
DATA &TabSortie;
SET &TabSortie;
IF stratum=1 THEN DO;
survival1=survival;
END;
ELSE IF stratum=2 THEN DO;
survival2=survival;
END;
DROP _censor_;
RUN;
PROC SQL NOPRINT;
SELECT MIN(survival1), MIN(survival2) into: survival1, :survival2
FROM &sortie_lifetest;
RUN; QUIT;
DATA &TabSortie;
SET &TabSortie;
IF ((survival1=. and stratum=1) and survival1^=0) THEN survival1= &survival1;
IF ((survival2=. and stratum=2) and survival2^=0) THEN survival2= &survival2;
RUN;
*---- Affichage du graphe - à modifier--***;
GOPTIONS RESET=all device=activex cback=white;
AXIS1 order=(0 to 200 by 12) minor=none label=(&LabelAxeX);
AXIS2 order=(0 to 1 by 0.2) minor=none label=(a=90 &LabelAxeY);
LEGEND1 LABEL=none value=(label_legend1 label_legend1);
ods rtf startpage=no;
PROC GPLOT data=&TabSortie;
SYMBOL1 l=1 v=dot c=couleur1 mode=include h=0.5 width=4 i=j;
SYMBOL2 l=1 v=dot c=coleur2 mode=include h=0.5 width=4 i=j;
PLOT survival1*&Delai survival2*&Delai /overlay haxis=axis1 vaxis=axis2 legend=legend1 noframe;
TITLE1 &Titre1;
RUN; QUIT; |
Partager