I asked the question on R-SIG-GEO and got two useful replies:

1. From Thierry Onkelinx

That allows for a hurdle model with an intercept only model for the zeroes. Your initial question suggests that you want a more complex model for the zeroes.

Op 17 feb. 2017 5:38 p.m. schreef "lancelot" :

Dear Thierry,

Thank you for your reply. I wonder whether it would be possible to do something with zero-inflated models of type 0: see http://www.math.ntnu.no/inla/r-inla.org/doc/likelihood/zeroinflated.pdf, first equation on page 1.

All the best

Renaud

Le 17/02/2017 à 15:29, Thierry Onkelinx a écrit :

Dear Renaud,

IMHO you can't. INLA currently only fits zero-inflated distributions

(with a single zero-inflation parameter). A hurdle model would require

to fit a zero-truncated distribution.

Best regards,

ir. Thierry Onkelinx

Instituut voor natuur- en bosonderzoek / Research Institute for Nature

and Forest

team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance

Kliniekstraat 25

1070 Anderlecht

Belgium

2. From Facundo Muñoz

 Are you aware of Quiroz et al. 2015? I think it is what you need, or

very close. And they provide data and code.

There is also Muñoz et al. 2016 (Disclaimer, that's me), which uses a

joint model Bernoulli for zeros and Gamma for >0, with spatial-temporal

effect via SPDE. The data and code are available as a companion R

package at https://github.com/famuvie/2016_RisingAshes.

Hope it helps.

ƒacu.-

Quiroz, Z. C., M. O. Prates, and H. Rue (2015, March). A bayesian

approach to estimate the biomass of anchovies off the coast of perú.

/Biometrics/ /71 /(1), 208-217. http://dx.doi.org/10.1111/biom.12227

Muñoz, F., B. Marçais, J. Dufour, and A. Dowkiw (2016, December). Rising

out of the ashes: Additive genetic variation for crown and collar

resistance to hymenoscyphus fraxineus in fraxinus excelsior.

/Phytopathology/ /106 /(12), 1535-1543.

http://dx.doi.org/10.1094/phyto-11-15-0284-r