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