I am looking for a solution to compute hypergeometric confidence interval for multivariable binary logistic regression models output using glm(). Easy to do in univariate analys