Multinomial returns nan

[junpenglao]

@Dominik

However one could achieve the same result by removing this category from the vector a and append x_i=0 to the result afterwards.

This is a good point. Would you like to send a pull request on Github to fix this?