lavaan
lavaan.Rmd
lavaan()
Results
library(lavaan)
#> This is lavaan 0.6-14
#> lavaan is FREE software! Please report any bugs.
mod <-
"
visual =~ x1 + x2 + x3
textual =~ x4 + x5 + x6
speed =~ x7 + x8 + x9
textual ~ c(a1, a2) * visual
speed ~ c(b1, b2) * textual
"
set.seed(870452)
id <- sample(nrow(HolzingerSwineford1939), 100)
fit <- sem(mod, data = HolzingerSwineford1939[id, ], group = "school")
est <- parameterEstimates(fit)
std <- standardizedSolution(fit)
est
#> lhs op rhs block group label est se z pvalue ci.lower
#> 1 visual =~ x1 1 1 1.000 0.000 NA NA 1.000
#> 2 visual =~ x2 1 1 0.267 0.213 1.251 0.211 -0.151
#> 3 visual =~ x3 1 1 0.720 0.281 2.565 0.010 0.170
#> 4 textual =~ x4 1 1 1.000 0.000 NA NA 1.000
#> 5 textual =~ x5 1 1 0.993 0.152 6.511 0.000 0.694
#> 6 textual =~ x6 1 1 0.965 0.130 7.440 0.000 0.711
#> 7 speed =~ x7 1 1 1.000 0.000 NA NA 1.000
#> 8 speed =~ x8 1 1 0.673 0.259 2.600 0.009 0.166
#> 9 speed =~ x9 1 1 0.450 0.187 2.410 0.016 0.084
#> 10 textual ~ visual 1 1 a1 0.682 0.269 2.538 0.011 0.155
#> 11 speed ~ textual 1 1 b1 0.180 0.153 1.180 0.238 -0.119
#> 12 x1 ~~ x1 1 1 0.211 0.233 0.907 0.364 -0.245
#> 13 x2 ~~ x2 1 1 1.220 0.233 5.235 0.000 0.763
#> 14 x3 ~~ x3 1 1 0.907 0.210 4.321 0.000 0.496
#> 15 x4 ~~ x4 1 1 0.378 0.109 3.452 0.001 0.163
#> 16 x5 ~~ x5 1 1 0.603 0.142 4.238 0.000 0.324
#> 17 x6 ~~ x6 1 1 0.219 0.087 2.534 0.011 0.050
#> 18 x7 ~~ x7 1 1 0.296 0.301 0.984 0.325 -0.294
#> 19 x8 ~~ x8 1 1 0.647 0.182 3.551 0.000 0.290
#> 20 x9 ~~ x9 1 1 0.605 0.130 4.659 0.000 0.350
#> 21 visual ~~ visual 1 1 0.683 0.282 2.420 0.016 0.130
#> 22 textual ~~ textual 1 1 0.603 0.201 3.008 0.003 0.210
#> 23 speed ~~ speed 1 1 0.819 0.356 2.304 0.021 0.122
#> 24 x1 ~1 1 1 5.083 0.126 40.213 0.000 4.836
#> 25 x2 ~1 1 1 5.746 0.151 38.168 0.000 5.450
#> 26 x3 ~1 1 1 2.556 0.150 17.027 0.000 2.262
#> 27 x4 ~1 1 1 2.917 0.152 19.153 0.000 2.618
#> 28 x5 ~1 1 1 4.018 0.164 24.461 0.000 3.696
#> 29 x6 ~1 1 1 2.038 0.139 14.693 0.000 1.766
#> 30 x7 ~1 1 1 4.545 0.143 31.776 0.000 4.265
#> 31 x8 ~1 1 1 5.662 0.136 41.712 0.000 5.396
#> 32 x9 ~1 1 1 5.500 0.118 46.713 0.000 5.270
#> 33 visual ~1 1 1 0.000 0.000 NA NA 0.000
#> 34 textual ~1 1 1 0.000 0.000 NA NA 0.000
#> 35 speed ~1 1 1 0.000 0.000 NA NA 0.000
#> 36 visual =~ x1 2 2 1.000 0.000 NA NA 1.000
#> 37 visual =~ x2 2 2 0.519 0.257 2.024 0.043 0.016
#> 38 visual =~ x3 2 2 0.969 0.321 3.019 0.003 0.340
#> 39 textual =~ x4 2 2 1.000 0.000 NA NA 1.000
#> 40 textual =~ x5 2 2 1.152 0.181 6.371 0.000 0.798
#> 41 textual =~ x6 2 2 1.002 0.160 6.264 0.000 0.689
#> 42 speed =~ x7 2 2 1.000 0.000 NA NA 1.000
#> 43 speed =~ x8 2 2 1.112 0.200 5.552 0.000 0.720
#> 44 speed =~ x9 2 2 0.942 0.185 5.084 0.000 0.579
#> 45 textual ~ visual 2 2 a2 0.665 0.236 2.814 0.005 0.202
#> 46 speed ~ textual 2 2 b2 0.376 0.162 2.318 0.020 0.058
#> 47 x1 ~~ x1 2 2 1.103 0.321 3.437 0.001 0.474
#> 48 x2 ~~ x2 2 2 1.200 0.272 4.409 0.000 0.667
#> 49 x3 ~~ x3 2 2 0.510 0.225 2.260 0.024 0.068
#> 50 x4 ~~ x4 2 2 0.461 0.121 3.824 0.000 0.225
#> 51 x5 ~~ x5 2 2 0.253 0.101 2.498 0.012 0.055
#> 52 x6 ~~ x6 2 2 0.235 0.082 2.850 0.004 0.073
#> 53 x7 ~~ x7 2 2 0.371 0.117 3.175 0.001 0.142
#> 54 x8 ~~ x8 2 2 0.238 0.117 2.038 0.042 0.009
#> 55 x9 ~~ x9 2 2 0.454 0.124 3.658 0.000 0.211
#> 56 visual ~~ visual 2 2 0.747 0.385 1.943 0.052 -0.007
#> 57 textual ~~ textual 2 2 0.423 0.175 2.417 0.016 0.080
#> 58 speed ~~ speed 2 2 0.555 0.191 2.908 0.004 0.181
#> 59 x1 ~1 2 2 5.045 0.205 24.605 0.000 4.644
#> 60 x2 ~1 2 2 6.608 0.178 37.023 0.000 6.258
#> 61 x3 ~1 2 2 2.159 0.166 13.010 0.000 1.834
#> 62 x4 ~1 2 2 3.515 0.166 21.159 0.000 3.190
#> 63 x5 ~1 2 2 4.812 0.169 28.516 0.000 4.482
#> 64 x6 ~1 2 2 2.669 0.150 17.779 0.000 2.375
#> 65 x7 ~1 2 2 3.925 0.153 25.615 0.000 3.625
#> 66 x8 ~1 2 2 5.507 0.155 35.545 0.000 5.203
#> 67 x9 ~1 2 2 5.422 0.154 35.238 0.000 5.121
#> 68 visual ~1 2 2 0.000 0.000 NA NA 0.000
#> 69 textual ~1 2 2 0.000 0.000 NA NA 0.000
#> 70 speed ~1 2 2 0.000 0.000 NA NA 0.000
#> ci.upper
#> 1 1.000
#> 2 0.685
#> 3 1.271
#> 4 1.000
#> 5 1.292
#> 6 1.220
#> 7 1.000
#> 8 1.181
#> 9 0.816
#> 10 1.208
#> 11 0.479
#> 12 0.668
#> 13 1.677
#> 14 1.319
#> 15 0.592
#> 16 0.882
#> 17 0.389
#> 18 0.887
#> 19 1.004
#> 20 0.859
#> 21 1.237
#> 22 0.997
#> 23 1.517
#> 24 5.331
#> 25 6.041
#> 26 2.850
#> 27 3.215
#> 28 4.340
#> 29 2.310
#> 30 4.825
#> 31 5.929
#> 32 5.731
#> 33 0.000
#> 34 0.000
#> 35 0.000
#> 36 1.000
#> 37 1.022
#> 38 1.599
#> 39 1.000
#> 40 1.507
#> 41 1.316
#> 42 1.000
#> 43 1.505
#> 44 1.306
#> 45 1.128
#> 46 0.694
#> 47 1.732
#> 48 1.734
#> 49 0.951
#> 50 0.698
#> 51 0.452
#> 52 0.396
#> 53 0.600
#> 54 0.466
#> 55 0.697
#> 56 1.501
#> 57 0.766
#> 58 0.930
#> 59 5.447
#> 60 6.958
#> 61 2.484
#> 62 3.841
#> 63 5.143
#> 64 2.963
#> 65 4.225
#> 66 5.810
#> 67 5.724
#> 68 0.000
#> 69 0.000
#> 70 0.000
std
#> lhs op rhs group label est.std se z pvalue ci.lower ci.upper
#> 1 visual =~ x1 1 0.874 0.150 5.820 0.000 0.580 1.168
#> 2 visual =~ x2 1 0.196 0.146 1.345 0.179 -0.089 0.481
#> 3 visual =~ x3 1 0.530 0.132 4.015 0.000 0.271 0.789
#> 4 textual =~ x4 1 0.842 0.054 15.514 0.000 0.736 0.949
#> 5 textual =~ x5 1 0.775 0.065 12.009 0.000 0.649 0.902
#> 6 textual =~ x6 1 0.892 0.048 18.612 0.000 0.798 0.986
#> 7 speed =~ x7 1 0.861 0.154 5.593 0.000 0.559 1.163
#> 8 speed =~ x8 1 0.611 0.136 4.494 0.000 0.345 0.877
#> 9 speed =~ x9 1 0.471 0.134 3.516 0.000 0.208 0.733
#> 10 textual ~ visual 1 a1 0.587 0.139 4.234 0.000 0.315 0.859
#> 11 speed ~ textual 1 b1 0.188 0.156 1.204 0.229 -0.118 0.493
#> 12 x1 ~~ x1 1 0.236 0.262 0.900 0.368 -0.278 0.751
#> 13 x2 ~~ x2 1 0.962 0.057 16.874 0.000 0.850 1.073
#> 14 x3 ~~ x3 1 0.719 0.140 5.136 0.000 0.445 0.993
#> 15 x4 ~~ x4 1 0.291 0.091 3.181 0.001 0.112 0.470
#> 16 x5 ~~ x5 1 0.399 0.100 3.987 0.000 0.203 0.595
#> 17 x6 ~~ x6 1 0.203 0.086 2.378 0.017 0.036 0.371
#> 18 x7 ~~ x7 1 0.259 0.265 0.976 0.329 -0.261 0.778
#> 19 x8 ~~ x8 1 0.627 0.166 3.772 0.000 0.301 0.952
#> 20 x9 ~~ x9 1 0.779 0.126 6.183 0.000 0.532 1.025
#> 21 visual ~~ visual 1 1.000 0.000 NA NA 1.000 1.000
#> 22 textual ~~ textual 1 0.655 0.163 4.023 0.000 0.336 0.974
#> 23 speed ~~ speed 1 0.965 0.058 16.512 0.000 0.850 1.079
#> 24 x1 ~1 1 5.374 0.525 10.235 0.000 4.345 6.403
#> 25 x2 ~1 1 5.100 0.500 10.198 0.000 4.120 6.081
#> 26 x3 ~1 1 2.275 0.253 8.988 0.000 1.779 2.772
#> 27 x4 ~1 1 2.559 0.276 9.263 0.000 2.018 3.101
#> 28 x5 ~1 1 3.269 0.337 9.713 0.000 2.609 3.928
#> 29 x6 ~1 1 1.963 0.229 8.587 0.000 1.515 2.412
#> 30 x7 ~1 1 4.246 0.423 10.041 0.000 3.417 5.075
#> 31 x8 ~1 1 5.574 0.543 10.258 0.000 4.509 6.639
#> 32 x9 ~1 1 6.242 0.605 10.321 0.000 5.057 7.428
#> 33 visual ~1 1 0.000 0.000 NA NA 0.000 0.000
#> 34 textual ~1 1 0.000 0.000 NA NA 0.000 0.000
#> 35 speed ~1 1 0.000 0.000 NA NA 0.000 0.000
#> 36 visual =~ x1 2 0.636 0.129 4.919 0.000 0.382 0.889
#> 37 visual =~ x2 2 0.379 0.156 2.433 0.015 0.074 0.685
#> 38 visual =~ x3 2 0.761 0.124 6.120 0.000 0.518 1.005
#> 39 textual =~ x4 2 0.788 0.068 11.599 0.000 0.654 0.921
#> 40 textual =~ x5 2 0.893 0.049 18.254 0.000 0.797 0.989
#> 41 textual =~ x6 2 0.874 0.052 16.796 0.000 0.772 0.975
#> 42 speed =~ x7 2 0.800 0.075 10.652 0.000 0.653 0.948
#> 43 speed =~ x8 2 0.880 0.066 13.361 0.000 0.751 1.009
#> 44 speed =~ x9 2 0.751 0.082 9.119 0.000 0.590 0.912
#> 45 textual ~ visual 2 a2 0.662 0.135 4.918 0.000 0.398 0.926
#> 46 speed ~ textual 2 b2 0.401 0.147 2.735 0.006 0.114 0.688
#> 47 x1 ~~ x1 2 0.596 0.164 3.629 0.000 0.274 0.918
#> 48 x2 ~~ x2 2 0.856 0.118 7.244 0.000 0.625 1.088
#> 49 x3 ~~ x3 2 0.420 0.189 2.220 0.026 0.049 0.792
#> 50 x4 ~~ x4 2 0.380 0.107 3.550 0.000 0.170 0.589
#> 51 x5 ~~ x5 2 0.202 0.087 2.311 0.021 0.031 0.373
#> 52 x6 ~~ x6 2 0.237 0.091 2.608 0.009 0.059 0.415
#> 53 x7 ~~ x7 2 0.359 0.120 2.988 0.003 0.124 0.595
#> 54 x8 ~~ x8 2 0.225 0.116 1.940 0.052 -0.002 0.452
#> 55 x9 ~~ x9 2 0.436 0.124 3.524 0.000 0.193 0.678
#> 56 visual ~~ visual 2 1.000 0.000 NA NA 1.000 1.000
#> 57 textual ~~ textual 2 0.562 0.178 3.150 0.002 0.212 0.911
#> 58 speed ~~ speed 2 0.839 0.118 7.134 0.000 0.609 1.070
#> 59 x1 ~1 2 3.709 0.423 8.765 0.000 2.880 4.539
#> 60 x2 ~1 2 5.581 0.614 9.093 0.000 4.378 6.784
#> 61 x3 ~1 2 1.961 0.258 7.609 0.000 1.456 2.466
#> 62 x4 ~1 2 3.190 0.372 8.576 0.000 2.461 3.919
#> 63 x5 ~1 2 4.299 0.482 8.911 0.000 3.353 5.244
#> 64 x6 ~1 2 2.680 0.323 8.297 0.000 2.047 3.314
#> 65 x7 ~1 2 3.862 0.438 8.809 0.000 3.002 4.721
#> 66 x8 ~1 2 5.359 0.591 9.070 0.000 4.201 6.517
#> 67 x9 ~1 2 5.312 0.586 9.065 0.000 4.164 6.461
#> 68 visual ~1 2 0.000 0.000 NA NA 0.000 0.000
#> 69 textual ~1 2 0.000 0.000 NA NA 0.000 0.000
#> 70 speed ~1 2 0.000 0.000 NA NA 0.000 0.000
The p-Values of Coefficients
text_coef_p(fit, "visual=~x2")
#> [1] "p = 0.211"
text_coef_p(fit, "visual=~x2", standardized = TRUE)
#> [1] "p = 0.179"
text_coef_p(fit, "visual=~x2", 4)
#> [1] "p = 0.2111"
text_coef_p(fit, "textual=~x5", 4)
#> [1] "p < 0.0001"
text_coef_p(fit, "visual=~x2", digits = 5, standardized = TRUE)
#> [1] "p = 0.17859"
Confidence Intervals
text_confint(fit, "visual=~x2")
#> [1] "[-0.151, 0.685]"
text_confint(fit, "visual=~x2", digits = 5)
#> [1] "[-0.15135, 0.68492]"
text_confint(fit, "visual=~x2", digits = 5, brackets = c("(", ")"))
#> [1] "(-0.15135, 0.68492)"