added a new function to get figure for BIC model selection with beta's on different scales do_randB_scl.r