Hello, Is there a way to obtain the variance-covariance matrix of the estimated parameters from GLM?
my.glm<-glm(mat ~X,family = binomial, data =myDATA) out1<-predict(my.glm,se.fit = TRUE) std<-out1$se.fit se.fit is for getting the standard errors of the estimated parameters (\betas). Is there a way to get the variance-covariance matrix of the estimated parameters? Many thanks, Barbara ______________________________________________ R-help@r-project.org 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.