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, is already up and running.

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

st: GLLAMM selection model with endogenous regressor

From   "Marcello Morciano (MED)" <>
To   "" <>
Subject   st: GLLAMM selection model with endogenous regressor
Date   Fri, 16 Nov 2012 05:46:05 +0000

Dear Statalisters,
I would like use GLLAMM for estimating a Heckman model (heckprob) with an endogenous covariate in the selection equation. I am still trying to understand the "language" of Gllamm to which I'm not used to.
Thanks to Sophia and Alfonso's paper ( ) I am able to reproduce results of a heckprob model with GLLAMM. Thanks to Sophia and Xiaohui's paper ( ) I am able to implement a MIMIC model with GLAMM.
I wondering whether I can combine these 2 pieces in a Full information maximum likelihood, instead of applying the Murphy-Topel estimator required in adopting a limited information maximum likelihood, two-step procedure for estimating the following model (in GLLAMM language):
*** MIMIC model
eq load: d1-d7
eq f1:  female age dhs2 dhs3 dhs4 dhs5
#delimit ;
gllamm score d2-d7, i(idauniq) eqs(load) geqs(f1)  link(logit) family(binom) pweight(wt)  adapt nip(30) iterate(100) difficult #delimit cr
where d1-d7 are binary indicators of the latent construct EndVAR_d, and f1 is the vector of covariates for which the conditional mean of EndVAR_d is assumed to vary (the first factor is constrained to unity for identification).
eq fac: end cv
constraint def 1 [id1_1]end=1
eq load: d1-d7
#delimit ;
gllamm resp  Ewalk_c Esit_c age_c female_c cons_c EndVAR_d female_d dhs2_d dhs3_d dhs4_d dhs5_d age_d cons_d, i(id)  weightf(wt)
constr(1) from(startv) long family(binom binom) nr(2,1) link(probit probit) fv(vartype) lv(vartype) nocons
eq(fac) adapt nip(16)      copy;
#delimit cr
Where *_c variables are covariates in the main model; *_d are covariates in the selection model (EndVAR_d included).
Is there a way to do so? How data should be organized, given the fact that to fit a model in GLLAMM, all responses must be stacked in one variable? How the model should be specified in GLLAMM?
Many thanks in advance for your attention.

Marcello Morciano
Research Fellow, Health Economics Group, Faculty of Medicine and Health Sciences, University of East Anglia, Norwich Research Park, Norwich NR4 7TJ

*   For searches and help try:

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