如何提取线性模型的系数并将其存储在 R 中的变量中？答案

【问题标题】：How to extract the coefficients of a linear model and store in a variable in R?如何提取线性模型的系数并将其存储在 R 中的变量中？
【发布时间】：2021-05-09 07:25:32
【问题描述】：

我有一个数据框，我做了一个线性模型。我想提取系数并将每个系数存储到使用 R 的变量中。这是我的数据框

df <- mtcars
fit <- lm(mpg~., data = df)

这就是我提取一个系数的方法

beta_0 = fit$coefficients[1]

我想对模型中的所有系数自动执行此操作。我尝试使用循环但不工作。我知道不是正确的代码，但这就是我发现的

for (i in fit$coefficients(1:11)) {
  d["s{0}".format(x)] = variable1
}

【问题讨论】：

什么是d。我认为您可以只做coef(fit) 来返回所有系数，而不是在全局环境中创建多个对象
这看起来你可能正在混合 Python 和 R 语法？？
嗨。感谢您的回答。不，这一切都在 R 中。我尝试使用该函数将值存储在不同的变量中。我需要将每个系数放在不同的变量中才能运行蒙特卡罗模拟。这就是我想要找到的，但我没有得到任何答案，所以我想一次解决一个问题。 stackoverflow.com/questions/66053962/…
也许你应该删除另一个问题，如果它本质上是这个问题的副本？上一个问题只有 2 小时的时间，您不应该立即获得答案。
您的其他问题实际上更相关。我认为您实际上正在寻找的答案是您应该将新的随机变量变成 matrix 并使用矩阵乘法来计算预测值：X %*% coef(fit)。然后，您永远不需要将系数向量解压缩成一堆单独的变量（并弄乱您的工作区等）

标签： r linear-regression

【解决方案1】：

df <- mtcars
fit <- lm(mpg~., data = df)

beta_0 = fit$coefficients[1]

#base R approach
coef_base <- coef(fit)
coef_base
#> (Intercept)         cyl        disp          hp        drat          wt 
#> 12.30337416 -0.11144048  0.01333524 -0.02148212  0.78711097 -3.71530393 
#>        qsec          vs          am        gear        carb 
#>  0.82104075  0.31776281  2.52022689  0.65541302 -0.19941925


#tidyverse approach with the broom package
coef_tidy <- broom::tidy(fit)
coef_tidy
#> # A tibble: 11 x 5
#>    term        estimate std.error statistic p.value
#>    <chr>          <dbl>     <dbl>     <dbl>   <dbl>
#>  1 (Intercept)  12.3      18.7        0.657  0.518 
#>  2 cyl          -0.111     1.05      -0.107  0.916 
#>  3 disp          0.0133    0.0179     0.747  0.463 
#>  4 hp           -0.0215    0.0218    -0.987  0.335 
#>  5 drat          0.787     1.64       0.481  0.635 
#>  6 wt           -3.72      1.89      -1.96   0.0633
#>  7 qsec          0.821     0.731      1.12   0.274 
#>  8 vs            0.318     2.10       0.151  0.881 
#>  9 am            2.52      2.06       1.23   0.234 
#> 10 gear          0.655     1.49       0.439  0.665 
#> 11 carb         -0.199     0.829     -0.241  0.812

for (i in coef_base) {
  #do work on i
  print(i)
}
#> [1] 12.30337
#> [1] -0.1114405
#> [1] 0.01333524
#> [1] -0.02148212
#> [1] 0.787111
#> [1] -3.715304
#> [1] 0.8210407
#> [1] 0.3177628
#> [1] 2.520227
#> [1] 0.655413
#> [1] -0.1994193

【讨论】：

谢谢。如何将每个系数存储到不同的变量中？
不确定我是否理解您的问题。如果您是 R 初学者，我强烈建议您略读本书的前几章adv-r.hadley.nz