Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: st: Trivariate ordered probit model in Stata?
From
Stas Kolenikov <[email protected]>
To
[email protected]
Subject
Re: st: Trivariate ordered probit model in Stata?
Date
Tue, 4 Jan 2011 09:17:47 -0600
On Tue, Jan 4, 2011 at 8:03 AM, Karabulut, Yigitcan
<[email protected]> wrote:
> I have a question regarding the feasibility of a trivariate ordered probit model in Stata 11: I would like to estimate a trivariate probit model where Y1,Y2 and Y3 are ordered variables all of which take the values 1 & 2 & 3:
>
> Y1=X*beta1+epsilon1
> Y2=X*beta2+gamma*Y1hat+epsilon2
> Y3=X*beta3+theta1*Y1hat +epsilon3
>
> However, while constructing the variance covariance matrix, I would like to allow correlation just for -Y1 and Y2- and -Y1 and Y3-. There is no correlation between Y2 and Y3. Therefore, I assume the relation between Y2 and Y3 is seemingly unrelated.
To clarify: is it really the predicted value of Y1 that goes into the
regression for Y2? Or do you mean the continuous analogue of Y1
(typically denoted by a star in econometric textbooks)?
And another clarification: you already have Y1 or Y1hat in the
regressions for Y2 and Y3; do you mean there's some unexplained
correlation left in epsilons?
Have you established identification of this system if all variables
were continuous, with the restrictions you impose? If it is not
identified, then there is no way the ordered version of it will be
identified.
> Do you think whether it makes sense to estimate two seperate bivariate ordered probit models and combine them using the command suest, i.e.:
>
> bioprobit Y1 Y2 X, vce(cluster id)
> estimates store p1
> bioprobit Y1 Y3 X, vce(cluster id)
> estimates store p2
> suest p1 p2
It might work. Better yet, I would try two or three different
approaches and see if they give the same results. -cmp- is the first
thing that jumps to my head when I see stuff like this; -gllamm- is
the second thing.
--
Stas Kolenikov, also found at http://stas.kolenikov.name
Small print: I use this email account for mailing lists only.
*
* 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/