| 12345678910111213141516171819202122 |
- % Local Regression and Likelihood, Figure 7.4.
- %
- % Censored Local Likelihood
- %
- % Author: Catherine Loader
- %
- % NEED: this is given problems for unknown reasons. something
- % to do with censored observations and the ibeta routine.
- % if i replace this by
- % res[ZLIK] = th*y-y*log(p)
- % res[ZDLL] = y*p
- % res[ZDDLL] = y*p*(1-p)
- % works fine. Memory leak somewhere??
- % also works if i ignore censoring.
- load border;
- fit = locfit(day,runs,'cens',no,'family','geom','alpha',0.7);
- figure('Name','fig7_4: censored local likelihood;' );
- lfplot(fit);
- xlabel('Date');
- ylabel('Runs');
|
