Andrzej Andrzej
Andrzej Andrzej

Reputation: 79

How to add a reference level for interaction in gtsummary and sjPlot

I hope my question makes sense and I want to understand what is the reference level for interaction term/characteristic/predictors, etc.

Here is my code:
library(gtsummary)
library(sjPlot)


set.seed(1000)
my_data <- rbind(
  data.frame(time = "Pre", treatment = "Control", response = rnorm(100, mean=1)),
  data.frame(time = "Pre", treatment = "Treatment", response = rnorm(100, mean=2)),
  data.frame(time = "Post", treatment = "Control", response = rnorm(100, mean=1)),
  data.frame(time = "Post", treatment = "Treatment", response = rnorm(100, mean=2))
) %>% mutate(time = factor(time, levels = c("Pre", "Post"))) %>%
mutate(treatment = factor(treatment, levels = c("Control", "Treatment")))

model3 <- lm(response ~ time * treatment, data = my_data)


gtsummary::tbl_regression(model3,
  pvalue_fun = ~ style_pvalue(.x, digits = 2),
  estimate_fun = ~ style_number(.x, digits = 4)
) %>%
  add_global_p() %>%
  bold_p(t = 0.10) %>%
  bold_labels() %>%
  italicize_levels()


sjPlot::tab_model(model3, show.reflvl = TRUE)

Which gives me:

enter image description here

I want to add a reference level (if possible) for interaction term to be explicitly visible and would be grateful for a bit of explanation how to interpret interaction effect here, I mean if it was statistically significant ?

In R usually those reference levels are hidden/omitted, but it would be good to include them in final summary table.


My desired table would be something like this:

enter image description here

Upvotes: 0

Views: 1000

Answers (0)

Related Questions