[R] analysis of data with observation weights
Thomas Lumley
tlumley at u.washington.edu
Fri Nov 15 19:13:38 CET 2002
On Fri, 15 Nov 2002, bojaniss wrote:
> Hello Thomas,
>
> Friday, November 15, 2002, 4:35:13 PM, you wrote:
>
> TL> The weights argument to lm and glm will give the right point estimates.
> TL> The standard errors will potentially be wrong. This can be fixed with
> TL> `sandwich' standard errors, so one option is to use gee() with each
> TL> observation being in a `group' on its own. Similarly, the `robust'
> TL> standard errors in coxph() will allow probability-weighted survival
> TL> analyses.
>
> TL> The sandwich standard errors used by gee() are not quite the same as the
> TL> ones used by survey samplers, but they are very similar and they are
> TL> consistent estimates of the same thing.
>
> TL> The usual linear model standard errors are often pretty good even for
> TL> probability weighting as long as important covariates aren't strongly
> TL> associated with the weights.
>
>
> TL> -thomas
>
> Where can I find the gee() function, it's not in base package nor in
> any packages I have installed.
>
In the gee package.
-thomas
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list