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


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

Re: st: Estout after xtmixed with ar(1) residuals


From   Tom Trikalinos <ttrikalin@gmail.com>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: Estout after xtmixed with ar(1) residuals
Date   Thu, 18 Mar 2010 12:23:53 -0400

Just guessing here -- you should be able to work out quickly if I'm
right or wrong.

Typically in an ML routine you would optimize for t=log( of an SD or a variance)
this guarantees that the exp(t) (the SD or the variance) is non negative.

So I'm guessing that to get your variance or SE it should be exp(r_lns2ose)

Similarly for correlations. If they are allowed to range between -1
and 1, one probably optimizes for z=tanh^-1 of the correlation, as the
tanh(z) is bounded between -1 and 1.
However, in an autoregressive model the correlation should be between
0 and 1, no? So I would program the optimization so that the
back-transformation is 0.5*tanh(z)+0.5  -- bounded between 0 and 1.  I
dunno what Stata does.

Please, send an e-mail to the list when you've figured out.



Thomas A Trikalinos MD, PhD

Co-Director Tufts Evidence-based Practice Center
Associate Director, Center for Clinical Evidence Synthesis

Institute for Clinical Research and Health Policy Studies
Tufts Medical Center | 800 Washington St | Boston, 02111 MA

Phone: +1 617 636 0734
Fax:   +1 617 636 8628
email: ttrikalinos@tuftsmedicalcenter.org




On Thu, Mar 18, 2010 at 8:17 AM, Dimitris Pavlopoulos
<Dimitris.Pavlopoulos@soc.kuleuven.be> wrote:
> Dear all,
>
> I am trying to use estout after xtmixed (STATA 11 /SE). I am using a model with random slopes and an ar(1) stucture for the residuals. Stata 11 allows the estimation of separate correlation and residual standard deviation according to the values of a categorical variable. I have used this feature for a binary variable.
>
> My problem is that I cannot find which function is estout using to transform the standard deviation of the residuals that corresponds to value 1 of the grouping variable. So, for the residuals estout present:
>
> lnsig_e: variance for residuals corresponding to value 0 of grouping variable
> r_atr1: autoregressive correlation corresponding to value 0 of grouping variable
> r_lns2ose: variance for residuals corresponding to value 1  of grouping variable
> r_atr2: autoregressive correlation corresponding to value 1 of grouping variable
>
> for lnsig_e, I need to exponantiate to get the standard deviation
> For r_atr1 and r_atr2, I am using the tanh()  function.
>
> Does anybody know what function do I have to use to get the standard deviation from r_lns2ose?
>
> Thanks in advance?
>
> best regards,
> Dimitris
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
>

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   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   |   Site index