R/get_predictions_glm.R

In javifar/ggeffects: Create Tidy Data Frames of Marginal Effects for 'ggplot' from Model Outputs

get_predictions_glm <- function(model, fitfram, ci.lvl, linv, value_adjustment, model_class, terms, vcov.fun, vcov.type, vcov.args, condition, interval, type, ...) {
  # does user want standard errors?
  se <- !is.null(ci.lvl) && !is.na(ci.lvl) && is.null(vcov.fun)

  if (type == "sim") {
    # compute ci, two-ways
    if (!is.null(ci.lvl) && !is.na(ci.lvl))
      ci <- (1 + ci.lvl) / 2
    else
      ci <- .975
    # simulate predictions
    .do_simulate(model, terms, ci, ...)
  } else {
    # for models from "robust"-pkg (glmRob) we need to
    # suppress warnings about fake models
    prdat <-
      suppressWarnings(stats::predict.glm(
        model,
        newdata = fitfram,
        type = "link",
        se.fit = se,
        ...
      ))

    # copy predictions
    .generic_prediction_data(model, fitfram, linv, prdat, se, ci.lvl, model_class, value_adjustment, terms, vcov.fun, vcov.type, vcov.args, condition, interval)
  }
}