R/parse-fit_survreg.R

In outreg: Regression Table for Publication

make_coef_part.survreg <- function(fit, modelname, robust = FALSE, ...)
{
  nvar <- nrow(fit$var)
  out <- data.frame(modelname = rep(modelname, nvar),
                    stringsAsFactors = FALSE)

  coef <- coefficients(fit)
  cname <- names(coef)

  nvar0 <- length(coef)
  if (nvar > nvar0) {
    coef <- c(coef, log(fit$scale))
    if ((nvar-nvar0) == 1) {
      cname <- c(cname, 'Log(scale)')
    } else {
      cname <- c(cname, paste0('Log(scale) ', names(fit$scale)))
    }
  }
  out$variable <- cname
  out$coef <- coef

  se <- rep(NA, nrow(out))
  vcov_func <- ifelse(robust, sandwich, vcov)
  se[!is.na(coef)] <- 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

  # separate scale part
  if (nvar > nvar0) {
    list(out[1:nvar0, ], out[(nvar0+1):nvar, ])
  } else {
    list(out)
  }
}


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