Dear Statalist,
I have a question about the scores post-estimation command after running an
ordinal probit. I am trying to calculate the value of the first derivative
for each observation of the ordered probit likelihood function with respect
to the intercept at the maximum likelihood estimates. It seems like the
"predict scores" command does this, however, when I try to double check the
derivatives manually, I get different results. My ordinal probit has 5
categories, and thus I get 4 scores. I am trying to double check the scores
by using the following formula:
(Pdf(cut j - linear predicted scores) - pdf(cut j+1 - linear predicted
scores))/ (cdf(cut j+1 - linear predicted scores) - cdf(cut j - linear
predicted scores)).
I am not sure why I can't get the same results as Stata.
Any information you could provide would be much appreciated.
Thanks,
Colin
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
