Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: multiple weights per person in GEE?


From   "Joseph Coveney" <jcoveney@bigplanet.com>
To   <statalist@hsphsun2.harvard.edu>
Subject   Re: st: multiple weights per person in GEE?
Date   Mon, 20 Jul 2009 09:35:33 +0900

Ariel Linden wrote (excerpted):

I have been using GLM with vce(cluster) with the IPTW weight, but the SE is
much larger than that produced in SAS using GEE. For example, with a beta
coeficient for a treatment variable of 2.47, GLM in stata gives me a SE of
0.484 (CI = 1.53, 3.43) while GEE in SAS gives me SE of 0.013 (CI = 2.45,
2.50).
 
This is a pretty meaningful difference, and in several models this can
change the treatment effect from being positive to one of non significance.

--------------------------------------------------------------------------------

I take it that the difference you're seeing in SEs with identical point
estimates is between

PROC GENMOD;
. . .
REPEATED SUBJECT = . . . / TYPE = IND;
SCWGT . . .;

and

glm . . . [<weight>= . . .], cluster(. . .) . . .

If so, then these are indeed larger differences than would be expected if the
two packages mean the same thing by "weight" in this context.  You've probably
already considered the following and more, but just in case:

1. What kind of weights are you declaring the IPTW to be in Stata?  Fewell et
al. (2004) used Stata's -pweight-.

2. Related to that, does PROC GENMOD need scaling of the weights so that they
sum to the number of observations?

3. Is it possible to cajole Stata into allowing the time-varying weights that
you want by viewing the observation time points in the same manner as waves of a
survey and setting the model up as a survey data analysis task?

Joseph Coveney

Z. Fewell, M. A. Hernán, F. Wolfe, K. Tilling, H. Choi, J. A. C. Sterne.  2004.
Controlling for time-dependent confounding using marginal structural models.
_The Stata Journal_ 4(4):402?420.




*
*   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/



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index