Statalist The Stata Listserver


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

Re: st: xtmixed how to fit a certain model which I can fit inSAS


From   "Jonathan Bartlett" <Jonathan.Bartlett@lshtm.ac.uk>
To   <statalist@hsphsun2.harvard.edu>
Subject   Re: st: xtmixed how to fit a certain model which I can fit inSAS
Date   Tue, 27 Jun 2006 11:04:58 +0100

Many thanks to Roberto G. Gutierrez, StataCorp for his response to fit
the model as:

   y(i12) =  u(i2) - u(i1) + e(i12) =  u(i2) + v(i12)
   y(i23) = -u(i2) + u(i3) + e(i23) = -u(i2) + v(i23) 

where

   v(i12) = -u(i1) + e(i12)
   v(i23) =  u(i3) + e(i23)

Unfortunately using this approach to fitting the model doesn't give the
same estimates as using proc mixed in SAS. I believe this is because
using the parametrization proposed by Roberto, var(v) is not constrained
to be at least as large as var(u), or equivalently, var(e) is allowed to
be negative.

In my small dataset unfortunately relaxing the constraint that
var(v)>=var(u) results in a model fit with var(v)<var(u), which implies
var(e)<0 and different variance estimates than I get in proc mixed. I
look forward to the new collinear option in xtmixed, or an alternative
workaround using xtmixed.

Many thanks
Jonathan
*
*   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