st: Generating Murphy-Topel standard errors in Stata

["Breanna Pellegrini" <breanna.pellegrini@rmit.edu.au>]

Dear Statalist,
 
I have a two-step estimation setup where the first stage is modelled using a binary probit specification. I then use the predicted values from the first stage model as an explanatory regressor in the second stage. The second stage is estimated using a linear regression model. Since I have included a generated regressor in the second stage, the standard errors from the second stage model will need to be adjusted using the Murphy-Topel variance estimator.
 
The Stata Journal paper by Hole (2006) titled "Calculating Murphy-Topel variance estimates in Stata: A simplified procedure" (http://www.stata-journal.com/sjpdf.html?articlenum=st0114) provides Stata code for calculating the Murphy-Topel standard errors for a two-step model with a logit model in the first stage and then a Poisson model in the second stage. Additional code changes are also suggested for a number of variations to the model specification at both stages of estimation.
 
I have been able to adapt the code from Hole (2006) for my particular model setup and the code runs fine and Murphy-Topel standard errors are produced. However, specific code amendments for the situation of having a linear regression model in the second stage were not discussed in the Hole (2006) paper. So I am a little unsure as to whether I have properly amended my code for the setup of a probit model in the first stage and a linear regression model in the second stage. In particular:

1) When calculating the score for the second stage linear regression model can I directly use the second stage code used in Hole (2006, pp. 523) and simply replace "poisson" with "regress"?
 
2) Would I need to also adjust the C-matrix and R-matrix calculations to take into account the residual variance parameter in the linear regression model? Hole (2006) discusses how to account for the auxiliary (dispersion) parameter if a negative binomial model is estimated in the second stage, but how would this be done in a linear regression setup?
 
3) In both stages of my estimation I am also using clustered standard errors to allow for intra-group correlation between twenty different regional sites. Are there any further code adjustments I would need to make to calculate the Murphy-Topel estimates in this case?
 
 
Any guidance from those who may have previous experience in dealing with and coding the Murphy-Topel variance estimates would be greatly appreciated.
 
 
Kind regards,
Breanna


 
 
 

Breanna Pellegrini
School of Economics, Finance and Marketing

RMIT University
379-405 Russell Street
Melbourne VIC 3000 Australia

Tel     +61 3 9925 0944
Fax    +61 3 9925 5986



*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/