If you have a single treatment variable (say, AUC with three levels 1, 2, or 3), you could do it in one PROC LOGISTIC step by fitting a generalized logit model. The following code does that and saves the individual level predicted probabilities for the three AUC levels in data set PS in the variables IP_1, IP_2, and IP_3 which you could then use in LIFEREG.
proc logistic data=test;
model AUC = var1 var2 var3 var4 / link=glogit;
output out=ps predprobs=i;
run;
... View more
