For BRMS zero-inflated / hurdle models, what does the output of `fitted` produce?


Are the values generated assuming we’ve “jumped” the hurdle? Or is the hurdle probability multiplied first? I couldn’t find anything in the documentation on this – apologies if I missed it.


By default, fitted() returns the predicted mean of the response distribution which is a combination of both parts of a hurdle model. If you want it to return only one of the two parts, use argument dpar = "mu" or `dpar = “hu”.