'stratified cox (PWP-TT) with interactions
im modeling recurrent events (hospital readmission) and i have some trouble with my pwp-tt model. My objective is to compare HR between strata. I have CP data format, it looks like this:
| id | start_date | stop_date | terap_dis | event | tstart | tstop | status |
|---|---|---|---|---|---|---|---|
| 1 | 16078 | 16377 | 1 | 1 | 0 | 299 | 1 |
| 1 | 18019 | 18213 | 0 | 2 | 299 | 2135 | 0 |
| 2 | 15719 | 16037 | 0 | 1 | 0 | 318 | 1 |
| 2 | 16049 | 16129 | 0 | 2 | 318 | 410 | 1 |
| 2 | 16212 | 16275 | 0 | 3 | 410 | 556 | 0 |
terap_dis is a binary variable that is 1 when a pacient reach a therapeutic discharge, event is the stratification variable, and status is the censoring variable (1 when and individual have a subsequent event.
coxph(Surv(tstart,tstop,status==1) ~ terap_dis+ terap_dis:strata(event)+cluster(id)+strata(event), data = data)
I think if i want to know the HR for terap_dis = 1 vs terap_dis = 0 for the second strata, i should calculate exp(-1.03320+2*0.74589) = 1.71. Is this approach correct? because it dont make sense that a pacient have and increased risk of readmission if they reach therapeutic discharge. One posibility is that previous admissions have greater risk of readmissions. I'm very confused.
Its there a way to plot survival and hazard plot for diferent strata? Pd: This is just a simplified model, because i have more covariates, but the effect is the same.
Thank you very much for your time
Sources
This article follows the attribution requirements of Stack Overflow and is licensed under CC BY-SA 3.0.
Source: Stack Overflow
| Solution | Source |
|---|
