[R-sig-ME] sample weights in lmer

Gasper Cankar gasper.cankar at ric.si
Wed Feb 11 12:31:12 CET 2009


Dear R users,

I found in this mailing list the post from Adam that's been unanswered.
Since I came across same problem, I'll join his kind request for any
help on using sample weights with lmer.

Thanks for any response,

___________________
Gasper Cankar. PhD 
researcher
National Examinations Centre (Ric)
gasper.cankar at ric.si
tel. +386 1 54 84 682
fax. +386 1 54 84 601

Adam's post:

>From aslez at ssc.wisc.edu  Sat Oct 18 00:09:31 2008
From: aslez at ssc.wisc.edu (aslez at ssc.wisc.edu)
Date: Fri, 17 Oct 2008 17:09:31 -0500 (CDT)
Subject: [R-sig-ME] sample weights in lmer
Message-ID:
<4739.128.104.27.172.1224281371.squirrel at webmail.ssc.wisc.edu>

Fellow R users,

I was wondering whether it is possible incorporate sample weights into
the
lmer command.  My understanding is that the existing weights option
follows the types of weights used with the lm command; that is, they are
the type of weights you would use when doing something like weighted
least
squares models.  To the best of my knowledge, sample or design weights
are
not equivalent to the types of model weights we use for WLS.  I have at
least quasi-confirmed this by comparing lmer results to results produced
using HLM.  Without weights, the results are identical. With weights,
they
differ.  I know for a fact HLM is using sample weights.

This maybe a relatively simple question, but I just made the switch to
R,
so the answer isn't completely obvious to me.  My thought is either that
there is a way of linking svydesign commands with the lmer routine, or
that theres a way of "tricking" lmer so that it produces sample-weighted
results using the WLS-type weights command.

Any help would be much appreciated.

Adam




More information about the R-sig-mixed-models mailing list