Quero calcular as taxas de incidência para apresentar as taxas de risco, a fim de apresentar medidas de risco relativas e absolutas. Vi em outros estudos que essas taxas de incidência podem ser calculadas usando modelos de Poisson com tempo de acompanhamento no modelo como um deslocamento. Então, eu tentei isso em R da seguinte maneira:
library(survival)
# Get example data
data(colon)
colon$status <- ifelse(colon$etype==1,0,1) # set to 0/1 (needed for poisson later on)
# Fit cox model for rx (age + sex adjusted)
coxph(Surv(time,status)~rx+sex+age, data=colon)
# HR (rxLev): 0.92
# HR (rxLev+5FU): 0.74
# Get incidence rates using poisson models with same terms and log(time) as offset
mod <- glm(status~offset(log(time))+rx+sex+age, data=colon, family=poisson)
# Get rates using predict-function
Obs <- predict(mod, data.frame(time=1, rx="Obs", age=mean(colon$age),
sex=mean(colon$sex)), type="response")
Lev <- predict(mod, data.frame(time=1, rx="Lev", age=mean(colon$age),
sex=mean(colon$sex)), type="response")
Lev5FU <- predict(mod, data.frame(time=1, rx="Lev+5FU", age=mean(colon$age),
sex=mean(colon$sex)), type="response")
# Calculate incidence rate ratio's:
Lev/Obs # 0.98
Lev5FU/Obs # 0.84
Eu esperaria que as taxas de incidência sejam semelhantes às taxas de risco do modelo Cox PH com os mesmos termos, mas de alguma forma elas diferem. Estou usando a abordagem correta para calcular taxas de incidência?
Qualquer ajuda seria muito apreciada!