View source for Models for time-to-event data

<!-- Menu for the Observations chapter -->
<sidebarmenu>
+[[Modeling the observations]]
*[[Modeling the observations| Introduction ]] | [[ Continuous data models ]] | [[Models for count data]]  | [[Model for categorical data]]  | [[Models for time-to-event data]] | [[Joint models]]  
</sidebarmenu>













Here, observations are the "times at which events occur". An event may be one-off (e.g., death, hardware failure) or repeated (e.g., epileptic seizures, metro strike).

<br>
==Single event==


To begin with, we will consider a one-off event.
Depending on the application, the length of time to this event  may be called the ''survival'' time (until death), ''failure'' time (until hardware fails), etc. To be general, we can just say ''event'' time.

The random variable representing the event time for subject $i$ is typically written $T_i$. Several situations are then possible to define the observations:


<ul>
* The event time is exactly observed.


::[[File:survival1.png|link=]]


: Then, the observation for individual $i$ is $y_i = t_i$, where $t_i$ is a realization of the random variable $T_i$.
<br>

* We may know the event has happened in an interval $I_i$ but not know the exact time $t_i$. This is ''interval censoring''. For example, at a routine check-up, cancer recurrence may be detected, and we only know that it has occurred at some point in time since the last check-up.


::[[File:survival3.png|link=]]


: The observation for individual $i$ is the event: $y_i = $ "$a_i < t_i \leq b_i$".
<br>

* If we assume that the trial ends at time $\tstop$, then the event may happen after the end of the trial period. This is ''right censoring''.


::[[File:survival2.png|link=]]


: There are several variations of this for defining what the observations are:
<br>

* If events (before $\tstop$) are exactly observed, then for $i=1,2,\ldots, N$,

