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

[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" <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]]