Hello Achim Zeileis, > We are trying to find out how to get the variance-covariance matrix of the > MLEs out of the glm function. Can anyone help?
Please try this : Suppose you have summary(glm(A ~ B,data=abc, family=Gamma)) then use summary(glm(A ~ B,data=abc, family=Gamma))$cov.unscaled I hope this will help. Cheers, Gaurav [[alternative HTML version deleted]] ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.