*Example 1: Lung data reprise from R;

	*Load Data;
DATA lung;
	input time status;
	cards;
306 2
455 2
1010 1
210 2
883 2
1022 1
310 2
361 2
218 2
166 2
170 2
654 2
728 2
71 2
567 2
144 2
613 2
707 2
61 2
88 2
301 2
81 2
624 2
371 2
394 2
520 2
574 2
118 2
390 2
12 2
473 2
26 2
533 2
107 2
53 2
122 2
814 2
965 1
93 2
731 2
460 2
153 2
433 2
145 2
583 2
95 2
303 2
519 2
643 2
765 2
735 2
189 2
53 2
246 2
689 2
65 2
5 2
132 2
687 2
345 2
444 2
223 2
175 2
60 2
163 2
65 2
208 2
821 1
428 2
230 2
840 1
305 2
11 2
132 2
226 2
426 2
705 2
363 2
11 2
176 2
791 2
95 2
196 1
167 2
806 1
284 2
641 2
147 2
740 1
163 2
655 2
239 2
88 2
245 2
588 1
30 2
179 2
310 2
477 2
166 2
559 1
450 2
364 2
107 2
177 2
156 2
529 1
11 2
429 2
351 2
15 2
181 2
283 2
201 2
524 2
13 2
212 2
524 2
288 2
363 2
442 2
199 2
550 2
54 2
558 2
207 2
92 2
60 2
551 1
543 1
293 2
202 2
353 2
511 1
267 2
511 1
371 2
387 2
457 2
337 2
201 2
404 1
222 2
62 2
458 1
356 1
353 2
163 2
31 2
340 2
229 2
444 1
315 1
182 2
156 2
329 2
364 1
291 2
179 2
376 1
384 1
268 2
292 1
142 2
413 1
266 1
194 2
320 2
181 2
285 2
301 1
348 2
197 2
382 1
303 1
296 1
180 2
186 2
145 2
269 1
300 1
284 1
350 2
272 1
292 1
332 1
285 2
259 1
110 2
286 2
270 2
81 2
131 2
225 1
269 2
225 1
243 1
279 1
276 1
135 2
79 2
59 2
240 1
202 1
235 1
105 2
224 1
239 2
237 1
173 1
252 1
221 1
185 1
92 1
13 2
222 1
192 1
183 2
211 1
175 1
197 1
203 1
116 2
188 1
191 1
105 1
174 1
177 1
;
PROC PRINT data=lung(obs=10);

	*Survival Plot and median survival time;
PROC LIFETEST data=lung plots=survival(atrisk cb );
	time time*status(1);

	*Survival time based on years;
PROC LIFETEST data=lung timelist=182.625 365.25 730.50;
	time time*status(1);

	*Hazard Function;
PROC LIFETEST data=lung plots=hazard;
	time time*status(1);
	
	*Culmulative Hazard Function;
PROC LIFETEST data=lung nelson method=pl;
	time time*status(1);
	ods output ProductLimitEstimates=lung_ple;

PROC PRINT data=lung_ple(obs=25);

PROC SGPLOT data=lung_ple;
	series x=time y=CumHaz;	
*---------------------------------------------------------------------------------------------;

*Example 2: Multiple Groups;

	*Get data;
DATA BMT;
	set sashelp.BMT;

PROC PRINT data=BMT(obs=15);

PROC FREQ data=BMT;
	tables Group*Status;

	*Overall Survival Plot;
PROC LIFETEST data=BMT plots=survival(atrisk cb);
	time T*Status(0);
	
	*Survival Function (low and high risk) and rank test;
PROC LIFETEST data=BMT plots=survival(atrisk cb);
	where Group not in ('ALL');
	time T*Status(0);
	strata Group / test=logrank;

	*Hazard Function;
PROC LIFETEST data=BMT plots=hazard;
	where Group not in ('ALL');
	time T*Status(0);
	strata Group;
	
	*Culmulative Hazard Function;
PROC LIFETEST data=BMT nelson method=pl;
	where Group not in ('ALL');
	time T*Status(0);
	strata Group;
	ods output ProductLimitEstimates=BMT_ple;
PROC PRINT data=BMT_ple;
PROC SGPLOT data=BMT_ple;
	series x=T y=CumHaz/ group=Group;

	*Adding a third group;
PROC LIFETEST data=BMT plots=survival(atrisk cb);
      time T*Status(0);
      strata Group / test=logrank adjust=sidak;

	*Plot curves separately;
PROC LIFETEST data=BMT plots=survival(cl cb=hw strata=panel);
      time T*Status(0);
      strata Group;