[R-sig-ME] Marginal effect plot for an interaction effect from GLMM following a beta-distribution (glmmTMB)

[Ben Bolker]

   I always get confused about the different meanings of the term 
'marginal', but it's possible that the emmeans package will do what you 
want ...

On 10/19/20 4:21 PM, Jan anye Velimsky via R-sig-mixed-models wrote:
> 
> Dear R project mixedmodels users,
> 
> I am struggling to estimateand visualize the marginal effects (interaction effect) from a GLMM Model.
> 
> We have been estimatinga GLMM model following a beta-distribution with the glmmTMB-package. The model consistsof factors influencing referendum turnout (range between 0-1) in german municipalities.The primary units of investigation are city districts nested in referendumsnested in cities.
> 
> Here an example model:
> 
> model1 <- glmmTMB (ref_turnout ~ unemployment + contestation+ experience_with_ref + (experience_with_ref *unemployment)+ (1| town/ referendum), family=list(family="beta", link ='logit'), data = ml)
> 
>   
> 
> Results example model:
> 
>                                                                                                    Estimate      Std. Error    zvalue     Pr(>|z|)
> 
> (Intercept)                                                                                -0.583             0.131      -4.455     8.4e-06***
> 
> unemployment                                                                        -0.067              0.002      -30.397    < 2e-16 ***
> 
> contestation                                                                             0.008              0.003        2.398     0.0165 *
> 
> experience_with_ref                                                               0.012              0.052       0.228      0.8200
> 
> unemployment: experience_with_ref                                 -0.001             0.001      -1.725      0.0845 .
> 
>   
> 
> The model contains an interaction effect with unemploymentper district and experience with referendums in the respective city.
> 
> The aim is to visualize the interaction effects plotting the marginal effects.
> 
> With the "marginal_effects" command from the (margins-package) it was possible to estimate marginal effects,but no further summaries which are needed for plotting the model.  The "margins" command does not work for the glmmTMBmodel.  The idea is to have a simplemarginal effects plot for the interaction effect like this:
> 
> 
> Thanks a lot for your help!
> Jan Velimsky
> 
> -- Jan Velimsky, MAResearch associate
> Ludwig-Maximilians-Universität Munich
> Geschwister-Scholl-Institut of Political Science (GSI)
> Oettingenstrasse 67
> D-80538 Munich
> 
> 
> 
> 
> 
> 
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>