首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >R中两个范畴变量的相互作用

R中两个范畴变量的相互作用
EN

Stack Overflow用户
提问于 2016-08-17 10:25:57
回答 1查看 3.7K关注 0票数 3

我们正在处理一个回归模型,其中包含两个分类变量-年龄组和性别。

我们希望在这两个范畴变量之间包含一个相互作用项,但结果模型只显示了女性与所有年龄组之间相互作用的效果。

我们如何调整代码,使其将“男性”年龄保持在“26-30岁”作为参考水平,并在其输出中显示所有其他群体的影响?

调整码

代码语言:javascript
运行
复制
count_med_op3 <- glm(Count_OP ~ Gender * age_group + otherfactors,
                     data = med, family = 'poisson')

想要的结果:

代码语言:javascript
运行
复制
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
EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2016-08-17 10:42:01

使用relevel

代码语言:javascript
运行
复制
# 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()
票数 3
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/38994034

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档