This is a small function Venables and Ripley provide in their MASS book. You don't need it anymore because vcov() has a method for the glm class. However, it is useful to see how to extract bits from a fitted model object.vcov.glm #return the variance-covariance matrix of a glm object #from p. 188 in Venables and Ripley....