Ich möchte Inzidenzraten berechnen, die entlang der Gefährdungsquoten dargestellt werden, um sowohl relative als auch absolute Risikomessgrößen darzustellen. Ich habe in anderen Studien gesehen, dass solche Inzidenzraten unter Verwendung von Poisson-Modellen mit Nachbeobachtungszeit im Modell als Offset berechnet werden können. Also habe ich das in R wie folgt versucht:
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
Ich würde erwarten, dass die Inzidenzratenquoten den Hazard Ratios des Cox PH-Modells mit denselben Begriffen ähnlich sind, sich aber irgendwie unterscheiden. Benutze ich den richtigen Ansatz zur Berechnung der Inzidenzraten?
Jede Hilfe wäre sehr dankbar!