Fitting NHPP-Lack of fitting statistics

Yazdan — Fri, 22 Jul 2016 16:05:39 GMT

Hello all,

I wrote the following code to fit a NHPP on my data. The problem is while it gives me the parameters and their relative errors, it does not provide a fitting error, like Anderson Darling. I want to see how far my fitted model is from the general data (e.g. mean square error). I want to know which model, power-law or log-linear, provides a better fit. I need your help, ASAP. Thank you!

code:
proc reliability data=app1;
         title 'Apple-Dataset';
         unitid id;
       distribution nhpp(pow);
       model time*day(-1);
        mcfplot time*day(-1) / fit=model noconf;
       run;

Regards,

Yazdan

Re: Fitting NHPP-Lack of fitting statistics

SteveDenham — Fri, 22 Jul 2016 19:34:01 GMT

This is kind of a kludgy way around it, and I am not sure that it will work. If you use the parameter estimates to calculate predicted values at each time point, you could get the residuals, and from there calculate PRESS values as a way of comparing. Another would be to look at Kolmogorov-Smirnov tests comparing the observed and predicted values.

Steve Denham

topic Re: Fitting NHPP-Lack of fitting statistics in Statistical Procedures

Fitting NHPP-Lack of fitting statistics

Re: Fitting NHPP-Lack of fitting statistics