R/adjustKRW.R

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

adjustKRW <-
function (K, acceptRate, q) {
    low <- 0.4
    upp <- 0.7
    if (acceptRate >= low && acceptRate <= upp)
        return(K)
    if (acceptRate > upp)
        K <- max(q, K * 0.7)
    if (acceptRate < low)
        K <- max(q, K * 1.2)
    K
}