I used the following output options but only missing values are generated for the fields "peares" and "devres". Could any experience programmer tell what is wrong with that? Thanks very much.
proc logistic data = myDATA descending;
model y = xbeta1-xbeta5 /lackfit rsq;
output out = MODEL_OUTDATA p=prob_def reschi=peares resdev=devres xbeta=xbeta dfbetas=_ALL_ difdev = diff_dev difchisq = diff_chisq;
ods select ModelInfo ResponseProfile ConvergenceStatus FitStatistics RSquare
GlobalTests ParameterEstimates OddsRatios Association LackFitPartition LackFitChiSq;
run;
... View more