FILENAME goptions URL "http://192.38.117.59/~linearpredictors/datafiles/goptions.sas";
%include goptions;
    
FILENAME url URL "http://192.38.117.59/~linearpredictors/datafiles/readFever.sas";
%INCLUDE url;

data fever;
set fever; 

if death<0 then delete;
/* Definition below makes intercept correspond to mage=30 */
mage30=mage-30;
RUN;

proc genmod data=fever;
    class parity;
    model episodes=parity mage30 / dist=poisson link=log;
    output out=fit p=predicted reschi=residual;
run;

proc rank data=fit out=rankfit groups=10;
    var predicted;
    ranks predicted_group;
run;

proc means nway noprint N min max sum data=rankfit;
    class predicted_group;
    var episodes predicted;
output out=chi2 N=number min=min_yobs min_predicted
                        max=max_yobs max_predicted
                        sum=sum_yobs sum_predicted;
run;

data chi2;
set chi2;

residual=(sum_yobs-sum_predicted)/sqrt(sum_predicted);
run;

proc print data=chi2;
var number min_predicted max_predicted _FREQ_ sum_yobs sum_predicted residual;
run;