new_DIRICHLET now accept to have a column with zero probability. If only one column have probability, returns 1 in that column and zero elsewhere
new_BETA_lci and new_BETA_lci2 are separated functions. The first keep the expected value but could be less accurate on CI, the second could be more accurate on the CI but the expected value may change. The second use a ML function to find the parameters. Equivalent changes has made to fitbeta and fitbeta_ml