https://github.com/cran/rstpm2
Tip revision: 4f6cced2730ffce02a21c60576beee964f690a74 authored by Mark Clements on 29 May 2018, 12:45:06 UTC
version 1.4.2
version 1.4.2
Tip revision: 4f6cced
aft.tex
\documentclass{article}
\usepackage{amsmath}
\begin{document}
\section{Smooth accelerated failure time models}
\subsection{Time-constant acceleration factors}
Let survival to time $t$ for covariates $x$ be modelled as an accelerated failure time model using
\begin{align*}
S(t|x) &= S_0(t \exp(-\eta(x;\beta)))
\end{align*}
where $\eta$ is a linear predictor. Moreover, let the baseline survival function be modelled as
\begin{align*}
S_0(t) &= \exp(-\exp(\eta_0(\log(t);\beta_0)))
\end{align*}
where $\eta_0$ is a linear predictor. Then the combined regression model is
\begin{align*}
S(t|x) &= \exp(-\exp(\eta_0(\log(t) -\eta(x; \beta);\beta_0)))
\end{align*}
We can calculate the hazard, such that
\begin{align*}
h(t|x) &= \frac{\partial}{\partial t} \left(-\log(S(t|x))\right) \\
&= \exp(\eta_0(\log(t) -\eta(x; \beta);\beta_0)) \eta_0'(\log(t) -\eta(x; \beta);\beta_0)/t
\end{align*}
\subsection{Time-dependent acceleration factors}
We can model survival as an accelerated failure time model with time-dependent effects as
\begin{align*}
S(t|x) &= S_0\left(\int_0^t\exp(-\eta_1(x,u;\beta)) \mathrm{d}u\right) = S_0\left(t\exp(-\eta(x,t;\beta))\right)
\end{align*}
for a time-specific linear predictor $\eta_1$ and where $\eta$ now models for cumulative time-dependent effects.
By differentiation with respect to time $t$, we have that
\begin{align*}
\exp(-\eta_1(x,t;\beta)) &= \frac{\partial}{\partial t} t\exp(-\eta(x,t;\beta)) \\
&= \exp(-\eta(x,t;\beta)) \left(1- t \frac{\partial}{\partial t} \eta(x,t;\beta)\right)
\end{align*}
so that $\eta_1(x,t;\beta)=\eta(x,t;\beta) - \log\left(1- t \frac{\partial}{\partial t} \eta(x,t;\beta)\right)$. This shows that we can recover the time-specific acceleration factors from $\eta$.
The combined regression model is then
\begin{align*}
S(t|x) &= \exp(-\exp(\eta_0(\log(t) -\eta(x,t; \beta);\beta_0)))
\end{align*}
The hazard is then
\begin{align*}
h(t|x) &= \frac{\partial}{\partial t} -\log(S(t|x)) \\
&= \exp(\eta_0(\log(t) -\eta(x,t; \beta);\beta_0))\, \eta_0'(\log(t) -\eta(x,t; \beta);\beta_0)\times \\
&\qquad \left(1/t-\frac{\partial}{\partial t}\eta(x,t; \beta)\right)
\end{align*}
\subsection{R implementation}
The linear predictor $\eta_0$ is modelled using natural splines. The linear predictor $\eta$ can be modelled freely provided the linear predictor is a smooth function of time.
Initial values for the time-constant log acceleration factors were calculated from a Weibull regression. The time-varying log acceleration factors were assumed to be zero. Cox regression and the Breslow estimator was used to calculate baseline survival; parameters for the baseline were estimated using linear regression with natural splines fitted to the log-times at the event times.
We included quadratic penalties to ensure that $\eta_0'(\log(t) -\eta(x,t; \beta);\beta_0)$ and $1/t-\frac{\partial}{\partial t}\eta(x,t; \beta)$ were positive. For speed, the model has been implemented in C++, with optimisation using the Nelder-Mead algorithm.
\end{document}
%%% Local Variables:
%%% mode: latex
%%% TeX-master: t
%%% End: