Re: [R-sig-Geo] How to fit a spatial hurdle model with the INLA package ?

2017-02-17 Thread 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 avai

Re: [R-sig-Geo] How to fit a spatial hurdle model with the INLA package ?

2017-02-17 Thread 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 somethi

Re: [R-sig-Geo] How to fit a spatial hurdle model with the INLA package ?

2017-02-17 Thread 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 Onkel

Re: [R-sig-Geo] How to fit a spatial hurdle model with the INLA package ?

2017-02-17 Thread Thierry Onkelinx
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

[R-sig-Geo] How to fit a spatial hurdle model with the INLA package ?

2017-02-17 Thread lancelot
Dear all, I am looking for an example (dataset and code) of fitting a spatial hurdle model (zero-inflated model with a single source of zeros) using the INLA package for R. With this approach, two models are jointly fitted: (1) the probability to observe at least one event (the single source