Sounds like you could also use a Multinomial distribution as likelihood?
Otherwise, you might have overflow issue with k, range from 0 to 1,000,000 - you can try scaling k by some constant (i think the Poisson likelihood still works)
1 Like