title: Confidence intervals of odds ratios of glmmPQL object (logistic regression)
body: <p>Computation of confidence intervals for odds ratios of a <code>glm</code> object is straight forward:</p>
  <pre><code>mydata &lt;- read.csv(&quot;http://www.ats.ucla.edu/stat/data/binary.csv&quot;)
  # set up logistic model
  mylogit &lt;- glm(admit ~ gre + gpa + rank, data = mydata, family = &quot;binomial&quot;)
  # calc odds ratios and confidence intervals
  exp(cbind(OR = coef(mylogit), confint(mylogit)))
  </code></pre>
  <p><a href="http://www.ats.ucla.edu/stat/r/dae/logit.htm" rel="nofollow noreferrer">Dataset and code credit</a></p>
  <p>However, the same approach does not work on <code>glmmPQL</code> objects (with random effects) returned by <code>MASS:glmmPQL()</code>.</p>
  <p>MWE from glmmPQL() help:</p>
  <pre><code>library(nlme) # will be loaded automatically if omitted
  mylogit &lt;- summary(glmmPQL(y ~ trt + I(week &gt; 2), random = ~ 1 | ID,
  +                 family = binomial, data = bacteria))
  exp(cbind(OR = coef(mylogit), confint(mylogit)))
  </code></pre>
  <p>Manually computing the confidence intervals (coefficient + 2* Std. Error) using the coefficients of the summary output also failed for me as I could not find them in the summary object (but they are there when printing the summary to console). Taking them out manually is no longterm approach...</p>
  <p>How does this work efficiently for <code>glmmPQL</code> objects?</p>