Hi Stata forum,
I want to estimate annual incidence rates of amputation (lea) in a dataset of people with kidney disease.
I have used a log poisson generalized linear model with robust standard errors and included an offset term with log exposure time. The code I have used is below. I am not sure, however, how to now estimate incidence rates (per year) as a post estimation command?
stset dox1, fail(lea1==1) origin(born) entry (entry) scale(365.25) id(id)
stsplit _year, after(time=d(1/1/2000)) at(0(1)15) trim
replace _year=2000 + _year
gen _y=_t-_t0
gen log_y = ln(_y)
glm lea1 _year i.age i.sex i.race, family(poisson) link(log) offset(log_y) vce(robust)
Thanks in advance
Jess
I want to estimate annual incidence rates of amputation (lea) in a dataset of people with kidney disease.
I have used a log poisson generalized linear model with robust standard errors and included an offset term with log exposure time. The code I have used is below. I am not sure, however, how to now estimate incidence rates (per year) as a post estimation command?
stset dox1, fail(lea1==1) origin(born) entry (entry) scale(365.25) id(id)
stsplit _year, after(time=d(1/1/2000)) at(0(1)15) trim
replace _year=2000 + _year
gen _y=_t-_t0
gen log_y = ln(_y)
glm lea1 _year i.age i.sex i.race, family(poisson) link(log) offset(log_y) vce(robust)
Thanks in advance
Jess
Comment