beta2m <- function(b) log2(b/(1-b))