How to plot a binomial GLM with a three way interaction?

Question

I have the following model:

all.fit <-  glm(Presence/Total~Season*ToD*Site, 
                family = binomial, weights = Total, 
                data = all.dt)

... which I would like to plot.

Looking here it seems to be fairly simple but when I try this I get an error:

plot_model(all.fit, type="pred", terms=c("Site","Season","ToD"))
Error in [.data.frame(get_data(x, verbose = verbose), , rn, drop = FALSE) : 
  undefined columns selected

Using code as cited by Gung:

Plot of raw values:

Output of model:


> summary(all.fit)
Call:
glm(formula = Presence/Total ~ Season * ToD * Site, family = binomial, 
    data = all.dt, weights = Total)
Deviance Residuals: 
    Min       1Q   Median       3Q      Max

-6.0379  -1.3245  -0.6221   0.3961   8.6824
Coefficients:
                                     Estimate Std. Error z value Pr(>|z|)

(Intercept)                          -4.66074    0.26852 -17.357  < 2e-16 ***
SeasonSpring                          0.60478    0.33133   1.825 0.067959 .

SeasonSummer                         -0.82682    0.48901  -1.691 0.090873 .

SeasonWinter                          1.26868    0.30149   4.208 2.58e-05 ***
ToDDay                                0.09865    0.31530   0.313 0.754381

ToDDusk                               2.48420    0.28183   8.814  < 2e-16 ***
ToDNight                              2.65294    0.27185   9.759  < 2e-16 ***
SiteKawau                             3.09495    0.27605  11.212  < 2e-16 ***
SiteNoises                            3.48048    0.27479  12.666  < 2e-16 ***
SiteTawharanui                        2.94694    0.27684  10.645  < 2e-16 ***
SiteTiritiri                          3.21127    0.27557  11.653  < 2e-16 ***
SeasonSpring:ToDDay                   0.90540    0.37560   2.411 0.015928 *

SeasonSummer:ToDDay                   2.16685    0.52151   4.155 3.25e-05 ***
SeasonWinter:ToDDay                  -0.12465    0.35416  -0.352 0.724863

SeasonSpring:ToDDusk                 -2.78750    0.40957  -6.806 1.00e-11 ***
SeasonSummer:ToDDusk                 -1.52114    0.55644  -2.734 0.006263 ** 
SeasonWinter:ToDDusk                 -0.90286    0.32121  -2.811 0.004941 ** 
SeasonSpring:ToDNight                -1.66162    0.34346  -4.838 1.31e-06 ***
SeasonSummer:ToDNight                -1.85826    0.51861  -3.583 0.000339 ***
SeasonWinter:ToDNight                -0.21758    0.30592  -0.711 0.476946

SeasonSpring:SiteKawau               -0.67217    0.34426  -1.953 0.050875 .

SeasonSummer:SiteKawau                1.24760    0.49690   2.511 0.012047 *

SeasonWinter:SiteKawau               -1.68237    0.31709  -5.306 1.12e-07 ***
SeasonSpring:SiteNoises              -0.79023    0.34223  -2.309 0.020940 *

SeasonSummer:SiteNoises               0.42570    0.49791   0.855 0.392560

SeasonWinter:SiteNoises              -0.73055    0.31139  -2.346 0.018972 *

SeasonSpring:SiteTawharanui          -0.25773    0.34405  -0.749 0.453794

SeasonSummer:SiteTawharanui           1.32456    0.49760   2.662 0.007770 ** 
SeasonWinter:SiteTawharanui          -0.51055    0.31368  -1.628 0.103604

SeasonSpring:SiteTiritiri            -0.48838    0.34274  -1.425 0.154178

SeasonSummer:SiteTiritiri             0.93694    0.49749   1.883 0.059656 .

SeasonWinter:SiteTiritiri            -1.27870    0.31406  -4.072 4.67e-05 ***
ToDDay:SiteKawau                      0.70933    0.32351   2.193 0.028335 *

ToDDusk:SiteKawau                    -2.35320    0.29549  -7.964 1.67e-15 ***
ToDNight:SiteKawau                   -5.33351    0.29943 -17.812  < 2e-16 ***
ToDDay:SiteNoises                     0.27852    0.32258   0.863 0.387906

ToDDusk:SiteNoises                   -2.53814    0.29389  -8.636  < 2e-16 ***
ToDNight:SiteNoises                  -4.04432    0.28257 -14.313  < 2e-16 ***
ToDDay:SiteTawharanui                -0.16455    0.32550  -0.506 0.613193

ToDDusk:SiteTawharanui               -2.76975    0.29927  -9.255  < 2e-16 ***
ToDNight:SiteTawharanui              -2.62701    0.28227  -9.307  < 2e-16 ***
ToDDay:SiteTiritiri                   0.77569    0.32307   2.401 0.016350 *

ToDDusk:SiteTiritiri                 -2.27149    0.29441  -7.715 1.21e-14 ***
ToDNight:SiteTiritiri                -4.15841    0.28495 -14.593  < 2e-16 ***
SeasonSpring:ToDDay:SiteKawau        -0.19831    0.38943  -0.509 0.610596

SeasonSummer:ToDDay:SiteKawau        -1.61316    0.53079  -3.039 0.002372 ** 
SeasonWinter:ToDDay:SiteKawau         0.57193    0.37047   1.544 0.122635

SeasonSpring:ToDDusk:SiteKawau        2.48229    0.43071   5.763 8.25e-09 ***
SeasonSummer:ToDDusk:SiteKawau        1.37786    0.56992   2.418 0.015622 *

SeasonWinter:ToDDusk:SiteKawau        0.81702    0.34933   2.339 0.019347 *

SeasonSpring:ToDNight:SiteKawau       1.41405    0.40613   3.482 0.000498 ***
SeasonSummer:ToDNight:SiteKawau       2.06227    0.54597   3.777 0.000159 ***
SeasonWinter:ToDNight:SiteKawau       0.08689    0.36989   0.235 0.814279

SeasonSpring:ToDDay:SiteNoises       -0.28007    0.38782  -0.722 0.470188

SeasonSummer:ToDDay:SiteNoises       -1.48112    0.53201  -2.784 0.005369 ** 
SeasonWinter:ToDDay:SiteNoises       -0.06044    0.36566  -0.165 0.868708

SeasonSpring:ToDDusk:SiteNoises       1.93847    0.43143   4.493 7.02e-06 ***
SeasonSummer:ToDDusk:SiteNoises       0.63232    0.57685   1.096 0.273006

SeasonWinter:ToDDusk:SiteNoises       0.83278    0.34003   2.449 0.014321 *

SeasonSpring:ToDNight:SiteNoises      0.69621    0.37365   1.863 0.062422 .

SeasonSummer:ToDNight:SiteNoises      0.06911    0.55930   0.124 0.901655

SeasonWinter:ToDNight:SiteNoises      0.48690    0.32194   1.512 0.130428

SeasonSpring:ToDDay:SiteTawharanui   -0.23150    0.39050  -0.593 0.553287

SeasonSummer:ToDDay:SiteTawharanui   -1.96624    0.53250  -3.692 0.000222 ***
SeasonWinter:ToDDay:SiteTawharanui    0.43487    0.36870   1.179 0.238216

SeasonSpring:ToDDusk:SiteTawharanui   2.04614    0.43633   4.689 2.74e-06 ***
SeasonSummer:ToDDusk:SiteTawharanui   0.59597    0.57680   1.033 0.301491

SeasonWinter:ToDDusk:SiteTawharanui   0.73280    0.34653   2.115 0.034456 *

SeasonSpring:ToDNight:SiteTawharanui  1.02467    0.36109   2.838 0.004544 ** 
SeasonSummer:ToDNight:SiteTawharanui  1.19957    0.52990   2.264 0.023589 *

SeasonWinter:ToDNight:SiteTawharanui -0.16571    0.32155  -0.515 0.606306

SeasonSpring:ToDDay:SiteTiritiri     -0.49278    0.38807  -1.270 0.204143

SeasonSummer:ToDDay:SiteTiritiri     -1.63113    0.53142  -3.069 0.002145 ** 
SeasonWinter:ToDDay:SiteTiritiri      0.37405    0.36773   1.017 0.309069

SeasonSpring:ToDDusk:SiteTiritiri     1.81283    0.43042   4.212 2.53e-05 ***
SeasonSummer:ToDDusk:SiteTiritiri     0.83455    0.57225   1.458 0.144744

SeasonWinter:ToDDusk:SiteTiritiri     0.50070    0.34463   1.453 0.146263

SeasonSpring:ToDNight:SiteTiritiri    1.49159    0.36761   4.057 4.96e-05 ***
SeasonSummer:ToDNight:SiteTiritiri    1.30689    0.53796   2.429 0.015127 *

SeasonWinter:ToDNight:SiteTiritiri    0.79478    0.32728   2.428 0.015164 *

Signif. codes:  0 ‘*’ 0.001 ‘’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 104046  on 22997  degrees of freedom

Residual deviance:  62740  on 22918  degrees of freedom
AIC: 86778
Number of Fisher Scoring iterations: 6

Pure R programmng questions is off-topic here, can you reformulate the Q so it is more statistical? — kjetil b halvorsen, Jun 25 '21 at 21:34
As a suggestion, I absolutely see a statistics question hiding in there about how to plot the four dimensions. — Dave, Jun 25 '21 at 21:35
I've had some good luck with plotting partials with the MRSea package. The figures are a bit basic but you can dive into the code to get the output and tweak that to look a bit nicer. https://www.creem.st-andrews.ac.uk/software/ or https://github.com/lindesaysh/MRSea/tree/stable — K.J. Palmer, May 10 '22 at 17:49

score 1 · Answer 1 · answered Jun 25 '21 at 22:00

1

This is not much different than your prior question. Because you have all categorical variables, and because they are nested within a three-way interaction, you simply have a set of proportions. Get the predicted proportions from your model and make some barplots. You can see such a plot in my answer to: Test logistic regression model using residual deviance and degrees of freedom.

answered Jun 25 '21 at 22:00

gung - Reinstate Monica

145,122

Thank you! I think I am just having difficulty understanding as there are so many possible combinations in my data, and they don't appear to match up with the output of my model, but I know that it is my understanding at fault. I am trying to plot the data so that I can understand the model. I followed your code as cited above to get the plot above, hopefully you can see what I mean. Also, I could not get the binom.confint function to work, is that specific to a package? It seems to be base R? – Louise Wilson Jun 27 '21 at 22:51
1

I'm not sure what you've got. You will have to adapt the code. Based on your previous question, you only have 36 combinations. I would pick comparisons (sets of combinations) that are of theoretical interest, & plot only those. Regarding the function, notice that I use library(binom). – gung - Reinstate Monica Jun 28 '21 at 01:00

How to plot a binomial GLM with a three way interaction?

1 Answers1