# NOT RUN {
mod <- lm(Sepal.Length ~ Sepal.Width, iris)
mod %>% model_compute_terms_contributions()
mod <- lm(hp ~ mpg + factor(cyl) + disp:hp, mtcars)
mod %>% model_compute_terms_contributions()
mod <- glm(
response ~ stage * grade + trt,
gtsummary::trial,
family = binomial,
contrasts = list(
stage = contr.sum,
grade = contr.treatment(3, 2),
trt = "contr.SAS"
)
)
mod %>% model_compute_terms_contributions()
mod <- glm(
response ~ stage * trt,
gtsummary::trial,
family = binomial,
contrasts = list(stage = contr.poly)
)
mod %>% model_compute_terms_contributions()
mod <- glm(
Survived ~ Class * Age + Sex,
data = Titanic %>% as.data.frame(),
weights = Freq, family = binomial
)
mod %>% model_compute_terms_contributions()
d <- dplyr::as_tibble(Titanic) %>%
dplyr::group_by(Class, Sex, Age) %>%
dplyr::summarise(
n_survived = sum(n * (Survived == "Yes")),
n_dead = sum(n * (Survived == "No"))
)
mod <- glm(cbind(n_survived, n_dead) ~ Class * Age + Sex, data = d, family = binomial)
mod %>% model_compute_terms_contributions()
# }
Run the code above in your browser using DataLab