#r #&lm
#r #&lm
Вопрос:
Я пытаюсь запустить некоторые GLM с нулевым значением, используя функцию zeroinfl, но многие из них дают мне отрицательные перехваты, которые не имеют смысла для моих переменных. Есть ли способ установить значение перехвата выше нуля? Спасибо! Например:
zeroinfl(formula = crab.burrows.m2 ~ Al&ae, data = ve&data)
Pearson residuals:
Min 1Q Median 3Q Max
-1.0557 -1.0557 -0.7236 0.1067 11.8967
Count model coefficients (poisson with lo& link):
Estimate Std. Error z value Pr(&&t;|z|)
(Intercept) 2.40350 0.03991 60.216 <2e-16 ***
Al&ae 0.03964 0.02067 1.918 0.0551 .
Zero-inflation model coefficients (binomial with lo&it link):
Estimate Std. Error z value Pr(&&t;|z|)
(Intercept) -0.30108 0.20093 -1.498 0.1340
Al&ae 0.18556 0.09047 2.051 0.0403 *
---
Si&nif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Number of iterations in BFGS optimization: 1
Lo&-likelihood: -525.5 on 4 Df```
Ответ №1:
Вам не нужно будет устанавливать значение перехвата равным нулю для части модели с «нулевой инфляцией».
Модель с нулевой инфляцией состоит из двух частей:
-
часть подсчета (верхний бит вывода), которая не должна иметь отрицательных перехватов.
-
часть с нулевой инфляцией, которая моделирует, является ли наблюдение нулевым или нет. Это оценивается с помощью логит- или пробит-модели, которая, в свою очередь, может иметь отрицательные перехваты.