Using GEE with sample weights

[Stefán Hrafn Jónsson]

Dear R community

I am using the gee package to run logistics regression on paired cases from
a panel sample.

We are getting request from a reviewer to use sample weights to account for
different level of attrition.

After searching the documentation I am unable to find a way to incorporate
sample weights into the gee formula. I found a way to incorporate precision
weight but I understand that I cant use that function for sample weight.

I appreciated all help on this issue.

Regards
Stefan Jonsson

        [[alternative HTML version deleted]]

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[Haynes, Maurice (NIH/NICHD) [E]-2]

See the function geeglm() in package:geepack

Maurice

O. Maurice Haynes, Ph.D.
Statistician, Child and Family Research Section
Laboratory for Comparative Ethology
Eunice Kennedy Shriver
National Institute of Child Health and Human Development
NIH, DHHS
6705 Rockledge Drive, Suite 8030
Bethesda, Maryland  20892
Phone: [phone]
E-mail: [hidden email]

-----Original Message-----
From: Stefán Hrafn Jónsson [mailto:[hidden email]]
Sent: Monday, April 30, 2012 6:05 AM
To: [hidden email]
Subject: [R] Using GEE with sample weights

Dear R community

I am using the gee package to run logistics regression on paired cases from a panel sample.

We are getting request from a reviewer to use sample weights to account for different level of attrition.

After searching the documentation I am unable to find a way to incorporate sample weights into the gee formula. I found a way to incorporate precision weight but I understand that I cant use that function for sample weight.

I appreciated all help on this issue.

Regards
Stefan Jonsson

        [[alternative HTML version deleted]]

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[Thomas Lumley-2]

On Mon, Apr 30, 2012 at 10:05 PM, Stefán Hrafn Jónsson <[hidden email]> wrote:

> Dear R community
>
> I am using the gee package to run logistics regression on paired cases from
> a panel sample.
>
> We are getting request from a reviewer to use sample weights to account for
> different level of attrition.
>
> After searching the documentation I am unable to find a way to incorporate
> sample weights into the gee formula. I found a way to incorporate precision
> weight but I understand that I cant use that function for sample weight.

For gee it doesn't matter whether you have sampling weights or
precision weights:
  -  the point estimate is the same (as usual for linear estimating functions)
  -  the standard errors are not model-based and they are also the
same (the standard error estimates are the same as the
Horvitz--Thompson + linearization estimates that you would use for
sampling)

If you wanted to use other features of a complex sampling design you'd
need different software, but not just for per-observation sampling
weights.

    -thomas

--
Thomas Lumley
Professor of Biostatistics
University of Auckland

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.