calculation of the initial distribution
prob_init <- matrix(0, 2^K, 1)
for(m in 1: 2^K){
prob_init[m] = (alpha_init_k[m] + count_init_k[m] - 1)/somainit}
write.csv(prob_init, file = "prob_init_chain.txt")