Как я могу изобразить взаимодействие в R при использовании отрицательной биномиальной модели? - PullRequest
Новая
       21

Как я могу изобразить взаимодействие в R при использовании отрицательной биномиальной модели?

0 голосов
/ 10 июля 2020

Я новичок в R. Я использовал отрицательную биномиальную модель, чтобы проверить влияние 2 переменных (1 двоичная категориальная переменная и 1 непрерывная переменная) и переменной ответа подсчета. Я также добавил их взаимодействие в модель.

Поскольку результаты glm.nb минимальны, я хотел бы каким-то образом отобразить результат, особенно взаимодействие.

Я сделал это для запуска модели: 

Y<- cbind(N_Cooperations)
Model8 <- glm.nb(Y ~ Condition + NR + Condition*NR)
summary(Model8)

Call:
glm.nb(formula = Y ~ Condition + NR + Condition * NR, init.theta = 2.012332023, 
    link = log)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.2063  -0.9508  -0.1757   0.3389   2.5682  

Coefficients:
             Estimate Std. Error z value Pr(>|z|)   
(Intercept)    1.5379     0.6920   2.222  0.02626 * 
Condition     -2.9514     1.0876  -2.714  0.00665 **
NR            -0.1470     0.2065  -0.712  0.47654   
Condition:NR   0.7771     0.3170   2.451  0.01423 * 
---

Затем я попытался использовать plot(allEffects(Model8)) из пакета effects для построения графика взаимодействия, но получил это сообщение: 

> plot(allEffects(Model8))
Error in mod.matrix[, components] : subscript out of bounds
In addition: Warning messages:
1: In factor.cols & stranger.cols :
  longer object length is not a multiple of shorter object length
2: In (!factor.cols) & stranger.cols :
  longer object length is not a multiple of shorter object length

Что мне не хватает?

Опять же, я новичок в R. Заранее извините, если это звучит глупо.

...