next up previous
Next: Multinomial Logistic Regression Up: Binomial Logistic Regression Previous: The Newton-Raphson Method

Caveats

There are two cautionary notes to consider during the iteration procedure. First, it is possible for a parameter estimate to tend to infinity. This is usually a sign that the model is either poorly specified or is behaving badly due to data sparseness in one or more populations. Obviously, a parameter that tends to infinity will never converge. However, it is sometimes useful to allow a model to converge even in the presence of infinite parameters. To accomplish this, each estimate can be tested against a threshold above which it is considered to be infinite. At that point, the iterations can continue while holding the infinite parameter constant, ignoring its new values in susbsequent iterations, and exempting it from the global test for convergence. The SAS System uses two criteria to test whether a parameter estimate is tending to infinity in PROC CATMOD with the /ML option. If (i) the absolute value of the estimate exceeds five divided by the range of the corresponding independent variable, and (ii) the standard error of the estimate is at least three times greater than the estimate itself.

A second cautionary note deals with a limitation of the Newton-Raphson method. Given certain conditions, it is possible for a given estimate to overshoot the true root in such a way that subsequent iterations enter into a repeating cycle that will never converge. To counter this possibility at each iteration, verify that the value for the likelihood function evaluated at that point is in fact higher than it was during the previous iteration. If at any point the likelihood decreases, this is a sign that the iterations have lost track of the true root and are in danger of converging to a local maximum or not converging at all. One strategy for dealing with this is to apply a ``step-halving'' function wherein half the distance between the current and prior estimates is tested. If the likelihood at that point is still lower than for the last iteration then half the distance again is tested. This continues for a reasonable number of ``sub-iterations'' until reaching a point where the likelihood does increase. This is again the method used by the SAS System in PROC CATMOD with the /ML option. We will look at these two caveats in more detail in the section on implementation.


next up previous
Next: Multinomial Logistic Regression Up: Binomial Logistic Regression Previous: The Newton-Raphson Method

Scott Czepiel
http://czep.net/contact.html