Utilize a função coef
:
dataset <- mtcars
modtest <- lm(mpg ~ cyl + hp + wt + vs, dataset)
coef(summary(modtest))
#> Estimate Std. Error t value Pr(>|t|)
#> (Intercept) 38.48404581 3.33740414 11.53113145 6.143114e-12
#> cyl -0.90501219 0.67894306 -1.33297213 1.936803e-01
#> hp -0.01785427 0.01224378 -1.45823142 1.563120e-01
#> wt -3.18355380 0.77368907 -4.11477159 3.264650e-04
#> vs 0.15457896 1.61535124 0.09569372 9.244706e-01
Entretanto, esse objeto não é um data frame:
is.data.frame(coef(summary(modtest)))
#> [1] FALSE
Portanto, é necessário convertê-lo para um data frame:
as.data.frame(coef(summary(modtest)))
#> Estimate Std. Error t value Pr(>|t|)
#> (Intercept) 38.48404581 3.33740414 11.53113145 6.143114e-12
#> cyl -0.90501219 0.67894306 -1.33297213 1.936803e-01
#> hp -0.01785427 0.01224378 -1.45823142 1.563120e-01
#> wt -3.18355380 0.77368907 -4.11477159 3.264650e-04
#> vs 0.15457896 1.61535124 0.09569372 9.244706e-01
Mas se o objetivo for formatar essa tabela para um relatório ou apresentação, sugiro utilizar a função knitr::kable
. Ela permite pegar esses resultados e transformá-los em uma tabela em markdown (como essa que colei a seguir), html ou LaTeX. O pacote officer
ajuda a criar tabelas similares para o Word.
library(knitr)
kable(coef(summary(modtest)))
|
Estimate |
Std. Error |
t value |
Pr(>|t|) |
(Intercept) |
38.4840458 |
3.3374041 |
11.5311314 |
0.0000000 |
cyl |
-0.9050122 |
0.6789431 |
-1.3329721 |
0.1936803 |
hp |
-0.0178543 |
0.0122438 |
-1.4582314 |
0.1563120 |
wt |
-3.1835538 |
0.7736891 |
-4.1147716 |
0.0003265 |
vs |
0.1545790 |
1.6153512 |
0.0956937 |
0.9244706 |
Created on 2023-05-09 with reprex v2.0.2