Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down on April 23, and its replacement, statalist.org is already up and running.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

RE: st: endogenous regression with OLS first stage and logit main regression


From   Cameron McIntosh <cnm100@hotmail.com>
To   STATA LIST <statalist@hsphsun2.harvard.edu>
Subject   RE: st: endogenous regression with OLS first stage and logit main regression
Date   Fri, 21 Oct 2011 16:29:24 -0400

Karen,
The error variance in logistic regression is fixed for identification purposes at pi-squared/3, so I believe that may be the reason for your perceived roadblock here:
Fielding, A. (2004). Scaling for residual variance components of ordered category responses in generalised linear mixed multilevel models. Quality and Quantity, 38, 425-433.
Bauer, D.J. (2009). A note on comparing the estimates of models for cluster-correlated or longitudinal data with binary or ordinal outcomes. Psychometrika, 74(1), 97-105.
Allison, P.D. (1999). Comparing logit and probit coefficients across groups. Sociological Methods & Research, 28(2), 186-208.
Mood, C. (2010). Logistic regression: Why we cannot do what we think we can do, and what we can do about it. European Sociological Review, 26(1), 67-82.
Williams, R. (2009). Using heterogeneous choice models to compare logit and probit coefficients across groups. Sociological Methods & Research, 37(4), 531-559.
Long, J.S. (2009). Group Comparisons in Logit and Probit Using Predicted Probabilities. Working paper draft 2009-06-25http://www.indiana.edu/~jslsoc/files_research/groupdif/groupwithprobabilities/groups-with-prob-2009-06-25.pdfhttp://www.indiana.edu/~jslsoc/research_groupdif.htm
Offhand, I'm not 100% sure if you can just substitute in that quantity in the code below and all will be well, but it may work. You may also want to consider a mediational model of the type Z--->X--->Y in -cmp- or gllamm, for example, where you can have both linear and logit link functions in a simultaneous equation framework:
Roodman, D. (2011). Fitting fully observed recursive mixed-process models with cmp. The Stata Journal, 11(2), 159-206.http://www.stata-journal.com/article.html?article=st0224http://ideas.repec.org/c/boc/bocode/s456882.html
www.gllamm.org
My two cents,
Cam
> Date: Fri, 21 Oct 2011 10:05:50 -0700
> From: ruckman@sfu.ca
> To: statalist@hsphsun2.harvard.edu
> Subject: st: endogenous regression with OLS first stage and logit main regression
> 
> Hello,
> 
> I am trying to run an instrumental variables regression with an OLS first stage and logit (or probit if necessary) in the main regression.  We have a system that is "triangular" which means that our endogenous regressor isn't a function of the dependent variable.  Therefore, our system looks like this:
> 
> main regression (logit or probit):
> y1 = y2 + x1 + e
> 
> first stage (OLS):
> y2 = instrument + error
> 
> The problem with using IVPROBIT is that the first stage will use all the exogenous regressors from the main regression as instruments, which we don't want or need as it won't bias the results to use only our one instrument.  There is a very useful discussion about how to use our one instrument in the first stage instead of all the exogenous regressors at the bottom of the FAQ below:
>    http://www.stata.com/support/faqs/stat/ivreg.html
> It says to simply regress the first stage and substitute y2hat for y2 into the main regression.  No problems there.  But we run into trouble when trying to correct the variance-covariance matrix with the "correct mean square error".  The code is written for a system where the first and main stages are both OLS.   Is there a mean square error for logit?  STATA doesn't produce it.  
> 
> Here is the code which recalculates the var-covar matrix after running an OLS main stage regression with y2hat substituted in.  
>  . scalar realmse = r(mean)*r(N)/e(df_r) 
>  . matrix bmatrix = e(b)
>  . matrix Vmatrix = e(V)
>  . matrix Vmatrix = e(V) * realmse / e(rmse)^2
>  . ereturn post bmatrix Vmatrix, noclear
>  . ereturn display
> 
> Any help with the issue with correcting the var-covar matrix would be most appreciated.  
> 
> Thanks!!
> 
> Karen Ruckman
> Associate Professor
> Beedie School of Business, SFU
> *
> *   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/
 		 	   		  
*
*   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/


© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index