title: How to plot predicted data of a GLMM from the package glmmTMB?
body: <p>i have the following data and created a model with the package glmmTMB in R for plant diameters ~ plant density (number of plants) with a random plot effect:</p>
  <pre><code>d &lt;- data.frame (diameter = c(17,16,15,13,11, 19,17,15,11,11, 19,15,14,11,8),
                        plant_density = c(1000,2000,3000,4000,5000, 1000,2000,3000,4000,5000, 1000,2000,3000,4000,5000),
                        plot = c(1,1,1,1,1, 2,2,2,2,2, 3,3,3,3,3))
  
  glmm.model &lt;- glmmTMB(diameter ~  plant_density + (1|plot),
                          data = d,
                          na.action = na.omit,
                          family=&quot;gaussian&quot;,
                          ziformula = ~ 0)
  </code></pre>
  <p>My intention was to create a plot with predicted diameter data for different plant densities with an included random plot effect. So i tried to predict the data:</p>
  <pre><code>new.dat &lt;- data.frame(diameter= d$diameter,
                        plant_density = d$plant_density,
                        plot= d$plot) 
  
  new.dat$prediction &lt;- predict(glmm.model, new.data = new.dat, 
                                type = &quot;response&quot;, re.form = NA)
  </code></pre>
  <p>Unfortunately I get an output for every plot but wanted a generalized prediction for the diameter ~ plant density.</p>
  <p>My goal is to create a plot like <a href="https://i.stack.imgur.com/8wwnN.png" rel="nofollow noreferrer">here</a>, but with a regression model from glmmTMB which consider the random effect.</p>
  <p>Thanks for ur help!</p>