matrices to store the values of the likelihood for each order and day
##########
prod_init <- matrix(0, 1, M)
for(K in 2:M){
count_init_k <- matrix(0, 1, (2^K))
mbase_init <- matrix(0, 1, anos)
alpha_init_k <- matrix(0, 1, 2^K)