title: GLM in R: Normal Q-Q Plot with the 95% confidence interval
body: <p>I'd like a function or package to plot the Normal Q-Q Plot with the 95% confidence interval, but I don't find for GLM, only GAM models and for response variables in package <code>car</code>. In my example:</p>
  <pre><code>#Data set example
  p &lt;- read.csv(&quot;https://raw.githubusercontent.com/Leprechault/PEN-533/master/bradysia-greenhouse.csv&quot;)
  
  #Quasi Poisson GLM
  m1 &lt;- glm(bradysia ~ area + mes, family=&quot;quasipoisson&quot;, data=p)
  
  #Normal Q-Q Plot 
  plot(m1, which = 2)
  </code></pre>
  <p><a href="https://i.stack.imgur.com/mCfIs.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/mCfIs.png" alt="pl1" /></a></p>
  <pre><code>#Normal Q-Q Plot with the confidence interval
  </code></pre>
  <p><a href="https://i.stack.imgur.com/AlqNV.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/AlqNV.png" alt="pl2" /></a></p>
  <p>I'll like a plot like above for GLM models, it is possible? There are any function or package?</p>