用Broom包求多元线性回归Anova表中固定值的P值
用Broom包求多元线性回归Anova表中固定值的P值
提问于 2021-10-03 13:47:43
在多元线性回归lm(FE_FCE2 ~ Trial + .x, data = DF_FCE3)中，有一个固定变量(试验)和多个x变量。我正在用试验作为固定效应来分析每个x变量对FE_FCE2的影响。然后，我使用boom软件包进行许多回归，并将结果绘制在一个表中。我已经得到了回归结果的结果。但是，不能将ANOVA表中的数据添加到带有map函数的Broom包中。
有可能吗？那么是怎样做的呢？
我使用以下公式从回归结果中获取数据：
DF_FCE3 %>%
select(-FE_FCE2, -Trial) %>% # exclude outcome, leave only predictors
map( ~lm(FE_FCE2 ~ Trial + .x, data = DF_FCE3)) %>%
map(summary) %>%
map_df(glance) %>%
round(3) -> rsme但是，我想从试验的方差分析表中获得P值(**4.26e-08 *****)。看看试验是否对x变量有显著影响。
**$x1
Analysis of Variance Table
**Response: FE_FCE2
Df Sum Sq Mean Sq F value Pr(>F)
Trial 3 0.84601 0.282002 15.0653 **4.26e-08 *****
.x 1 0.00716 0.007161 0.3826 0.5377
Residuals 95 1.77827 0.018719**
---**是否可以使用带有map函数的broom包来获得包含anova回归的所有p值的表？
回答 1
Stack Overflow用户
回答已采纳
发布于 2021-10-03 15:16:29
像这样(使用mpg)？这将返回一个数据帧，其中包含原始列，一行包含除结果和目标之外的p值(本例中为hwy和cyl，本例中为FE_FCE2和Trial )。
mpg %>%
select(-hwy, -cyl) %>% # exclude outcome, leave only predictors
map( ~lm(hwy ~ cyl + .x, data = mpg)) %>%
map(anova) %>%
map(broom::tidy) %>%
map_df(~.$p.value[1])页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/69425435
复制相关文章
相似问题
腾讯云开发者
