Help with mixed logistic model ( beta-binomial)

["Hunsicker, Lawrence" <lawrence-hunsicker@uiowa.edu>]

Greetings, all:

I have run a mixed logistic model to explore the variation in center practice in ordering a diagnostic test in patients initiating hemodialysis, corrected for baseline patient level covariates (fixed effects), using glmmPQL in the MASS library. 

The call:  glmmPQL(fixed = any~FIRST.SE, random = 1|PROVUSRD, …)

The model runs successfully and I have the output.  But for graphical purposes I’d like to be able to graph the estimated likelihood of testing by center (PROVUSRD).  Is there a neat way to get from the model object a list of the fitted values for each center, rather than for each patient? 

Many thanks in advance for any help that you can give me on this.

Larry Hunsicker