How can I knit coeftest() regression outputs?

Question

I first created a probit model like so:

fit1 <-  glm(Y ~ X1, data = dat1, family = binomial (link="probit"))

and used huxtable t