如何解决在线性模型结果中保留因子标签
如何在线性模型中保留标签?如何在下面的系数表中将斜率标记为 Bar
或 xBar
?
dat <- data.frame(
y = c(1:3,10:12),x = gl(2,3,labels = c('Foo','Bar'))
)
coef(summary(mod <- lm(y ~ x,dat)))
# Estimate Std. Error t value Pr(>|t|)
# (Intercept) 6.5 0.4082483 15.92168 0.0000909619
# x1 -4.5 0.4082483 -11.02270 0.0003850677
实际用例更复杂,几十个线性模型有几十个变量。我确实意识到我可以编写自己的函数来检索标签并拟合模型,但我想知道我是否遗漏了一些更简单的东西。
解决方法
它对我来说已经是这样了:
dat <- data.frame(
y = c(1:3,10:12),x = gl(2,3,labels = c('Foo','Bar'))
)
coef(summary(mod <- lm(y ~ x,dat)))
# Estimate Std. Error t value Pr(>|t|)
# (Intercept) 2 0.5773503 3.464102 0.0257214207
# xBar 9 0.8164966 11.022704 0.0003850677
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。