26 mai 2013 · Built into the R survival package since survival 2 37-1 ▷ (Available in other R status Multiple event analysis does not always lead to gains
Previous PDF | Next PDF |
[PDF] Appraisal of Several Methods to Model Time to Multiple Events per
This is often analyzed with time-to-first-event survival analysis which ignores possible subsequent events, such as several HF-related hospitalizations Accounting for multiple events provides more detailed information on the disease-control process, and allows a more precise understanding of the prognosis of patients
[PDF] Application of Survival Analysis in Multiple Events - PharmaSUG
There is a growing interest in analyzing recurrent or repeated events data in clinical studies It has been proposed to extend survival models based on Cox proportional hazard approach to handle recurrent multiple events data
[PDF] Multiple and Correlated Events
SURVIVAL ANALYSIS Traditional Introduction Multiple events are allowed An extension of the proportional hazards model is to allow for multiple strata in
[PDF] Recurrent-event analysis with Stata: methods and applications
15 nov 2018 · Recurrent-event analysis with Stata: Recurrent Events in Survival Analysis 2 hospital admission using multiple-failure survival models: a
[PDF] Multiple Outcome Survival Models - MBSW
26 mai 2013 · Built into the R survival package since survival 2 37-1 ▷ (Available in other R status Multiple event analysis does not always lead to gains
[PDF] SUGI 28: Survival Analysis Using Cox Proportional Hazards - SAS
model time until single and multiple hospitalizations using PROC PHREG and tools Originally, the event of interest was death and the analysis consisted of
[PDF] 255-2010: Analysis of Survival Data with Recurrent Events Using
The analysis based on the first recurrence times cannot be used to examine if there is a treatment effect over time from Z on multiple events So, recurrent event
[PDF] suspended drywall ceiling details
[PDF] sustainability in cosmetics industry
[PDF] sustainable cities
[PDF] sustainable city concept/pdf
[PDF] sustainable city ideas
[PDF] sustainable development goals 2019
[PDF] sustainable development goals 2030 pdf
[PDF] sustainable development goals pdf 2018
[PDF] sustainable development goals pdf india
[PDF] sustainable development goals report 2019
[PDF] sustainable development mechanism unfccc
[PDF] sustainable personal care products
[PDF] sustainable solutions corporation
[PDF] sustaining antibodies
Multiple Outcome Survival Models
Terry Therneau
May 2013
Liver Transplant in 1999
IRegion
IStatus
I1: acute failure
I2A: chronic, near death
I2B: chronic, severe
I3: other
I Score I0, 5, or 10 points for compatable blood type
I0{10 points for waiting time
w w w w w w w w w w w w889092949698
YearNumber (1000s)
t t t t t t t t t t t t d d d d d d d d d d d d 0 2 4 6 8 10 12 14Transplants
Waiting
Died on list
Impact of Waiting Time
I All Mayo patients listed from Feb 1990 to Aug 1999 I815 subjects: 636 OLT, 66 death, 37 withdraw, 76 censored
I age, sex, MELD, blood typeI Primary question: Did increased waiting time harm survival?For whom?
I Simple (wrong) analysis: Kaplan-Meier curves of death by epoch, censoring at OLT or withdrawalImpact of Waiting Time
I All Mayo patients listed from Feb 1990 to Aug 1999 I815 subjects: 636 OLT, 66 death, 37 withdraw, 76 censored
I age, sex, MELD, blood typeI Primary question: Did increased waiting time harm survival?For whom?
I Simple (wrong) analysis: Kaplan-Meier curves of death by epoch, censoring at OLT or withdrawal012345
0.0 0.2 0.4 0.6 0.8 1.0Years from listing
DeathA KM is the wrong curve
IAnswers the wrong question
I \What do we think that death rates would be, if the other outcomes ceased to exist?"IWho cares?I
Treats transplant aslost to follow-up
I \Uninformative censoring" assumption IKM censoring other causes
IAnswers a question we rarely want
IThe answer is almost certainly wrong
A KM is the wrong curve
IAnswers the wrong question
I \What do we think that death rates would be, if the other outcomes ceased to exist?"IWho cares?I
Treats transplant aslost to follow-up
I \Uninformative censoring" assumption IKM censoring other causes
IAnswers a question we rarely want
IThe answer is almost certainly wrong
A KM is the wrong curve
IAnswers the wrong question
I \What do we think that death rates would be, if the other outcomes ceased to exist?"IWho cares?I
Treats transplant aslost to follow-up
I \Uninformative censoring" assumption IKM censoring other causes
IAnswers a question we rarely want
IThe answer is almost certainly wrong
Prevalence curves
IPj(t) = Pr(in statejat timet)
IPij(t) = Pr(in statejatt, starting in statei)I
ij(t) = hazard function rate of change from stateitoj, given you are in statei hazard functions don't care about the past I ij(t) =Rt0ij(s) = cumulative hazard
IP=exp() in a 2-state modelI
incidence = synonym for hazard (epidemiology) I \Cumulative Incidence Function" = synonym for prevalence.One of the worst labels in our literature
I confusingImisleading
Iextremely well established
Prevalence curves
IPj(t) = Pr(in statejat timet)
IPij(t) = Pr(in statejatt, starting in statei)I
ij(t) = hazard function rate of change from stateitoj, given you are in statei hazard functions don't care about the past I ij(t) =Rt0ij(s) = cumulative hazard
IP=exp() in a 2-state modelI
incidence = synonym for hazard (epidemiology) I \Cumulative Incidence Function" = synonym for prevalence.One of the worst labels in our literature
I confusingImisleading
Iextremely well established
Prevalence curves
IPj(t) = Pr(in statejat timet)
IPij(t) = Pr(in statejatt, starting in statei)I
ij(t) = hazard function rate of change from stateitoj, given you are in statei hazard functions don't care about the past I ij(t) =Rt0ij(s) = cumulative hazard
IP=exp() in a 2-state modelI
incidence = synonym for hazard (epidemiology) I \Cumulative Incidence Function" = synonym for prevalence.One of the worst labels in our literature
I confusingImisleading
Iextremely well established
0 20 4060
80
Years
Prevalence
0123Transplant
DiedWithdrew
0 20 4060
80
Years
Prevalence
0123Transplant
DiedWithdrew
"KM" of deathDrawing the curves
I Built into the R survival package since survival 2.37-1 I (Available in other R packages prior to this) I Leteventbe a factor variable (class variable) with \censoring" as it's rst level. I Use this as the status variable in thesurvfitfunction I Standard print, plot, etc of usual survival curves the default for plots is uphill rather than downhill. > table(cdata$event) censored death ltx withdraw76 66 636 37
> cfit1 <- survfit(Surv(years, event)1, cdata) > plot(cfit1, mark.time=F, lwd=2, col=c(1,2,4), xlab='Years', ylab='Prevalence') > legend(1.5, .5, c("Transplant", "Died", "Withdrew"), col=c(1,2,4), lty=1, bty='n', cex=1.2, lwd=2) > abline(v=1, lty=2)R and SAS
INo, I don't know how to do this in SAS.I
Two sets of wrenches. Use them both.
Anyone who refuses to use anything but SAS is running a race with leg irons. INeither is easy to learn in toto
IWell crafted, sensible, readable languages.
(SAS macro excepted)IThe problem is that they are so very BIG.
IFind a good book and use only what you need.
I Most of what I will show in the course can be done in either. II will usually show R output
Iit is more compact and easier to t onto the slideR and SAS
INo, I don't know how to do this in SAS.I
Two sets of wrenches. Use them both.
Anyone who refuses to use anything but SAS is running a race with leg irons. INeither is easy to learn in toto
IWell crafted, sensible, readable languages.
(SAS macro excepted)IThe problem is that they are so very BIG.
IFind a good book and use only what you need.
I Most of what I will show in the course can be done in either. II will usually show R output
Iit is more compact and easier to t onto the slideTransplants (death in red)0
20 4060
80
Years
Prevalence
01231990-92
1993-95
1996-97
1998-99
Death0
2 4 6 8 Years Death 012390-92
93-95
96-97
98-99
cfit2 <- survfit(Surv(years, event)period, cdata) plot(cfit2[, 1:2], mark.time=F, col=c(2,1,2,3,2,4,2,6), lwd=1:2, xmax=3, xaxt='n', xlab="Years", ylab="Prevalence", yscale=100) axis(1, 0:3, 0:3) text(c(.1, 1.4, 1, 2), c(.9, .88, .7, .6), levels(period), col=c(1,3,4,6), cex=1.2) cfit2[,1:2] records n.max n.start events median 0.95LCL 0.95UCL death, period=1990-92 170 170 170 170 NA NA NA ltx, period=1990-92 170 170 170 170 0.123 0.151 0.101 death, period=1993-95 244 244 244 243 NA NA NA ltx, period=1993-95 244 244 244 243 0.381 0.435 0.337 death, period=1996-97 210 210 210 194 NA NA NA ltx, period=1996-97 210 210 210 194 0.542 0.805 0.460 death, period=1998-99 191 191 191 132 NA NA NA ltx, period=1998-99 191 191 191 132 1.221 1.525 0.997 I