Блог

Главная - Вопросы по программированию - Модуль Пуассона GLM: формульное представление взаимодействия

Модуль Пуассона GLM: формульное представление взаимодействия

#r #glm #poisson

#r #glm #poisson

Вопрос:

У меня есть вопрос о Пуассоне GLM и представлении формулы:

Учитывая набор данных:

 p <- read.csv("https://raw.githubusercontent.com/Leprechault/PEN-533/master/bradysia-greenhouse.csv")
  1. Без учета взаимодействия:
 
    m1 <- glm(bradysia ~ area   mes, family="quasipoisson", data=p)
        
    summary(m1)

 
    #(Intercept)  4.36395    0.12925  33.765  < 2e-16 ***
    #areaCV      -0.19696    0.12425  -1.585    0.113    
    #areaMJC     -0.71543    0.08553  -8.364 3.11e-16 ***
    #mes         -0.08872    0.01970  -4.503 7.82e-06 ***

Окончательная формула: bradysia = exp(4.36395*CS-0.19696*CV-0.71543-0.08872*mes)

  1. Учитывая взаимодействие:
     m2 <- glm(bradysia ~ area*mes, family="quasipoisson", data=p)
    summary(m2)
      
    #(Intercept)  4.05682    0.15468  26.227  < 2e-16 ***
    #areaCV       0.15671    0.35219   0.445   0.6565    
    #areaMJC      0.54132    0.31215   1.734   0.0833 .  
    #mes         -0.03943    0.02346  -1.680   0.0933 .  
    #areaCV:mes  -0.05724    0.05579  -1.026   0.3052    
    #areaMJC:mes -0.22609    0.05576  -4.055 5.57e-05 **

Окончательная формула такова: bradysia = exp(?????) и любая помощь, пожалуйста?

Комментарии:

1. y ~ x1*x2 в R означает то же самое, что y ~ x1 x2 x1:x2 или математически, что Y = a b1*x1 b2*x2 b3*x1*x2.

2. Спасибо @DanY, но и проблема с незначимыми коэффициентами? Потому что очень сложно найти категориальный и количественный фактор, смешанный в GLM со всеми коэффициентами, имеющими значение.