Stata The Stata listserver
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: RE: Returned results in gllamm


From   "Mitchell, Michael" <mnm@ucla.edu>
To   <statalist@hsphsun2.harvard.edu>
Subject   st: RE: Returned results in gllamm
Date   Thu, 25 Aug 2005 13:04:46 -0700

Greetings
  I was searching for an answer to this very same question and found
Miguel's query, but no answer. I ran a random intercept, random slope
model with gllamm like this...

. eq cons : one
. eq x    : x
. generate one = 1
. gllamm y x , nrf(2) eqs(x cons) i(id) link(logit) family(bin) adapt
///
  nip(10)

and then got this output...

number of level 1 units = 2000
number of level 2 units = 500
 
Condition Number = 12.83709
 
gllamm model
 
log likelihood = -175.56749
 
------------------------------------------------------------------------
--
           y |      Coef.   Std. Err.      z    P>|z|   [95% Conf.
Interval]
-------------+----------------------------------------------------------
--
           x |     1.4471   .7456141     1.94   0.052  -.0142766
2.908477
       _cons |   3.038546   .6268254     4.85   0.000   1.809991
4.267101
------------------------------------------------------------------------
----
 
 
Variances and covariances of random effects
----------------------------------------------------------------------- 
***level 2 (id)
 
    var(1): .94825463 (.56235789)
    cov(2,1): -.84546937 (.95886005) cor(2,1): -1
    var(2): .75382543 (1.5231523)
----------------------------------------------------------------------

  I think I found the variance of the intercept this way... 

mat list e(chol)

e(chol)[2,2]
            c1          c2
r1   .97378367           0
r2  -.86823121   3.586e-07

. di .97378367 ^2 
.94825464

But this does not work for the variance of the slope...

. di 3.586e-07^2
1.286e-13

Any suggestions would be greatly appreciated!

Best regards,

Michael Mitchell
UCLA ATS Statistical Consulting Group
http://www.ats.ucla.edu/stat/ 


> -----Original Message-----
> From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-
> statalist@hsphsun2.harvard.edu] On Behalf Of miguel foguel
> Sent: Monday, July 25, 2005 10:32 AM
> To: statalist@hsphsun2.harvard.edu
> Subject: st: Returned results in gllamm
> 
> Dear all,
> 
> Could anyone tell me how I can retrieve the estimated variances and
> std-errors of the random "intercepts" from an estimation using
-gllamm-?
> They don't seem to be saved in the e(V) matrix.
> 
> Thank you in advance.
> 
> Miguel
> 
> 
> *
> *   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/



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