2つのカテゴリ変数間の相互作用R

Question

2つのカテゴリ変数と年齢グループと性別を含む回帰モデルを扱っています。

2つのカテゴリ変数の間にインタラクション項を含める必要がありますが、結果のモデルには、すべての年齢層の女性とのやりとりの影響のみが表示されます。

「男性」を「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 

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()