How to fit a spatial hurdle model with the INLA package ?
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" <renaud.lancelot at cirad.fr>: 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
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
<mailto: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 <tel:%2B33%204%2067%2059%2037%2017> - Fax
+33 4 67 59 37 98 <tel:%2B33%204%2067%2059%2037%2098>
Secr. +33 4 67 59 37 37 <tel:%2B33%204%2067%2059%2037%2037> - Cell.
+33 6 77 52 08 69 <tel:%2B33%206%2077%2052%2008%2069>
Twitter: @RenaudLancelot
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org <mailto:R-sig-Geo at r-project.org>
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
<https://stat.ethz.ch/mailman/listinfo/r-sig-geo>
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
[[alternative HTML version deleted]]