Suppose I have the following model:
$Q=1(x'\beta+e>0)$,
$D=1(x'\alpha+\gamma Q+u>0)$
and I want to estimate $\gamma$, the error terms $e,u$ are jointly normal.
If $e$ and $u$ are correlated, can I still estimate $\gamma$ consistently by running a probit using the second equation only?
Why or why not?
If I cannot, how could I estimate $\gamma$ in this case?