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

Re: st: not getting an expected result with gllamm


From   Joseph Coveney <jcoveney@bigplanet.com>
To   Statalist <statalist@hsphsun2.harvard.edu>
Subject   Re: st: not getting an expected result with gllamm
Date   Mon, 03 Oct 2005 11:18:28 +0900

Josť Maria Pacheco de Souza wrote (excerpted):

I am trying an artificial data set in order to understand gllamm from
scratch, but the resulting output is not in accord with the MlWiN  output,
which seems correct to me:
[deleted]
which gives five paralel straight lines (one for each group), each one with
slope=1 and intercepts 4, 2, 0 -2, 4.
I used the commands
         eq inter:cons
.         eq slope:x
       . gllamm  y x, i(group) eqs(inter slope) nrf(2) ip(m)
and got the following output:

snip...

--------------------------------------------------------------------------------

I take it that MLwiN either failed to converge, or gave residual variance
and covariance as zero and dropped the fixed effects intercept term and gave
a group variance of something like 10.  Could it be that -gllamm- is giving
you these results because of its use of numerical differentiation in its
calls to -ml-?

By the way, what is the option -ip(m)-?  m is not among the options
for -ip()- in either the online help file or the pdf user's manual that I
could find.  And why does -ip(g)- give an error? (version 2.3.10 SRH 4 jan
2005)

Joseph Coveney

*
*   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