-
Notifications
You must be signed in to change notification settings - Fork 0
/
project.sas
66 lines (56 loc) · 1.62 KB
/
project.sas
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
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
proc import out=cirrhosis
datafile="C:/Users/ntlmp/Desktop/Datasets/cirrhosis.csv"
dbms=csv replace;
run;
data cirrhosis;
set cirrhosis;
if Status = "D" then censored = 0;
else if Status = "C" or Status = "CL" then censored = 1;
if Stage = 1 then Stage_level = "Stage1";
else if Stage = 2 then Stage_level = "Stage2";
else if Stage = 3 then Stage_level = "Stage3";
else if Stage = 4 then Stage_level = "Stage4";
run;
proc print data=cirrhosis; run;
/* fitting overall survival curve */
proc lifetest data=cirrhosis plots=(survival);
time N_Days*censored(1);
run;
/* KM curves and log-rank test by Drug v. Placebo */
proc lifetest plots=(survival);
time N_Days*censored(1);
strata Drug;
run;
/* KM curves and log-rank test by Sex */
proc lifetest plots=(survival);
time N_Days*censored(1);
strata Sex;
run;
/* KM curves and log-rank test by Hepatomegaly */
proc lifetest plots=(survival);
time N_Days*censored(1);
strata Hepatomegaly;
run;
/* KM curves and log-rank test by Ascites */
proc lifetest plots=(survival);
time N_Days*censored(1);
strata Ascites;
run;
/* fitting Cox model */
/* estimating beta coefficients */
proc phreg data=cirrhosis outest=betas;
class Drug(ref="Placebo")
Sex(ref="F")
Ascites(ref="N")
Hepatomegaly(ref="N")
Spiders(ref="N")
Edema(ref="N")
Stage_level(ref="Stage1");
model N_Days*censored(1) = Drug Age Sex Ascites Hepatomegaly Spiders Edema Bilirubin
Cholesterol Albumin Copper Alk_Phos SGOT Tryglicerides Platelets Prothrombin Stage_level;
baseline out=outdata survival=Sbar;
run;
proc print data=betas;
run;
proc print data=outdata;
run;