Может быть, вы можете уточнить, что вы подразумеваете под хи-квадрат и exp (B). Вы можете сделать следующее:
da <- MASS::Pima.tr
model <- glm(type ~ .,data=da,family = binomial(link = cloglog))
results <- data.frame(coefficients(summary(model)),check.names=FALSE)
# some random values
results$chisq = rchisq(nrow(results),1)
results$expB = exp(results$Estimate)
Или вы можете использовать приборку из метлы:
library(broom)
results = tidy(model)
results$expB = exp(results$Estimate)