Как установить нулевое значение GLM-перехвата положительным в R?

#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:

Вам не нужно будет устанавливать значение перехвата равным нулю для части модели с «нулевой инфляцией».

Модель с нулевой инфляцией состоит из двух частей:

  • часть подсчета (верхний бит вывода), которая не должна иметь отрицательных перехватов.

  • часть с нулевой инфляцией, которая моделирует, является ли наблюдение нулевым или нет. Это оценивается с помощью логит- или пробит-модели, которая, в свою очередь, может иметь отрицательные перехваты.