data bis679; input time status treatment; cards;
5 1 1
6 1 1
12 1 1
15 0 1
20 1 1
25 1 1
26 1 1
30 0 1
31 0 1
4 1 0
5 1 0
10 0 0
18 1 0
25 1 0
26 1 0
29 0 0
30 1 0
33 1 0
; run;
*proc lifetest gives kaplan-meier stats and log rank test;
proc lifetest plots=(s) graphics;
title ‘Comparing Survival Experience of Two Treatments’;
time time*status(0);
strata treatment;
*proc phreg will do Cox proportional hazards model;
proc phreg;
model time*status(0) = treatment;
run;
*proc sgplot gives some plots;
*https://blogs.sas.com/content/sascom/2011/08/22/how-to-make-a-cluster-grouped-bar-chart-graph-using-sas-sg-procedures/;
proc sgplot;
vbar time/group=treatment /*groupdisplay=cluster*/;
run;