Dear StataList,
Referring to the paper by Xu in Statistics in Medicine 2003,
"Measuring explained variation in linear mixed models"
(22:3527-3541), has this measure of explained variation been
implemented in Stata as a postestimation command to xtmixed or
anywhere else?
From xtmixed, Stata gives the residual variance (sigma^2), that is
the variance in the response variable conditional on all the fixed
and random effects. How do I calculate, or convince Stata to produce,
the variance in the response conditional only on the clustering
(sigma_0^2)?
Xu's formula for Omega^2 as a measure of the explained variation is
the ratio of the variance of the response conditional on all fixed
and random effects (sigma^2) and the variance of the response
conditional on the clustering only (sigma_0^2) subtracted from 1. It
is the proportion of the variance conditional on the clustering that
is explained by the covariates. That is:
Omega^2 = 1 - (sigma^2 / sigma_0^2)
This seems to be simple partitioning of the variance that I believe
SAS gives you automatically, but I'm not sure how to get this from
Stata.
Can anyone help me here?
Many thanks in advance,
______________________________
Patrick Phillips Statistician MRC Clinical Trials Unit 222 Euston
Road London NW1 2DA
Email: [email protected] Tel: 0207 670 4637 Fax: 0207 670 4829
*********************************************************************************
This communication is confidential and may contain privileged
information intended solely for the named recipient(s). It may not be
used or disclosed except for the purpose for which it has been sent.
If you are not the intended recipient, you must not copy, distribute,
take any action or reliance on it. If you have received this
communication in error, do not open any attachments but please notify
the MRC Euston Campus ICT Service Desk by e-mailing [email protected]
quoting the sender details, and then delete this message along with
any attached files. E-mail messages are not secure and attachments
could contain software viruses which may damage your computer system.
Whilst every reasonable precaution has been taken to minimise this
risk, The MRC cannot accept any liability for any damage sustained as
a result of these factors. You are advised to carry out your own
virus checks before opening any attachments. Unless expressly stated,
opinions in this message are those of the! e-mail author and not of
the Medical Research Council (MRC).
*********************************************************************************
* * 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/
