Extracting model equation and other data from 'glm' function in R
In order to extract some data from the fitted glm model object, you need to figure out where that data resides (use documentation and str() for that). Some data might be available from the summary.glm object, while more detailed data is available from the glm object itself. For extracting model parameters, you can use coef() function or direct access to the structure.

UPDATE:

From Princeton's* introduction to R course's website, GLM section - see for details &amp; examples:


  The functions that can be used to extract results from the fit
  include

- 'residuals' or 'resid', for the deviance residuals
- 'fitted' or 'fitted.values', for the fitted values (estimated probabilities)
- 'predict', for the linear predictor (estimated logits)
- 'coef' or 'coefficients', for the coefficients, and
- 'deviance', for the deviance. 

  
  Some of these functions have optional arguments; for example, you can
  extract five different types of residuals, called "deviance",
  "pearson", "response" (response - fitted value), "working" (the
  working dependent variable in the IRLS algorithm - linear predictor),
  and "partial" (a matrix of working residuals formed by omitting each
  term in the model). You specify the one you want using the type
  argument, for example residuals(lrfit,type="pearson").


*) More accurately, this website is by Germán Rodríguez from Princeton University.