/* had to make liver transplant a 1 to get what Tibshirani (1997) got */
options nodate nonumber ps=2000 ls=100;
data pbc;
infile 'K:\pbc.dat.txt';
    input id futime status drug age sex ascites hepato spiders edema
            bili chol albumin copper alk_phos sgot trig platelet
            protime stage;
    age = age/365.25;
if _n_ > 312 then delete;
if status=1 then status=0;
%include 'K:\fsr_phreg91.txt';
%fsr_phreg91(dataset=pbc,model=
drug|age|sex|ascites|hepato|spiders|edema|bili|chol|albumin|copper|alk_phos|sgot|trig|platelet|protime|stage @q,
            gamma=.05,y=futime,censor=status,censorvalue=0,method=1,terms=20);
run;