    # Re: st: gllamm and marginal effects

 From David Jaeger To statalist@hsphsun2.harvard.edu Subject Re: st: gllamm and marginal effects Date Mon, 6 Mar 2006 14:55:55 +0100

Ok, here's the answer to my own question. One can use gllamm to estimate the random effects probit
(and random effects tobit) models (among many others, of course), and then send the parameter estimates
through xtprobit or xttobit to get the marginal effects that one can calculate after estimation with those commands.

For random effects probit, the following should work:

matrix a=e(b)
local n=colsof(a)
matrix a[1,`n']=ln(a[1,`n'])
xtprobit y x1 x2 x3 , re i(id) from(a,copy) intpoints(30) iterate(0)
mfx compute, predict(pu0)

And for random effects Tobit (where the dependent variable is left- censored only -- trivial to modify for
right and left censoring):

gen var=cond(movedist==0,2,1)

gllamm ycens x1 x2 x3 , i(id) fam(gauss binom) link(ident sprobit) lv (var) fv(var) adapt
matrix aa=e(b)
local n=colsof(aa)
matrix a=aa
matrix a[1,`n'-1]=aa[1,`n']
matrix a[1,`n']=exp(aa[1,`n'-1])
xttobit ycens x1 x2 x3 , re i(id) ll(0) intpoints(30) from(a,copy) iterate(0)
mfx compute, predict(pr0(0, .))
mfx compute, predict(e0(0,.))
mfx compute, predict(ys(0,.))

Hope someone finds this useful!

--dj

On 3 Mar 2006, at 15:00, David Jaeger wrote:

I'm using gllamm to estimate a random effects probit model. Anybody know if there's
an easy way to get marginal effects with gllamm, equivalent to using mfx compute, predict(pu0) after
using xtprobit?

__________________________________________________
Associate Professor of Economics and Public Policy
College of William and Mary
Williamsburg, VA 23187-8795 USA
net: www.djaeger.org
ph: +1 757 221 2375
__________________________________________________

*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
```*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/
```