R/extractInits.R

In JMbayes: Joint Modeling of Longitudinal and Time-to-Event Data under a Bayesian Approach

extractInits <-
function (object) {
    if (!inherits(object, "JMbayes"))
        stop("Use only with 'JMbayes' objects.\n")
    init <- object$postMeans
    if (!is.null(init$sigma)) {
        init$tau <- 1 / object$postMeans$sigma^2
        init <- init[names(init) != "sigma"]
    } else {
        init$tau <- NA
    }
    init$invD <- solve(object$postMeans$D)
    init
}