How to fit a spatial hurdle model with the INLA package ?
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 To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of. ~ Sir Ronald Aylmer Fisher The plural of anecdote is not data. ~ Roger Brinner The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. ~ John Tukey 2017-02-17 13:10 GMT+01:00 lancelot <renaud.lancelot at cirad.fr>:
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 of zeros) (2) the (strictly > 0) count of events The spatial correlation might be represented by a conditional auto-regressive model. I am aware of the R-!NLA documentation available at http://www.r-inla.org/ All the best Renaud -- Renaud Lancelot UMR Animaux, Sant?, Territoires, Risques et ?cosyst?mes Animals, Health, Territories, Risks & Ecosystems (ASTRE) Adjoint au Directeur / Deputy director http://umr-astre.cirad.fr/ CIRAD, UMR ASTRE, Campus International de Baillarguet TA A-DIR / B F34398 Montpellier Tel. +33 4 67 59 37 17 - Fax +33 4 67 59 37 98 Secr. +33 4 67 59 37 37 - Cell. +33 6 77 52 08 69 Twitter: @RenaudLancelot
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo