topic How to compute incidence rate in Statistical Procedures

How to compute incidence rate

Ximena_O — Tue, 19 Dec 2017 20:29:42 GMT

Hello,

This is my first time writing here. I have mortality data from a cohort study. I've been trying to calculate incidence rates for mortality with no success. I used proc genmod with lsmeans to estimate the IR using Poisson, but the results are very very different from those obtained by my colleague using STATA.This is the code I've been using:

proc genmod data=muj;
class muerto diabetes(ref=first) educacion(ref=first) sexo(ref=last);
model muerto=edad educacion diabetes/dist=poisson
link=log
offset=logpyears;
lsmeans diabetes/ilink exp cl diff OM;
run;

The stata code used for this is the following:

poisson muerto i.diabetes i.cat_edad i.educacion, exposure(pyears) irr
margins diabetes, predict(ir)

Is there some procedure or code that gives me adjusted incidence rates?

Thanks to all in advance

Re: How to compute incidence rate

Miracle — Mon, 18 Dec 2017 14:26:29 GMT

Hi, it maybe helpful to post the Stata code as well?

Re: How to compute incidence rate

Kurt_Bremser — Mon, 18 Dec 2017 14:28:34 GMT

I moved this to the SAS/STAT community, so that the experts for that may be notified of your question.

Re: How to compute incidence rate

Ximena_O — Mon, 18 Dec 2017 14:30:55 GMT

Thanks 😄

Re: How to compute incidence rate

Ximena_O — Mon, 18 Dec 2017 15:03:53 GMT

Thanks for the tip, just included the stata code in the question. 🙂

Re: How to compute incidence rate

StatDave — Mon, 18 Dec 2017 15:43:08 GMT

The SAS code looks correct and follows what is discussed in this note. But your Stata code looks like the model only has one predictor while your SAS code has three, so that will clearly cause a difference. Additionally, note that the predictive margin computed in Stata is not the same as an LS-mean.

Re: How to compute incidence rate

Ximena_O — Mon, 18 Dec 2017 15:45:31 GMT

What would be the equivalent in SAS fro the margin computed in STATA?
(we've been varying the covariates, so it is possible that the codes are a little different, but the results we get are different anyway)

Re: How to compute incidence rate

StatDave — Mon, 18 Dec 2017 16:08:24 GMT

You would need to fit exactly the same model in both SAS and Stata since any difference in the fitted model can cause a difference, large or small, in estimated values. A predictive margin for a level of your predictor is the average predicted value when all observations in your data have the predictor set to that level. So, you could easily compute the point estimate of the margin by creating a copy of your data set with your predictor set to one level and then compute all the predicted values and average them.

Re: How to compute incidence rate

Ximena_O — Tue, 19 Dec 2017 20:30:20 GMT

Just updated the code for STATA...thanks