Modello di effetti misti binomiali negativi a zero inflazionati in R

17

Esiste un tale pacchetto che prevede la stima del modello di effetti misti binomiali negativi a zero inflazionati in R?

Con ciò intendo:

Inflazione zero in cui è possibile specificare il modello binomiale per l'inflazione zero, come nella funzione zeroinfl nel pacchetto pscl:
```
zeroinfl (y ~ X | Z, dist = "negbin")
```
dove Z è la formula per il modello di inflazione zero;
Distribuzione binomiale negativa per la parte di conteggio del modello;
Effetti casuali specificati simili alla funzione lmer del pacchetto lme4.

Capisco che glmmADMB può fare tutto ciò, tranne che la formula per l'inflazione zero non può essere specificata (è solo un'intercettazione, cioè Z è solo 1). Ma ci sono altri pacchetti che possono fare tutto?

Sarò molto grato per il tuo aiuto!

— Nikita Samoylov
fonte

Sto cercando la stessa cosa. Penso che glmmADMB sia quello che voglio, ma non riesco a farlo funzionare.

— gregmacfarlane,

Mi chiedo se i pacchetti ZIM o aod possano fare quello che vuoi fare?

— Graeme Walsh,

2

Come aggiornamento, il pacchetto glmmTMB di Ben Bolker supporta un modello misto lineare generalizzato a gonfiaggio zero (ZIGLMM).

— JWilliman,

5

Penso che questo sia il pacchetto che ti serve: glmmADMB. L'ho scaricato qui: http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html

Ma avevo ancora alcuni problemi per farlo funzionare, quindi ho seguito le istruzioni fornite in questo link e ora funziona bene http://glmmadmb.r-forge.r-project.org/

Spero che sia di aiuto!

— Rocioc
fonte

Si noti che questo pacchetto consente solo l'adattamento di un termine costante per la parte a zero inflazione del modello.

— user2390246,

1

Il pacchetto pscl prevede un modello Poisson a gonfiaggio zero. Non penso che possa fare un modello binomiale negativo, ma potrebbe essere un punto di partenza. L'articolo JSS collegato tratta anche i pacchetti correlati, che possono portare a ciò che stai cercando.

— Charlie
fonte

1

Il pacchetto pscl (ora) consente modelli binomiali negativi usandozeroinfl(..., dist = "negbin", ...)

— user2390246

0

A seconda di cosa stai cercando di fare, potresti voler guardare il pacchetto aster . I modelli di aster consentono l'analisi congiunta di più variabili con diverse distribuzioni di probabilità e recentemente sono state aggiornate per consentire effetti casuali . Sono stati progettati per l'analisi della storia della vita e funzioneranno in situazioni in cui è possibile dividere la risposta in parti distinte con diverse distribuzioni, (ad esempio sopravvivenza = Bernoulli, riproduzione = Poisson). Sono in grado di gestire "l'inflazione zero" modellando la maggior parte degli zero come bernoulli e il resto della risposta come binomiale negativo.

Troverai molta documentazione qui:

http://www.stat.umn.edu/geyer/aster/

— JohnSG
fonte