Compute the variance covariance matrix (i.e. inverse of the information) for logistic regressions.
.vcov.logit(object, indiv, center)
a glm object corresponding to a logistic regression.
[logical] should the individual contribution be output instead of the total variance-covariance?
[logical] should the individual contribution be centered around the average?