x <- ifelse(x < 0 | x > 1, NA, x ) if( any(is.na(x)) ) warning('x not between 0 and 1') log(x / (1 - x) )