Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.


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

st: Dealing with exposure differences using logarithmic offsets versus rescaling


From   "Dimitriy V. Masterov" <[email protected]>
To   Statalist <[email protected]>
Subject   st: Dealing with exposure differences using logarithmic offsets versus rescaling
Date   Wed, 14 Aug 2013 16:07:12 -0700

I am trying to test my intuition that logarithmic offsets should yield
the same estimates as dividing the outcome by exposure and using
exposure to weight.

This seems to work for count variable BI with the regular poisson or
its glm equivalent. However, when I try it with a continuous variable
GMB (measuring expenditures) using GLM with gamma family and log link,
the equivalence goes away. I am at a loss for why that is the case.
Any ideas or references would be much appreciated:

Here's my code:

gen lnt1 = ln(t1);
lab var lnt1 "Treatment Offset (ln of Days)";

gen bipd1 = bi1/t1;
gen gmbpd1 = gmb1/t1;

/* Simple Poisson Model for BI */
poisson bi1 i.treat, offset(lnt1) robust;
poisson bipm1 i.treat [iw=t1], robust;

/* Poisson Using GLM for BI */
glm bi1 i.treat, offset(lnt1) link(log) family(poisson) robust;
glm bipm1 i.treat [iw=t1], link(log) family(poisson) robust;

/* GMB using GLM */
glm gmb1 i.treat, offset(lnt1) link(log) family(gamma) robust;
glm gmbpm1 i.treat [iw=t1], link(log) family(gamma) robust;
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index