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”.