R/parse-fit_glm.R

In kota7/outreg: Regression Table for Publication

make_coef_part.glm <- function(fit, modelname, robust = FALSE, ...)
{
  out <- data.frame(modelname = modelname,
                    variable = names(coefficients(fit)),
                    coef = coefficients(fit),
                    stringsAsFactors = FALSE)

  se <- rep(NA, nrow(out))
  vcov_func <- ifelse(robust, sandwich, vcov)
  se[!is.na(coefficients(fit))] <- sqrt(diag(vcov_func(fit)))
  out$se <- se

  out$zv <- out$coef / out$se
  out$pv <- pnorm(-abs(out$zv))*2

  # define t value, identical with z value
  # this may help when displaying LM and GLM together
  out$tv <- out$zv

  rownames(out) <- NULL
  out
}


make_stat_part.glm <- function(fit, modelname, ...)
{
  data.frame(modelname = modelname,
             variable = '',
             nobs = length(fit$fitted.values),
             aic = AIC(fit),
             stringsAsFactors = FALSE)
}