term corresponding to the initial distribution
```

```
prod_1_k1 = 1
for(j in 1:2){
prod_1_k1 = prod_1_k1*((gamma(count_init_1[j]+alpha_init_1[j])/
gamma(alpha_init_1[j])))*((gamma(sum(alpha_init_1))/
gamma(sum(count_init_1+ alpha_init_1))))}
term_init_like[2] = prod_1_k1