Dear Communities:
Greetings all!
I have data in the following format.
Patient_ID status level prediction
165201 0 0 0.0090204072
165201 0 1 0.4682205524
165201 0 2 0.8122016298
165201 0 3 0.9751736188
165205 1 0 0.0081253999
165205 1 1 0.4420893893
165205 1 2 0.7955942710
165205 1 3 0.9724900318
165208 2 0 0.0011220970
.
. . .
165210 3 0 0.0980113921
165210 3 1 0.3479984169
165210 3 2 0.8289891254
165210 3 3 0.9509891254
...etc
Here I have two questions:
1. The "prediction" variable is the predicted cumulative probability for ordinal response 0 to 4(but 4 is reference due to the fact that it is not displayed in the table). But I want create a variable that change the cumulative probability into non-cumulative probability, that is, example for the first patient_Id with the 0 level: cumulative probability at level 1 minus at level 0, for the first patient_Id with the second level: cumulative probability at level 2 minus at level 1, for the first patient_Id with the third level: cumulative probability at level 3 minus at level 2, and finally one minus cumulative probability at level 3. The same procedure for all Patient Id.
Finally
2. Based on the created variable in "1" I want to select a level with maximum probability so that I can get one probability value per patient_Id.
Thank you for your help!
I'm happy to justify it if it is not clear.
data want;
set have;
by patient_id;
prev_pred=lag(prediction);
if first.patient_id then do;
non_cumulative_prediction=prediction;
output;
end;
else if not last.patient_id then do;
non_cumulative_prediction=prediction-prev_pred;
output;
end;
else if last.patient_id then do;
non_cumulative_prediction=prediction-prev_pred;
output;
level=4;
non_cumulative_prediction=1-prev_pred;
output;
end;
drop prev_pred;
run;
proc summary nway data=want;
class patient_id;
var non_cumulative_prediction;
output out=maximum max=max_prediction
maxid(non_cumulative_prediction(level))=level_at_max_pred;
run;
data want;
set have;
by patient_id;
prev_pred=lag(prediction);
if first.patient_id then do;
non_cumulative_prediction=prediction;
output;
end;
else if not last.patient_id then do;
non_cumulative_prediction=prediction-prev_pred;
output;
end;
else if last.patient_id then do;
non_cumulative_prediction=prediction-prev_pred;
output;
level=4;
non_cumulative_prediction=1-prev_pred;
output;
end;
drop prev_pred;
run;
proc summary nway data=want;
class patient_id;
var non_cumulative_prediction;
output out=maximum max=max_prediction
maxid(non_cumulative_prediction(level))=level_at_max_pred;
run;
Dear PaigeMiller:
Appreciated!!! Thank you very much in advance!
Registration is open! SAS is returning to Vegas for an AI and analytics experience like no other! Whether you're an executive, manager, end user or SAS partner, SAS Innovate is designed for everyone on your team. Register for just $495 by 12/31/2023.
If you are interested in speaking, there is still time to submit a session idea. More details are posted on the website.
Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.
Find more tutorials on the SAS Users YouTube channel.