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]]
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: 301-496-8180 E-mail: mh192j at nih.gov -----Original Message----- From: Stef?n Hrafn J?nsson [mailto:shj107 at gmail.com] Sent: Monday, April 30, 2012 6:05 AM To: r-help at r-project.org 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]] ______________________________________________ R-help at r-project.org 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.
On Mon, Apr 30, 2012 at 10:05 PM, Stef?n Hrafn J?nsson <shj107 at gmail.com> 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