The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.
The R package JLPM implements in the jointLPM function the estimation of a joint shared random effects model.
The longitudinal data (\(y_1, y_2, \dots, y_K\)) can be continuous or ordinal and are modelled using a mixed model.
For the continuous Gaussian case:
\[ \forall k \in 1,\dots,K \hspace{1cm} y_k(t_{ijk}) = X_i(t_{ijk})\beta + Z_i(t_{ijk})u_i + \varepsilon_{ijk} \]
For the continuous non Gaussien case, a transformation \(H_k\) is estimated for each outcome:
\[ \forall k \in 1,\dots,K \hspace{1cm} H_k(y_k(t_{ijk}), \eta_k) = X_i(t_{ijk})\beta + Z_i(t_{ijk})u_i + \varepsilon_{ijk} \]
For the ordinal case, mixed models are combined to the Item Response Theory:
\[ \forall k \in 1,\dots,K \hspace{1cm} \mathbb{P}( y_k(t_{ijk}) = m) \Leftrightarrow \mathbb{P}(\eta_{m-1} < X_i(t_{ijk})\beta + Z_i(t_{ijk})u_i + \varepsilon_{ijk} < \eta_m) \]
where \(X_i(t_{ijk})\) and \(Z_i(t_{ijk})\) are vectors of covariates measured at time \(t_{ijk}\) for subject \(i\), \(\beta\) is the vector of fixed effects, \(u_i \sim \mathcal{N}(0,B)\), \(\varepsilon_{ijk} \sim \mathcal{N}(0,\sigma_k^2)\), \(\eta_k\) are the parameters of the link function \(H_k\) or the thresholds associated to the outcome \(y_k\).
Note that even with multiple longitudinal outcomes, a univariate mixed model is estimated.
The time-to-event data are modelled in a proportional hazard model where different associations between the longitudinal outcome and the event can be included.
With an association through the random effects of the longitudinal model:
\[ \alpha_i(t) = \alpha_0(t) \exp(\tilde{X}_i \gamma + \delta u_i) \]
With an association through the current level of the longitudinal model:
\[ \alpha_i(t) = \alpha_0(t) \exp(\tilde{X}_i \gamma + \delta (X_i(t)\beta + Z_i(t)u_i)) \]
With an association through the current slope of the longitudinal model:
\[ \alpha_i(t) = \alpha_0(t, \omega) \exp(\tilde{X}_i \gamma + \delta ( \frac{d}{dt}X_i(t)\beta + \frac{d}{dt}Z_i(t)u_i)) \]
with \(\alpha_0(t, \omega)\) the baseline risk function at time \(t\), parameterized with \(\omega\),\(\tilde{X}_i\) a vector of covariates, \(\gamma\) the fixed effects, \(\delta\) the association parameter.
The parameters \(\beta, B, \sigma_k, \eta_k, \omega, \gamma, \delta\) are estimated using a Marquardt-Levenberg algorithm.
These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.