R中两个分类变量之间的交互作用

Question

我们正在处理一个包含两个分类变量（年龄组和性别）的回归模型。

我们想在这两个分类变量之间包括一个交互项，但是得到的模型只显示了女性与所有年龄组之间的交互作用效应。

我们该如何调整代码，以使其保留男性“26-30”岁作为参考水平，并显示所有其他组的效应输出？

count_med_op3 <- glm(Count_OP ~ Gender * age_group + otherfactors,
                     data = med, family = 'poisson')

需要的结果：

GenderMale:age_group"0-1" 
GenderMale:age_group"2-6"
GenderMale:age_group"7-18"
GenderMale:age_group"19-25"
GenderMale:age_group"31-36"
Genderfemale:age_group"0-1"
Genderfemale:age_group"2-6"
Genderfemale:age_group"7-18"
Genderfemale:age_group"19-25"
Genderfemale:age_group"26-30"
other factors

- Ali Turab Lotia

Answer 1

使用relevel函数：

# simulate some data
df_foo = data_frame(
  age = as.factor(sample(seq(10, 90, 10), 100, replace = TRUE)),
  y = rnorm(100),
  gender = as.factor(sample(c("Male", "Female"), 100, replace = TRUE))
)

# female as omitted level
df_foo %>% 
  lm(y ~ age*gender, data = .) %>% 
  summary()

# male as omitted level
df_foo %>% 
  lm(y ~ age*relevel(gender, ref = "Male"), data = .) %>% 
  summary()

- tchakravarty