Regression example
# Load libraries
library(lavaan)
library(lavaanExtra)
# Define our regression terms
regression <- list(mpg = names(mtcars)[2:5],
disp = names(mtcars)[4:7])
# Write the model, and check it
mtcars.model <- write_lavaan(regression = regression)
cat(mtcars.model)
## ##################################################
## # [---------Regressions (Direct effects)---------]
##
## mpg ~ cyl + disp + hp + drat
## disp ~ hp + drat + wt + qsec
# Fit the model with `lavaan`
fit.reg <- sem(mtcars.model, data = mtcars)
# Get regression parameters only
lavaan_reg(fit.reg)
## Outcome Predictor B p
## 1 mpg cyl -0.242 0.244
## 2 mpg disp -0.287 0.123
## 3 mpg hp -0.264 0.128
## 4 mpg drat 0.192 0.159
## 5 disp hp 0.175 0.177
## 6 disp drat -0.177 0.034
## 7 disp wt 0.614 0.000
## 8 disp qsec -0.186 0.061
# We can get it prettier with the `rempsyc::nice_table` integration
lavaan_reg(fit.reg, nice_table = TRUE, highlight = TRUE)
Outcome |
Predictor |
β |
p |
---|---|---|---|
mpg |
cyl |
-0.24 |
.244 |
mpg |
disp |
-0.29 |
.123 |
mpg |
hp |
-0.26 |
.128 |
mpg |
drat |
0.19 |
.159 |
disp |
hp |
0.18 |
.177 |
disp |
drat |
-0.18 |
.034 |
disp |
wt |
0.61 |
< .001 |
disp |
qsec |
-0.19 |
.061 |