R/esp.tr.R

In JRM: Joint Regression Modelling

esp.tr <- function(vrb.st, margin){
 
 
mub <- c("BE") 
 

if( !(margin %in% mub) ){
 
   vrb.st <- ifelse( vrb.st > 20,   20, vrb.st )  # it was 28
   vrb.st <- ifelse( vrb.st < -14.5, -14.5, vrb.st ) 
   
   vrb    <- exp(vrb.st)
    
}


if( margin %in% mub ){
 
   vrb.st <- ifelse( vrb.st >  16,  16, vrb.st )  
   vrb.st <- ifelse( vrb.st < -16, -16, vrb.st ) 
   vrb    <- plogis(vrb.st)
    
}


    
    
 list(vrb = vrb, vrb.st = vrb.st )  
 
}