Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: st: providing raw weights for multivariate meta-analysis
From
Austin Nichols <[email protected]>
To
[email protected]
Subject
Re: st: providing raw weights for multivariate meta-analysis
Date
Sun, 10 Jul 2011 14:02:21 -0400
Nick Darson <[email protected]>:
Have you checked out -gllamm- yet (on SSC and at http://www.gllamm.org/)?
Official commands in Stata 11 does not allow different weights at
different levels, but -xtmixed- in Stata 12 will:
http://stata.com/stata12/multilevel-models-with-survey-data/
My own inclination is not to impose normality using RE/xtmixed but to
run a fixed effects model instead, as it requires weaker assumptions.
You can absorb one dimension of fixed effects with -areg- if you don't
want to include a large number of dummies.
In that case, you can multiply your precision weights times any other
weights, rescale, and specify the product as pweights.
The point estimates should be identical whether specifying pweights or
aweights; only SEs differ (pw=aw+_robust).
See also
http://www.stata.com/statalist/archive/2005-10/msg00079.html
http://www.stata.com/statalist/archive/2007-11/msg00812.html
http://www.stata.com/statalist/archive/2009-05/msg01011.html
and look up Sean Reardon's -hlm- package on SSC.
On Sun, Jul 10, 2011 at 1:43 PM, Nick Darson <[email protected]> wrote:
> Thanks for your help, Stas.
>
> The book I referred to is: Joop J. Hox 2010: Multilevel analysis:
> techniques and applications, 2nd edition, Routledge
>
> To my understanding, "aweights" normalizes the weights (such that the
> sum of the weights =N, the number of observations).
>
> However, it is important for the multivariate meta-analysis approach
> that the weights are not normalized (so called precision weights or
> raw weights).
>
> This is one of the two mentioned requirements by Hox (2010). The
> second requirement for multivariate meta-analyses is the constraint of
> the lowest-level variance to 1 - what is given in Stata - please
> correct me if I am wrong on this).
> Therefore, I was wondering whether the weight issue is the reason that
> Stata cannot be used for multivariate meta-analyses.
>
> Cheers,
> Nick
>
>
>
> On Mon, Jul 11, 2011 at 3:12 AM, Stas Kolenikov <[email protected]> wrote:
>> I personally think that -help weights- explains it all. Variance
>> weights are aweights: they state that the measurement error in a given
>> observation is such and such. The pweights are probability weights,
>> and when you invoke these, the calculation of not only the point
>> estimates, but also the standard errors is different (according to
>> probability sampling theory rather than the likelihood theory). There
>> are also frequency weights, which is essentially the result of
>> -collapse (count) ... , by(*)-. The iweights are "all other weights,
>> whichever these may mean for you" -- there is no strict definition,
>> and the program that allows them usually has its own idea what to do
>> with them.
>>
>> Knowing very little about meta-analysis beyond the fancy name, I have
>> a feeling that you should be looking towards -gllamm- rather than
>> -xtmixed- to find the flexibility with the weights that you need. You
>> might still want to run (the much faster) -xtmixed- first to get good
>> starting values that you'd feed to -gllamm-. Joop Hox (a proper
>> reference is in place, according to the rules of Statalist; I have a
>> vague idea that this is a famous multilevel author) is a somewhat
>> opinionated guy, to my impression (although I am much more opinionated
>> than he is :)).
>>
>> On Sun, Jul 10, 2011 at 7:02 AM, Nick Darson <[email protected]> wrote:
>>> Dear Statalisters,
>>>
>>> I attempt to carry out a multivariate meta-analysis using "xtmixed". I
>>> have three outcome measures (the first level).
>>>
>>> When using a multilevel approach for the multivariate meta-analysis, I
>>> read that it is important to provide raw-weights (non-normalized) of
>>> the inverse sampling variances for the first level (containing the
>>> outcome measures).
>>>
>>> I know that HLM6 has a special function for this ("V-known") for two
>>> and three level models.
>>>
>>> Now, I wanted to know whether "pweights" is an equivalent solution for
>>> xtmixed in Stata?
>>>
>>> I must admit that I am a bit confused by the various descriptions of
>>> the weight functions I found online (aweights, iweights, etc - though
>>> this are not allowed for xtmixed) .
>>>
>>> Moreover, surprisingly, I read in Hox (2010,p. 230) that public
>>> domain software for multilevel analysis does not support the required
>>> options for a multivariate model so far (i.e. providing raw weights
>>> and being able to constrain the lowest-level variance to 1 - however,
>>> both should be feasible).
*
* 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/