term corresponding to the rest of the days (rows)
prod_2_k1 <- matrix(0,(dias - 1),1)
prod_3_k1 <- matrix(0,(dias - 1),1)
prod_k1 <- matrix(0,(dias - 1),1)
for(i in 1:(dias - 1)){
for(j in 1:2){
prod_2_k1[i] = (gamma(count_0_k1[i,j] + alpha_0_k1[i,j])/gamma(alpha_0_k1[i,j]))*
(gamma(count_1_k1[i,j] + alpha_1_k1[i,j])/gamma(alpha_1_k1[i,j]))}
for(j in 1:2){
prod_3_k1[i] = (gamma(sum(alpha_0_k1[i,])+1)*gamma(sum(alpha_1_k1[i,])))/
(gamma(sum(count_0_k1[i,]+alpha_0_k1[i,])+sum(count_1_k1[i,]+alpha_1_k1[i,])+1))}
prod_k1[i] = prod_2_k1[i] * prod_3_k1[i]}
for(i in 1:(dias-1)){
term_general_like[i, 2] = prod_k1[i]}
####