Модель Z Zeroinfl - PullRequest
       16

Модель Z Zeroinfl

1 голос
/ 13 июля 2011

Я выполняю нулевой раздуваемый отрицательный биномиальный GLM для некоторых данных о количестве насекомых в R. Моя проблема заключается в том, чтобы заставить R считывать данные о моем виде в виде одного столбца с накоплением, чтобы сохранить нулевую инфляцию. Если я подвожу итоги и импортирую его в R как одну строку под названием «Изобилие», я теряю нули, и модель не работает. Я уже пытался:

  1. сами складывают данные (есть 80 столбцов * 47 строк), поэтому с 3760 строками после суммирования вручную вы можете представить, насколько медленным становится R при использовании команды pscl zeroinfl() (требуется 20 минут на моем компьютере !, он все еще работал)

  2. Следующая проблема касается пространственной корреляции. Определенные пробоотборники отбираются с одного и того же носителя, чтобы нарушить независимость. Могу ли я просто указать средний фактор в модели?

1 Ответ

0 голосов
/ 01 сентября 2012

3760 строк занимают 20 минут с PSCL? Боже мой, у меня битва 30.000 рядов :), поэтому мой подсчет pscl не закончился ...

Однако затем я работал с GLMM, включая вложенные случайные эффекты (lme / gamm) и отрицательное биномиальное распределение, устанавливая тета на низкое значение, чтобы распределение обрабатывало нулевую инфляцию. Я думаю, что это зависит от степени нулей. В моем случае это было 44%, а остатки выглядели довольно хорошо.

...