calculating the value of the likelihood L(Y | K = 0)
prod_1_k0 <- matrix(0,dias,1)
prod_2_k0 <- matrix(0,dias,1)
prod_k0 <- matrix(0,dias,1)
for(i in 1:dias){
for(j in 1:2){
prod_1_k0[i] = (gamma(count_init_0[i,j]+alpha_init_0[i,j])/gamma(alpha_init_0[i,j])) }
prod_2_k0[i] = (gamma(sum(alpha_init_0[i,]))/gamma(sum(count_init_0[i,] +
 alpha_init_0[i,])))
prod_k0 [i] = prod_1_k0[i]*prod_2_k0[i]}
for(i in 1:dias){
term_general_like[i, 1] = prod_k0[i]}
#######