{{Equation1|
equation=<math>
y_i = \left\{
\begin{array}{ll}
t_i & {\rm if \quad} t_i \leq \tstop \\
{\rm t_i > \tstop \quad} & {\rm otherwise. \quad}
\end{array} \right.
</math>}}
</ul>


{{ExampleWithText&Table
|title1=Example:
|title2=
|equation=
Assume  that a trial starts at $\tstart=0$ and ends at $\tstop=5$, and that we obtain the following observations from 4 individuals: 

$y_1 = 3.2$ 

$y_2=$ "$t_2>5$"

$y_3= 2.7$ 

$y_4 =$ "$t_4>5$"


These observations can be stored in a data file as shown in the table on the right.

Here, "event=0" at time $t$ means that the event happened after $t$ while "event=1" means that the event happened at time $t$. 

The lines with $t=0$ are used to state the trial start time $\tstart=0$.

|table=
{{{!}} class="wikitable" align="center" style="width: 75%"
!{{!}} ID {{!}}{{!}} TIME  {{!}}{{!}} EVENT  
{{!}}-
{{!}}1 {{!}}{{!}} 0 {{!}}{{!}} 0 
{{!}}-
{{!}}1 {{!}}{{!}} 3.2 {{!}}{{!}} 1 
{{!}}-
{{!}}2 {{!}}{{!}} 0 {{!}}{{!}} 0 
{{!}}-
{{!}}2 {{!}}{{!}} 5 {{!}}{{!}} 0 
{{!}}-
{{!}}3 {{!}}{{!}} 0 {{!}}{{!}} 0 
{{!}}-
{{!}}3 {{!}}{{!}} 2.7 {{!}}{{!}} 1 
{{!}}-
{{!}}4 {{!}}{{!}} 0 {{!}}{{!}} 0 
{{!}}-
{{!}}4 {{!}}{{!}} 5 {{!}}{{!}} 0 
{{!}}} 
}}


<ul>
* If events before $\tstop$ are interval censored, then for $i=1,2,\ldots, N$,

{{Equation1
|equation=<math>
y_i = \left\{
\begin{array}{ll}
{\rm a_i < t_i \quad \leq \quad b_i} & {\rm if \quad} t_i\leq \tstop \\
{\rm t_i  >  \tstop \quad} & {\rm otherwise.}
\end{array}
\right.
</math> }}
</ul>


{{ExampleWithText&Table
|title1=Example:
|title2=
|equation=
Assume that we have censoring intervals of length 1: 


$(0,1],(1,2],\ldots,(4,5]$.


For the same four individuals as the previous example, we now have the following observations: 


$y_1=$  "$3 < t_1 \leq 4$", 

$y_2=$  "$t_2>5$", 

$y_3=$ "$2< t_3 \leq 3$", 

$y_4=$  "$t_4>5$". 


These observations can be stored in a data file as shown in the table on the right.

Here "event=0" at time $t$ means that the event happened after $t$ while "event=1" means that the event happened before time $t$.
|table=
{{{!}} class="wikitable" align="center" style="width: 75%"
!{{!}} ID {{!}}{{!}} TIME  {{!}}{{!}} EVENT  
{{!}}-
{{!}}1 {{!}}{{!}} 0 {{!}}{{!}} 0 
{{!}}-
{{!}}1 {{!}}{{!}} 3 {{!}}{{!}} 0 
{{!}}-
{{!}}1 {{!}}{{!}} 4 {{!}}{{!}} 1 
{{!}}-
{{!}}2 {{!}}{{!}} 0 {{!}}{{!}} 0 
{{!}}-
{{!}}2 {{!}}{{!}} 5 {{!}}{{!}} 0 
{{!}}-
{{!}}3 {{!}}{{!}} 0 {{!}}{{!}} 0 
{{!}}-
{{!}}3 {{!}}{{!}} 2 {{!}}{{!}} 0 
{{!}}-
{{!}}3 {{!}}{{!}} 3 {{!}}{{!}} 1 
{{!}}-
{{!}}4 {{!}}{{!}} 0 {{!}}{{!}} 0 
{{!}}-
{{!}}4 {{!}}{{!}} 5 {{!}}{{!}} 0 
{{!}}}
}}


<br>

== Probability distributions == 


Several functions play key roles in time-to-event analysis: the [http://en.wikipedia.org/wiki/Survival_function survival function] the [http://en.wikipedia.org/wiki/Survival_analysis#Hazard_function_and_cumulative_hazard_function hazard function] and the [http://en.wikipedia.org/wiki/Survival_analysis#Hazard_function_and_cumulative_hazard_function cumulative hazard function].
We are still working under a population approach here and so these functions, detailed below, are therefore individual functions, i.e., each subject has its own. As we are using parametric models, this means that these functions depend on individual parameters $(\psi_i)$.


<ul>
* The '''survival function''' $S(t; \psi_i)$ gives the probability that the event happens to individual $i$ after time $t>t_{start}$:

{{Equation1
|equation=<math> 
S(t; \psi_i) \ \ \eqdef \ \ \prob{T_i>t ; \psi_i} .
</math> }}



* The '''hazard function''' $\hazard(t;\psi_i)$ is defined for individual $i$ as the instantaneous rate of the event at time $t$, given that the event has not already occurred:

{{Equation1
|equation=<math>  
\hazard(t;\psi_i) \ \ \eqdef \ \  \lim_{dt\to 0} \displaystyle{\frac{S(t;\psi_i) - S(t + dt;\psi_i)}{ S(t;\psi_i) \, dt} }. 
</math> }}

: This is equivalent to: 

{{Equation1
|equation=<div id="HazardSurvival" ><math> 
\hazard(t;\psi_i) \ \ = \ \ -\displaystyle{ \frac{d}{dt} } \log{S(t;\psi_i)}. 
</math></div>
|reference=(1)
}} 


* Another useful quantity is the '''cumulative hazard function''' $\cumhaz(a,b;\psi_i)$,  defined for individual $i$ as:

{{Equation1
|equation=<math>
\cumhaz(a,b;\psi_i) \ \ \eqdef \ \ \displaystyle{\int_a^b \hazard(t;\psi_i) \, dt }.
</math>}}

: Note that [[#HazardSurvival|(1)]] implies that:

{{Equation1
|equation=<math>
S(t;\psi_i) \ \ = \ \ e^{-\cumhaz(t_{start},t;\psi_i)}.
</math> }}
</ul>


Equation [[#HazardSurvival|(1)]] shows that the hazard function $\hazard(t;\psi_i)$ characterizes the problem, because knowing it is the same as knowing the survival function $S(t;\psi_i)$. The probability distribution of survival data is therefore completely defined by the hazard function.
Let $\qcyipsii$ be the conditional distribution of the  observation $y_i$ given the vector of individual parameters $\psi_i$. Its pdf can be easily computed for the various censoring situations discussed above:


<ol>
<li>If the event is exactly observed with $y_i=t_i$, the density is the derivative of the cumulative density function, i.e., the derivative of $1 - S(t_i;\psi_i)$:</li>

{{Equation1
|equation=<math>
\begin{eqnarray}\pcyipsii(y_i {{!}} \psi_i) &=& \frac{d}{dt_i}\left(1 - e^{-\cumhaz(t_{start},t_i;\psi_i)}\right)\\
%&=& \left(\frac{d}{dt_i} \int_{t_{start} }^{t_i} \hazard(u;\psi_i) \, du   \right)  e^{-\cumhaz(t_{start},t_i;\psi_i)}\\
&=&\hazard(t_i;\psi_i)e^{-\cumhaz(t_{start},t_i;\psi_i)} .
\end{eqnarray}</math> }}


<li>If the event is interval-censored with $y_i=\,$  "$a_i<t_i\leq b_i$":</li>

{{Equation1
|equation=<math>\begin{eqnarray}
\pcyipsii(y_i {{!}} \psi_i) &=& \prob{T_i \in (a_i,b_i]\,{{!}} \,\psi_i}  \\
%&=& \prob{T_i \leq b_i {{!}} \psi_i} - \prob{T_i \leq a_i {{!}} \psi_i}  \\
%&=& (1-S( b_i ; \psi_i)) - (1-S( a_i ; \psi_i)) \\
&=&  e^{-\cumhaz(t_{start},a_i;\psi_i)} - e^{-\cumhaz(t_{start},b_i;\psi_i)} .
\end{eqnarray}</math> }}


<li>If the event is right-censored with $y_i= \,$ "$t_i>t_{stop}$":</li>

{{Equation1
|equation=<math>\begin{eqnarray}
\pcyipsii(y_i {{!}} \psi_i) &=& \prob{T_i > t_{stop}   {{!}} \psi_i}  \\
%&=& S( t_{stop} ; \psi_i) \\
&=&  e^{-\cumhaz(t_{start},t_{stop};\psi_i)} .
\end{eqnarray}</math> }}
</ol>


<br><br>

==Repeated events==



Sometimes, an event can potentially happen again and again, e.g., [http://en.wikipedia.org/wiki/Epileptic_seizure epileptic seizures], heart attacks,  etc.
For any given hazard function $\hazard$, the survival function $S$ for individual $i$  now represents survival since the previous event  at $t_{i,j-1}$, written here in terms of the cumulative hazard from $t_{i,j-1}$ to $t_{i,j}$:

{{Equation1
|equation=<math>\begin{eqnarray}
S(t_{i,j} {{!}} t_{i,j-1};\psi_i) &=& \prob{T_{i,j} > t_{i,j}\, {{!}} \,T_{i,j-1} = t_{i,j-1};\psi_i} \\
&=&  e^{-\cumhaz(t_{i,j-1},t_{i,j};\psi_i)}  \\
&=&  \exp\left({-\int_{t_{i,j-1} }^{t_{i,j} } \hazard(t;\psi_i) \, dt}\right) .
\end{eqnarray}</math> }}

<!--%In the most simple case, $y_i$ is a vector of known event times: $y_i = (t_{i1},t_{i2},\ldots,t_{i\,n_i}).$ -->


<br>

==Censoring and probability distributions==


Taking into account censoring for repeated events is slightly more complicated than for one-off events.
First, let us assume that a trial starts at time $t_{start}$ and ends at time $t_{stop}$. Let $(T_{i1}, T_{i2}, \ldots )$ be random event times after $t_{start}$. Then, we can distinguish between the two following situations:



<ul>
1. ''Exactly observed events:'' A sequence of $n_i$ event times  is precisely observed before $t_{stop}$, i.e., ${\rm y_i = (t_{i,1},t_{i,2},\ldots,t_{i,n_i}, \quad t_{i,n_i+1}>\tstop)}$. 

: The conditional pdf of $y_i$ is given by:

{{EquationWithRef
|equation=<div id="repeatcensor" ><math> 
\pcyipsii(y_i {{!}} \psi_i) = \left(\prod_{j=1}^{n_i}\hazard(t_{ij};\psi_i)e^{-\cumhaz(t_{i,j-1},t_{i,j};\psi_i)}  \right)e^{-\cumhaz(t_{n_i},\tstop;\psi_i)} ,
</math></div>
|reference=(1) }}

: where $t_{i0}=\tstart$.
</ul>

{{ExampleWith2Tables
|title1=Example
|title2=
|text=
Suppose that for individual $i=1$ we know there were 8 events but only 7 of them occurred before $\tstop$. Here is a graphic showing the events that were exactly observed:


::[[File:survival4.png|link=]]


This data is then stored in the table on the left below. We see that the 8th and final event is noted "event = 0" with time $\tstop = 18$, indicating that the event was not observed at the end of the time period $\tstop$. In the table on the right, we show the contributions of each observation to the conditional pdf of $y_1$. Indeed, equation [[#repeatcensor|(1)]] means that the pdf of $y_1=(y_{1,1}, \ldots, y_{1,8})$ is the product of the conditional pdfs given in the right table.


|table1=
{{{!}} class="wikitable"  align="center" style="width:120%; margin-left:10%;margin-right:10%"
!{{!}} ID  {{!}}{{!}} TIME  {{!}}{{!}} EVENT  
{{!}}-
{{!}}1 {{!}}{{!}} 0 {{!}}{{!}} 0 
{{!}}-
{{!}}1 {{!}}{{!}} 1.4 {{!}}{{!}} 1 
{{!}}-
{{!}}1 {{!}}{{!}} 3.5 {{!}}{{!}} 1 
{{!}}-
{{!}}1 {{!}}{{!}} 4.4 {{!}}{{!}} 1 
{{!}}-
{{!}}1 {{!}}{{!}} 5.6 {{!}}{{!}} 1 
{{!}}-
{{!}}1 {{!}}{{!}} 9.7 {{!}}{{!}} 1 
{{!}}-
{{!}}1 {{!}}{{!}} 11.4 {{!}}{{!}} 1 
{{!}}-
{{!}}1 {{!}}{{!}} 15.8 {{!}}{{!}} 1 
{{!}}-
{{!}}1 {{!}}{{!}} 18 {{!}}{{!}} 0 
{{!}}}

|table2 =
{{{!}} class="wikitable" align="center" style="width:200%; margin-right:10%; margin-left:10%"
!{{!}} pdf 
{{!}}-
{{!}} 1 
{{!}}-
{{!}} $\hazard(1.4;\psi_1)e^{-\cumhaz(0,1.4;\psi_1)}$
{{!}}-
{{!}} $\hazard(3.5;\psi_1)e^{-\cumhaz(1.4,3.5;\psi_1)}$
{{!}}-
{{!}} $\hazard(4.4;\psi_1)e^{-\cumhaz(3.5,4.4;\psi_1)}$
{{!}}-
{{!}} $\hazard(5.6;\psi_1)e^{-\cumhaz(4.4,5.6;\psi_1)}$
{{!}}-
{{!}} $\hazard(9.7;\psi_1)e^{-\cumhaz(5.6,9.7;\psi_1)}$
{{!}}-
{{!}} $\hazard(11.4;\psi_1)e^{-\cumhaz(9.7,11.4;\psi_1)}$
{{!}}-
{{!}} $\hazard(15.8;\psi_1)e^{-\cumhaz(11.4,15.8;\psi_1)}$
{{!}}-
{{!}} $e^{-\cumhaz(15,18;\psi_1)}$
{{!}}}
}}


<ul>
2. ''Interval-censored events:'' Let $(b_{0}, b_1], (b_{1}, b_2], \ldots , (b_{K-1}, b_K]$ be a sequence of successive intervals with $\tstart=b_0<b_1<b_2 < \ldots <b_K  = \tstop$. We do not know the exact event times, but a sequence $(m_{ik}; \, 1 \leq k  \leq  K)$ is observed, where $m_{ik}$ is the number of events that occurred for individual $i$ in interval  $(b_{k-1}, b_k]$.

: We can show that the conditional pdf of $y_i$ is given by:

{{EquationWithRef
|equation=<div id="pdf_mult_int" ><math>
\pcyipsii(y_i {{!}} \psi_i) =  \prod_{k=1}^{K} e^{-\cumhaz(b_{k-1}, b_k;\psi_i)} \displaystyle{\frac{\cumhaz^{m_{ik} }(b_{k-1}, b_k;\psi_i)}{m_{ik}!} } .
</math></div>
|reference=(2) }}

: In other words, the number of events per interval for individual $i$ is a (possibly non-homogeneous) Poisson process with intensity $\cumhaz(b_{k-1}, b_k;\psi_i)$ in interval $(b_{k-1}, b_k]$.


{{ExampleWith2Tables
|title1=Example
|title2=

|text= Here is a graphic that shows an example of the interval boundaries and the number of events that occurred in each interval for individual $i=1$.


::[[File:survival5.png|link=]]


The table on the left below shows the same data. Using [[#pdf_mult_int|(2)]] we see that the conditional pdf of $y_1=(y_{1,1}, \ldots, y_{1,6})$ is the product of the conditional pdfs given in the table on the right.


|table1=
{{{!}} class="wikitable" align="center" style="width:120%; margin-left:10%;margin-right:10%"
!{{!}} ID  {{!}}{{!}} TIME  {{!}}{{!}} EVENT  
{{!}}-
{{!}}1 {{!}}{{!}} 0 {{!}}{{!}} 0 
{{!}}-
{{!}}1 {{!}}{{!}} 3 {{!}}{{!}} 1 
{{!}}-
{{!}}1 {{!}}{{!}} 6 {{!}}{{!}} 3 
{{!}}-
{{!}}1 {{!}}{{!}} 9 {{!}}{{!}} 0 
{{!}}-
{{!}}1 {{!}}{{!}} 12 {{!}}{{!}} 2 
{{!}}-
{{!}}1 {{!}}{{!}} 15 {{!}}{{!}} 0 
{{!}}-
{{!}}1 {{!}}{{!}} 18 {{!}}{{!}} 1 
{{!}}}

|table2=
{{{!}} class="wikitable" align="center" style="width:200%; margin-right:10%; margin-left:10% "
!{{!}} pdf 
{{!}}-
{{!}} 1 
{{!}}-
{{!}} $e^{-\cumhaz(0,3;\psi_1)}\cumhaz(0,3;\psi_1)  $
{{!}}-
{{!}} $e^{-\cumhaz(3,6;\psi_1)} {\cumhaz^{3}(3,6;\psi_1)}/{6}  $
{{!}}-
{{!}} $e^{-\cumhaz(6,9;\psi_1)}$
{{!}}-
{{!}} $e^{-\cumhaz(9,12;\psi_1)} {\cumhaz^{2}(9,12;\psi_1)}/{2}  $
{{!}}-
{{!}} $e^{-\cumhaz(12,15;\psi_1)}$
{{!}}-
{{!}} $e^{-\cumhaz(15,18;\psi_1)}\cumhaz(15,18;\psi_1)  $
{{!}}}
}}


{{Remarks
|title=Remark
|text= if the total number $n_i$ of (observed and unobserved) events for individual $i$ is known to be finite, then formula [[#pdf_mult_int|(2)]] is slightly modified when the last event occurs before $\tstop$ ($t_{n_i}<\tstop$).
Assume that the last event for individual $i$ occurs in the $K_i$-th interval.  Let $s_{i} = \sum_{i=1}^{k_i-1} m_{ik}$ be the number of events that occurred before this interval. Then, we can show that

{{EquationWithRef_Special
|equation=<div id="pdf_mult_int2"><math>
\pcyipsii(y_i {{!}} \psi_i) = \prod_{k=1}^{K_i-1} \left( \displaystyle{ \frac{\cumhaz^{m_{ik} }(b_{k-1}, b_k;\psi_i)}{m_{ik}!}  }e^{-\cumhaz(b_{k-1}, b_k;\psi_i)} \right)
\!\times \!\left(1 - \sum_{\ell=0}^{n_i-s_{i} }   \displaystyle{ \frac{\cumhaz^{\ell}(b_{k_i -1},b_{k_i};\psi_i)}{\ell!} } e^{-\cumhaz(b_{k_i -1},b_{k_i};\psi_i)}\right) . </math></div>
|reference=(3) }}
}}


<br>

== Examples of hazard functions==



<ul>
* ''Constant hazard model:'' 
: The most simple case is that of a constant hazard function: $\hazard(t;\psi_i) = \hazard_i \in \Rset$. Here, $\psi_i=\hazard_i$. 
<br>


* ''Proportional hazards model:''

{{Equation1
|equation=<math>
\hazard(t;\psi_i) = \hazard_0(t;\alpha_i) \, e^{  \langle \beta , c_i  \rangle}.
</math>}}

: Here, the hazard is decomposed into two terms: a baseline function $\hazard_0$ of $t$, and an "individual" term, function of some individual covariates $c_i$. $ \langle \beta , c_i  \rangle$ means a scalar product, i.e., a linear function of $c_i$. In a proportional hazards model, a unit increase in the value of a covariate has a multiplicative effect on the hazard.

: In the usual proportional hazard model, $\alpha_i$ is a population constant ($\alpha_i=\alpha$). Then, $\psi_i$ can be decomposed into a set of population parameters $\alpha$ and an individual parameter $ \langle \beta , c_i  \rangle$. A straightforward extension consists in assuming that $\alpha_i$ is also an individual parameter.
<br>


* ''Extended proportional hazards model:''

: Another possible extension assumes that the hazard function is a (possibly nonlinear) function $u$ of  a regression variable $x_i$:

{{Equation1
|equation=<math>
\hazard(t;\bpsi_i) = \hazard_0(t;\alpha_{i}) \, e^{  u(\beta_i,x_i(t))} .
</math> }}

:Consider for example that $x_i(t)$ is the plasmatic concentration of a drug at time $t$ for individual $i$. Then, $u(\beta_i,x_i(t))$ is the term that represents (i.e., models) the effect of the drug  on the hazard, while $\hazard_0(t;\alpha_i)$  might model the effect of disease progression on the hazard.
<!--%We consider here parametric functions that possibly depend on individual parameters.-->

: In this example, $x_i(t)$ is the "true" plasmatic concentration for subject $i$ at time $t$, and it is a continuous function of time. However, in practice it is only measured at precise times, so a longitudinal model for plasmatic concentration is needed to give a concentration value for each $t$.
:Therefore, in practice we need to develop a ''joint model'' in order to simultaneously model time-to-events data and longitudinal data. Such an approach is introduced in the [[Joint models]] section.
<br>
 

* ''Accelerated failure time (AFT) model:''

:Unlike proportional hazards models, the AFT model supposes that a change in a covariate has a multiplicative effect not on the hazard but the ''predicted event time''. This can be written as:

{{Equation1
|equation=<math>
\log(T_i) =  \langle \psi_i , c_i  \rangle + \xi_i
</math>
}}

: where $\xi_i$ is a zero-mean random variable, e.g., a centered normal distribution. Usually, parameters are fixed effects: $\psi_i=\psi$ for each subject $i$.
: To calculate the hazard function, let us first denote  $p_{\xi_i}$ the density and $F_{\xi_i}$ the cdf of $\xi_i$, and to simplify, denote $\mu_i = \langle \psi_i , c_i  \rangle$ the mean of $\log(T_i)$. We begin by calculating the survival function:

{{Equation1
|equation=<math>\begin{eqnarray}
S(t;\psi_i) &=& \prob{\log{T_i} > \log{t} ; \bpsi_i} \\
&=& \int_{\log{t}-\mu_i}^{\infty} p_{\xi_i}(u; \psi_i) \, du \\
&=& 1 - F_{\xi_i}(\log{t}-\mu_i ; \psi_i) .
\end{eqnarray}</math> }}

:Calculating [[#HazardSurvival|(1)]] then gives the hazard function:

{{Equation1
|equation=<math>
\hazard(t;\psi_i) = \displaystyle{ \frac{p_{\xi_i}(\log{t} - \mu_i; \psi_i)}{t(1- F_{\xi_i}(\log{t} - \mu_i; \psi_i))} }\,
</math> }}

<br><br>
-------
<br><br>

{{Summary 
|title=Summary
|text=
For a given vector of individual parameters $\psi_i$, a model for (repeated)  time-to-event data is completely defined by


<ol>
<li> the hazard function $\hazard(t ; \psi_i)$, or the survival function $S(t ; \psi_i)$ </li>

<li> (possibly) the interval and/or right censoring process </li>

<li> (possibly) the maximum number of possible events </li>
</ol> }}


<br>

<!--
==$\mlxtran$ for time-to-event data models==
-->

<br>

==Bibliography==

<bibtex>
@book{aalen2008,
author = {Aalen, O. and Borgan, O. and Gjessing, H.},
title  = {Survival and Event History Analysis. },
publisher = {Springer},
address = {New York},
year = {2008}
}
</bibtex>
<bibtex>
@book{andersen2006survival,
  title={Survival analysis},
  author={Andersen, P. K.},
  year={2006},
  publisher={Wiley Online Library}
}
</bibtex>
<bibtex>
@article{diggle1994,
author = {Diggle, P. and Kenward, M. G.},
title  = {Informative drop-out in longitudinal data analysis.},
journal = {Appl. Stats},
volume = {43},
number = {},
pages = {49-93},
year = {1994}
}
</bibtex>

<bibtex>
@book{duchateau2008,
author = {Duchateau, L. and Janssen, P.},
title  = {The Frailty Model. Statistics for Biology and Health },
publisher = {Springer.},
volume = {},
pages = {},
year = {2008},
series = {},
address = {New York},
edition = {},
month = {}
}
</bibtex>
<bibtex>
@book{fleming2011counting,
  title={Counting processes and survival analysis},
  author={Fleming, T. R. and Harrington, D. P.},
  volume={169},
  year={2011},
  publisher={Wiley}
}
</bibtex>
<bibtex>
@article{huang2007,
author = {Huang, X. and Liu, L.},
title  = {A joint frailty model for survival and gap times between recurrent events.},
journal = {Biometrics},
volume = {63},
number = {},
pages = {389-397},
year = {2007}
}
</bibtex>
<bibtex>
@book{ibrahim2005bayesian,
  title={Bayesian survival analysis},
  author={Ibrahim, J. G. and Chen, M.-H. and Sinha, D.},
  year={2005},
  publisher={Wiley Online Library}
}
</bibtex>
<bibtex>
@book{kalbfleisch2011statistical,
  title={The statistical analysis of failure time data},
  author={Kalbfleisch, J. D. and Prentice, R. L.},
  year={2011},
  publisher={Wiley-Interscience}
}
</bibtex>
<bibtex>
@article{kelly2000,
author = {Kelly, P. J. and Jim, L. L.},
title  = {Survival analysis for recurrent event data: an application to childhood infectious disease.},
journal = {Statistics in Medicine},
volume = {19},
number = {1},
pages = {13-33},
year = {2000}
}
</bibtex>
<bibtex>
@book{klein2003survival,
  title={Survival analysis: techniques for censored and truncated data},
  author={Klein, J. P. and Moeschberger, M. L.},
  year={2003},
  publisher={Springer}
}
</bibtex>
<bibtex>
@book{klein1997,
author = {Klein, J. P. and Moeschberger, M. L.},
title  = { Survival Analysis - Techniques for Censored and Truncated Data. },
publisher = {Springer-Verlag},
volume = {},
pages = {},
year = {1997},
series = {},
address = {New York},
edition = {},
month = {}
}
</bibtex>
<bibtex>
@book{kleinbaum2011survival,
  title={Survival analysis},
  author={Kleinbaum, D. G.},
  year={2011},
  publisher={Springer}
}
</bibtex>
<bibtex>
@book{littell2006sas,
  title={SAS for mixed models},
  author={Littell, R. C.},
  year={2006},
  publisher={SAS institute}
}
</bibtex>
<bibtex>
@book{miller2011survival,
  title={Survival analysis},
  author={Miller Jr, R. G.},
  year={2011},
  publisher={Wiley-Interscience}
}
</bibtex>
<bibtex>
@book{wienke2010frailty,
  title={Frailty models in survival analysis},
  author={Wienke, A.},
  volume={37},
  year={2010},
  publisher={Chapman & Hall}
}
</bibtex>



{{Back&Next
|linkBack=Model for categorical data
|linkNext=Joint models }}