I will discuss my experience with Stata's ml method d2 when coding and estimator for a multi-state Markov transition model with unobserved heterogeneity. When analytical derivatives are available, programming a "d2" estimator is in principle straightforward and offers potentially huge rewards in terms of convergence and speed of convergence: When the likelihood is flat, method "d0" may fail to converge (after a many iterations) as numerical derivatives cannot be computed, whereas convergence is often achieved quickly with method "d2". However, when the likelihood function is non-standard, programming a "d2" estimator may be complicated by Stata's limited range of matrix commands. In these cases, the researcher has to be inventive and may have to take a significant "diversion" to compute blocks of the Hessian that should have been straightforward with enhanced matrix capabilities. These "diversions" may be difficult to code and increase evaluation time significantly. With large datasets, this may also push the memory requirements beyond the available limit.
Download Info
To download:
If you experience problems downloading a file, check if you have the
proper application to
view it first. Information about this may be contained
in the File-Format links below. In case of further problems read
the IDEAS help
page. Note that these files are not on the IDEAS
site. Please be patient as the files may be large.
Did you know? Each page is provided with a technical contact, in case something is not right with the supplied information. See under "publisher info".