Hi Everyone,
I have the following simple section of data from the clinic visits:
I computed the incidence rate (per 1000) of the Y, outcome, after counting it as follows and using to the total population of the study area obtained from the census (i.e., 2000) due to a convincing reason provided from the team:
count if Y==1
gen incidence_rate=.
replace incidence_rate= (5/2000)*1000
My first inquiry is I kindly wanted your feedback on whether what I did was right or not.
The second question is in relation to fitting the poisson regression for the above data to see the incidence rate ratio by the available covariates (as logistic regression may not be appropriate- I am open if you have any other modelling suggestion for the data I have). My questions are:
1) How could I create a count outcome variable (let us say this variable Z) using the variables I have in the table above?
2) Since the denominator I used to calculate the incidence risk ratio was obtained from the local area census (i.e. 2000), not the total number of study participants, is there a need to consider this denominator while fitting the poisson regression? If that is the case, I would appreciate your help in how to consider this denominator in the poisson regression Stata code.
Thanks for your help!
I have the following simple section of data from the clinic visits:
| Person_id | year | Outcome (Y) (No/Yes) |
| 001 | 2010 | 0 |
| 001 | 2010 | 1 |
| 001 | 2011 | 1 |
| 002 | 2010 | 0 |
| 002 | 2012 | 1 |
| 003 | 2013 | 1 |
| 004 | 2010 | 0 |
| 004 | 2011 | 0 |
| 004 | 2013 | 0 |
| 005 | 2012 | 1 |
| 005 | 2013 | 0 |
| 005 | 2015 | 0 |
count if Y==1
gen incidence_rate=.
replace incidence_rate= (5/2000)*1000
My first inquiry is I kindly wanted your feedback on whether what I did was right or not.
The second question is in relation to fitting the poisson regression for the above data to see the incidence rate ratio by the available covariates (as logistic regression may not be appropriate- I am open if you have any other modelling suggestion for the data I have). My questions are:
1) How could I create a count outcome variable (let us say this variable Z) using the variables I have in the table above?
2) Since the denominator I used to calculate the incidence risk ratio was obtained from the local area census (i.e. 2000), not the total number of study participants, is there a need to consider this denominator while fitting the poisson regression? If that is the case, I would appreciate your help in how to consider this denominator in the poisson regression Stata code.
Thanks for your help!
