added faster implementation based on direct computation of multinomial likelihood