For delta_p calculation, try the following code:
proc logistic data=&data ;
model &depvar (event="1") = ... ;
ods output parameterestimates=_parms (where=(variable ne "Intercept"));
run;
proc sql noprint; select sum(&depvar)/count(*) into :prior from &data; quit;
data delta_p;
set _parms; retain P0 &prior;
L0=log(P0/(1-P0));
L1=L0+estimate;
P1=exp(L1)/(1+exp(L1));
Delta_P=P1-P0; drop L0 L1;
run;
Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.
Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.
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.