R/yj_adjust_lambda.R

In mdmb: Model Based Treatment of Missing Data

## File Name: yj_adjust_lambda.R
## File Version: 0.18

yj_adjust_lambda <- function( lambda, lambda0, eps=1E-3 )
{
    if (is.na(lambda) ){
        lambda <- lambda0 + 2*eps
    }
    if ( abs( lambda - lambda0 ) < eps ){
        sign_lambda <- 1 * ( sign(lambda) >=0  )
        lambda <- lambda0 +  sign_lambda * eps
    }
    return(lambda)
}