时间戳

首页 > 解决方案 > 如何将对应的统计数据从线性回归模型导出到 R 中的表格？

r - 如何将对应的统计数据从线性回归模型导出到 R 中的表格？

问题描述

我正在拟合具有 1 个因变量和 45 个自变量的多元线性回归模型。我想将相应的效果、95% CI 和 p 值放在表中。

我的数据集如下：

a  b  x1 x2 x3 x4 .... x45
23 15 1 34 4  45       8
10 45 2 15 2  55       18

有没有更方便的方法来编写回归语法而不是编写lm(y ~ x1 + x2 + x3 + x4......+x45, data=df)因为我有 45 个不同名称的变量？

我试过了

mod1 <- lm(a ~ x1 + x2 + ....+x45, data=df)
mod2 <- lm(b ~ x1 + x2 + .... +x45, data=df)

但我不知道如何将结果放入表格中，如下所示：

model variable     effect       95%CI      p
 mod1    x1       177.93  79.16- 276.71   0.003
 mod1    x2       -75.13 -116.46 - -33.8  0.003 
 ...
 mod1    x45      118.61  53.09-184.13    0.005
 mod2    x1       79.53   36.94 - 122.13  0.004
 mod2    x2       201.93  60.48 - 343.38  0.01 
 ...
 mod2    x45      61.56   20.87 - 102.25   0.005

如果有人可以提供帮助，那就太好了。非常感谢！

标签： rstatisticsexportregression

解决方案

当我使用 lmer 函数进行回归时，此代码可以提取结果表。请参考代码！

RESULT_model1 <- summary(model1<-lmer(LOGTHSMON_SELNG_AMT~ .-BLK_CD-LOGPERTHSMON_SELNG_AMT
                                      - STDR_YM_CD - bub_dong-M201507-M201508-M201509
                                      -STOR_CO
                                      +(1|BLK_CD)
                                      ,data=dta_fin_WEST))


RESULT_model1_cof <- as.data.frame(tidy(RESULT_model1$coefficients))
RESULT_model1_cof$p.value <- 2 * (1 - pnorm(abs(RESULT_model1_cof$t.value)))
RESULT_model1_cof$VAR <- RESULT_model1_cof$.rownames
RESULT_model1_cof$STDERR <- RESULT_model1_cof$Std..Error

推荐阅读