Hi Tiago,
Yeah, it’s actually a good thing it crashes: the concentration vector of the Dirichlet has to be strictly positive, but you entered zeros at the second and third positions, so the log probability is -inf, making the model impossible to sample from.
Here is a post about Dirichlet and Multinomial distributions by @clausherther – that should help you
Hope this helps