Fit negative binomial regression models across the multiverse

Source: R/linear_models.R

glm.nb_mverse fits MASS::glm.nb across the multiverse according to model specifications provided by formula_branch. At least one formula_branch must have been added.

glm.nb_mverse(.mverse)

Arguments

.mverse

a mverse object.

Value

A mverse object with glm.nb fitted.

See also

Other model fitting functions: glm_mverse(), lm_mverse()

Examples

# \donttest{
# Displaying the multiverse table with \code{glm.nb} models fitted.
hurricane_outliers <- filter_branch(
  !Name %in% c("Katrina", "Audrey", "Andrew"),
  TRUE # include all
)
model_specifications <- formula_branch(alldeaths ~ MasFem)
mv <- create_multiverse(hurricane) %>%
  add_filter_branch(hurricane_outliers) %>%
  add_formula_branch(model_specifications) %>%
  glm.nb_mverse()
summary(mv)
#> # A tibble: 4 × 12
#>   universe hurricane_outliers_branch model_specifications_branch term   estimate
#>   <fct>    <fct>                     <fct>                       <chr>     <dbl>
#> 1 1        hurricane_outliers_1      model_specifications_1      (Inte…   2.38  
#> 2 1        hurricane_outliers_1      model_specifications_1      MasFem   0.0855
#> 3 2        hurricane_outliers_2      model_specifications_1      (Inte…   2.23  
#> 4 2        hurricane_outliers_2      model_specifications_1      MasFem   0.201 
#> # ℹ 7 more variables: std.error <dbl>, statistic <dbl>, p.value <dbl>,
#> #   conf.low <dbl>, conf.high <dbl>, hurricane_outliers_branch_code <fct>,
#> #   model_specifications_branch_code <fct>
# }