### Bootstrapped mediation of latent constructs
### First, set the working directory where R will find the data
setwd("H:/data/User/Classes/multiv_13/Class11")
getwd()
## [1] "H:/data/User/Classes/multiv_13/Class11"
### Read the data
library(haven)
Wk12<- read_sav("class11_homework_2020.sav")
head(Wk12)
## # A tibble: 6 x 15
## Fixadur Fixanum Fixlocs Energy FAS FSS HRQOL BDI2 GAD7 SixMinWalk
## <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 50.1 49.6 51.7 48.5 48.0 50.5 40.7 46.4 41.3 50.8
## 2 37.1 38.5 46.0 58.4 56.2 49.3 55.0 56.4 53.4 57.9
## 3 45.4 45.2 46.5 48.7 47.3 43.1 46.9 55.5 54.2 44.2
## 4 49.3 49.2 48.2 42.9 42.1 49.0 50.2 52.1 51.2 54.7
## 5 47.1 41.7 46.3 47.0 48.9 49.1 49.9 50.7 46.9 50.2
## 6 58.8 53.3 54.0 44.7 54.1 42.9 53.3 62.7 57.9 50.4
## # ... with 5 more variables: TenMWalk <dbl>, BalConf <dbl>, Disclose <dbl>,
## # Hopeful <dbl>, TherComp <dbl>
names(Wk12)
## [1] "Fixadur" "Fixanum" "Fixlocs" "Energy" "FAS"
## [6] "FSS" "HRQOL" "BDI2" "GAD7" "SixMinWalk"
## [11] "TenMWalk" "BalConf" "Disclose" "Hopeful" "TherComp"
#Pre-SEM: Basic (CFA) Measurement model
library(lavaan)
## This is lavaan 0.6-7
## lavaan is BETA software! Please report any bugs.
#Confirmatory factor model with raw data
Wk10.mod<- ' Attention =~ Fixadur + Fixanum + Fixlocs
ZoomFatigue =~Energy + FAS + FSS
QOL =~ HRQOL + BDI2 + GAD7
Mobility =~ SixMinWalk + TenMWalk + BalConf
TherAlli =~ Disclose + Hopeful+ TherComp '
fit.Wk10 <- cfa(Wk10.mod, data=Wk12[1:15])
summary(fit.Wk10, fit.measures=TRUE, standardized=TRUE)
## lavaan 0.6-7 ended normally after 181 iterations
##
## Estimator ML
## Optimization method NLMINB
## Number of free parameters 40
##
## Number of observations 200
##
## Model Test User Model:
##
## Test statistic 76.026
## Degrees of freedom 80
## P-value (Chi-square) 0.605
##
## Model Test Baseline Model:
##
## Test statistic 2556.997
## Degrees of freedom 105
## P-value 0.000
##
## User Model versus Baseline Model:
##
## Comparative Fit Index (CFI) 1.000
## Tucker-Lewis Index (TLI) 1.002
##
## Loglikelihood and Information Criteria:
##
## Loglikelihood user model (H0) -9062.990
## Loglikelihood unrestricted model (H1) -9024.977
##
## Akaike (AIC) 18205.980
## Bayesian (BIC) 18337.912
## Sample-size adjusted Bayesian (BIC) 18211.188
##
## Root Mean Square Error of Approximation:
##
## RMSEA 0.000
## 90 Percent confidence interval - lower 0.000
## 90 Percent confidence interval - upper 0.035
## P-value RMSEA <= 0.05 0.997
##
## Standardized Root Mean Square Residual:
##
## SRMR 0.028
##
## Parameter Estimates:
##
## Standard errors Standard
## Information Expected
## Information saturated (h1) model Structured
##
## Latent Variables:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## Attention =~
## Fixadur 1.000 6.885 0.893
## Fixanum 1.044 0.055 18.837 0.000 7.189 0.932
## Fixlocs 0.880 0.056 15.733 0.000 6.062 0.830
## ZoomFatigue =~
## Energy 1.000 6.634 0.900
## FAS 1.054 0.059 17.932 0.000 6.993 0.909
## FSS 0.911 0.059 15.397 0.000 6.046 0.824
## QOL =~
## HRQOL 1.000 6.632 0.896
## BDI2 1.083 0.054 20.089 0.000 7.181 0.928
## GAD7 0.961 0.054 17.831 0.000 6.370 0.874
## Mobility =~
## SixMinWalk 1.000 6.911 0.889
## TenMWalk 1.031 0.052 19.755 0.000 7.127 0.929
## BalConf 0.976 0.054 18.026 0.000 6.747 0.885
## TherAlli =~
## Disclose 1.000 6.641 0.890
## Hopeful 1.040 0.055 19.033 0.000 6.909 0.921
## TherComp 0.908 0.054 16.818 0.000 6.032 0.860
##
## Covariances:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## Attention ~~
## ZoomFatigue -10.122 3.604 -2.809 0.005 -0.222 -0.222
## QOL 5.553 3.506 1.584 0.113 0.122 0.122
## Mobility 19.500 3.984 4.895 0.000 0.410 0.410
## TherAlli 10.358 3.602 2.876 0.004 0.227 0.227
## ZoomFatigue ~~
## QOL -10.504 3.478 -3.020 0.003 -0.239 -0.239
## Mobility -16.934 3.785 -4.474 0.000 -0.369 -0.369
## TherAlli -15.890 3.637 -4.369 0.000 -0.361 -0.361
## QOL ~~
## Mobility 21.379 3.920 5.453 0.000 0.466 0.466
## TherAlli 26.020 4.010 6.488 0.000 0.591 0.591
## Mobility ~~
## TherAlli 0.731 3.498 0.209 0.835 0.016 0.016
##
## Variances:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## .Fixadur 12.098 1.918 6.308 0.000 12.098 0.203
## .Fixanum 7.792 1.786 4.363 0.000 7.792 0.131
## .Fixlocs 16.616 2.043 8.132 0.000 16.616 0.311
## .Energy 10.294 1.805 5.703 0.000 10.294 0.190
## .FAS 10.227 1.936 5.283 0.000 10.227 0.173
## .FSS 17.332 2.154 8.046 0.000 17.332 0.322
## .HRQOL 10.807 1.551 6.967 0.000 10.807 0.197
## .BDI2 8.312 1.514 5.492 0.000 8.312 0.139
## .GAD7 12.565 1.643 7.650 0.000 12.565 0.236
## .SixMinWalk 12.607 1.757 7.174 0.000 12.607 0.209
## .TenMWalk 8.104 1.495 5.420 0.000 8.104 0.138
## .BalConf 12.659 1.726 7.333 0.000 12.659 0.218
## .Disclose 11.599 1.691 6.860 0.000 11.599 0.208
## .Hopeful 8.598 1.555 5.528 0.000 8.598 0.153
## .TherComp 12.854 1.663 7.729 0.000 12.854 0.261
## Attention 47.404 6.013 7.884 0.000 1.000 1.000
## ZoomFatigue 44.009 5.534 7.952 0.000 1.000 1.000
## QOL 43.979 5.485 8.018 0.000 1.000 1.000
## Mobility 47.766 6.030 7.922 0.000 1.000 1.000
## TherAlli 44.097 5.585 7.896 0.000 1.000 1.000
fitmeasures(fit.Wk10)
## npar fmin chisq df
## 40.000 0.190 76.026 80.000
## pvalue baseline.chisq baseline.df baseline.pvalue
## 0.605 2556.997 105.000 0.000
## cfi tli nnfi rfi
## 1.000 1.002 1.002 0.961
## nfi pnfi ifi rni
## 0.970 0.739 1.002 1.002
## logl unrestricted.logl aic bic
## -9062.990 -9024.977 18205.980 18337.912
## ntotal bic2 rmsea rmsea.ci.lower
## 200.000 18211.188 0.000 0.000
## rmsea.ci.upper rmsea.pvalue rmr rmr_nomean
## 0.035 0.997 1.575 1.575
## srmr srmr_bentler srmr_bentler_nomean crmr
## 0.028 0.028 0.028 0.030
## crmr_nomean srmr_mplus srmr_mplus_nomean cn_05
## 0.030 0.028 0.028 269.013
## cn_01 gfi agfi pgfi
## 296.502 0.953 0.930 0.636
## mfi ecvi
## 1.010 0.780
#install.packages("lavaanPlot")
library(lavaanPlot)
## Warning: package 'lavaanPlot' was built under R version 4.0.3
lavaanPlot(model = fit.Wk10, sig=.05,stars=c("regress","latent","covs"),
node_options = list(shape = "box", fontname = "Arial"),
edge_options = list(color = "black"),
coefs = TRUE, covs = TRUE, stand=TRUE)
#This is the SEM based on the same data
#Includes bootstrapped tests of indirect effects
base.mod<- '
# measurement model
Attention =~ Fixadur + Fixanum + Fixlocs
ZoomFatigue =~Energy + FAS + FSS
QOL =~ HRQOL + BDI2 + GAD7
Mobility =~ SixMinWalk + TenMWalk + BalConf
TherAlli =~ Disclose + Hopeful+ TherComp
# regressions
TherAlli ~ a1*Attention + a2*ZoomFatigue + a3*Mobility
QOL ~ b1*TherAlli + Attention + ZoomFatigue + Mobility
# exogenous and mediator correlations
Attention ~~ ZoomFatigue + Mobility
ZoomFatigue ~~ Mobility
# indirect effects
ind_Attention_TherAlli := a1*b1
ind_ZoomFatigue_TherAlli := a2*b1
ind_Mobility_TherAlli := a3*b1
total_ind := ind_Attention_TherAlli + ind_ZoomFatigue_TherAlli + ind_Mobility_TherAlli
'
base.out <- sem(base.mod, data = Wk12[1:15], meanstructure = TRUE,
se = "boot", bootstrap = 1000)
summary(base.out, fit.measures=TRUE, standardized=TRUE)
## lavaan 0.6-7 ended normally after 146 iterations
##
## Estimator ML
## Optimization method NLMINB
## Number of free parameters 55
##
## Number of observations 200
##
## Model Test User Model:
##
## Test statistic 76.026
## Degrees of freedom 80
## P-value (Chi-square) 0.605
##
## Model Test Baseline Model:
##
## Test statistic 2556.997
## Degrees of freedom 105
## P-value 0.000
##
## User Model versus Baseline Model:
##
## Comparative Fit Index (CFI) 1.000
## Tucker-Lewis Index (TLI) 1.002
##
## Loglikelihood and Information Criteria:
##
## Loglikelihood user model (H0) -9062.990
## Loglikelihood unrestricted model (H1) -9024.977
##
## Akaike (AIC) 18235.980
## Bayesian (BIC) 18417.387
## Sample-size adjusted Bayesian (BIC) 18243.141
##
## Root Mean Square Error of Approximation:
##
## RMSEA 0.000
## 90 Percent confidence interval - lower 0.000
## 90 Percent confidence interval - upper 0.035
## P-value RMSEA <= 0.05 0.997
##
## Standardized Root Mean Square Residual:
##
## SRMR 0.026
##
## Parameter Estimates:
##
## Standard errors Bootstrap
## Number of requested bootstrap draws 1000
## Number of successful bootstrap draws 1000
##
## Latent Variables:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## Attention =~
## Fixadur 1.000 6.886 0.893
## Fixanum 1.044 0.063 16.490 0.000 7.189 0.932
## Fixlocs 0.880 0.058 15.263 0.000 6.062 0.830
## ZoomFatigue =~
## Energy 1.000 6.634 0.900
## FAS 1.054 0.061 17.358 0.000 6.993 0.909
## FSS 0.911 0.058 15.818 0.000 6.046 0.824
## QOL =~
## HRQOL 1.000 6.632 0.896
## BDI2 1.083 0.055 19.854 0.000 7.182 0.928
## GAD7 0.960 0.058 16.645 0.000 6.370 0.874
## Mobility =~
## SixMinWalk 1.000 6.911 0.889
## TenMWalk 1.031 0.050 20.588 0.000 7.127 0.929
## BalConf 0.976 0.055 17.823 0.000 6.748 0.885
## TherAlli =~
## Disclose 1.000 6.639 0.890
## Hopeful 1.041 0.057 18.335 0.000 6.909 0.921
## TherComp 0.908 0.056 16.292 0.000 6.032 0.860
##
## Regressions:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## TherAlli ~
## Attention (a1) 0.223 0.086 2.591 0.010 0.231 0.231
## ZoomFatig (a2) -0.392 0.076 -5.138 0.000 -0.392 -0.392
## Mobility (a3) -0.215 0.076 -2.827 0.005 -0.224 -0.224
## QOL ~
## TherAlli (b1) 0.707 0.069 10.205 0.000 0.708 0.708
## Attention -0.245 0.064 -3.857 0.000 -0.255 -0.255
## ZoomFatig 0.193 0.067 2.865 0.004 0.193 0.193
## Mobility 0.605 0.076 7.932 0.000 0.631 0.631
##
## Covariances:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## Attention ~~
## ZoomFatigue -10.121 3.728 -2.715 0.007 -0.222 -0.222
## Mobility 19.499 3.798 5.134 0.000 0.410 0.410
## ZoomFatigue ~~
## Mobility -16.936 3.695 -4.584 0.000 -0.369 -0.369
##
## Intercepts:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## .Fixadur 50.168 0.537 93.397 0.000 50.168 6.503
## .Fixanum 50.118 0.537 93.350 0.000 50.118 6.499
## .Fixlocs 49.333 0.507 97.237 0.000 49.333 6.753
## .Energy 49.967 0.530 94.317 0.000 49.967 6.781
## .FAS 49.845 0.563 88.505 0.000 49.845 6.482
## .FSS 50.401 0.536 94.104 0.000 50.401 6.866
## .HRQOL 49.803 0.519 95.939 0.000 49.803 6.728
## .BDI2 50.091 0.532 94.120 0.000 50.091 6.473
## .GAD7 49.828 0.507 98.238 0.000 49.828 6.835
## .SixMinWalk 50.149 0.556 90.222 0.000 50.149 6.454
## .TenMWalk 50.024 0.541 92.456 0.000 50.024 6.518
## .BalConf 50.298 0.538 93.488 0.000 50.298 6.594
## .Disclose 49.948 0.538 92.892 0.000 49.948 6.694
## .Hopeful 50.159 0.528 94.999 0.000 50.159 6.683
## .TherComp 49.568 0.491 100.918 0.000 49.568 7.064
## Attention 0.000 0.000 0.000
## ZoomFatigue 0.000 0.000 0.000
## .QOL 0.000 0.000 0.000
## Mobility 0.000 0.000 0.000
## .TherAlli 0.000 0.000 0.000
##
## Variances:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## .Fixadur 12.097 2.036 5.941 0.000 12.097 0.203
## .Fixanum 7.792 2.051 3.799 0.000 7.792 0.131
## .Fixlocs 16.616 2.043 8.134 0.000 16.616 0.311
## .Energy 10.294 1.908 5.395 0.000 10.294 0.190
## .FAS 10.227 1.966 5.202 0.000 10.227 0.173
## .FSS 17.332 2.339 7.412 0.000 17.332 0.322
## .HRQOL 10.806 1.553 6.956 0.000 10.806 0.197
## .BDI2 8.313 1.669 4.982 0.000 8.313 0.139
## .GAD7 12.565 1.712 7.340 0.000 12.565 0.236
## .SixMinWalk 12.607 2.023 6.230 0.000 12.607 0.209
## .TenMWalk 8.104 1.204 6.733 0.000 8.104 0.138
## .BalConf 12.659 1.681 7.532 0.000 12.659 0.218
## .Disclose 11.600 1.735 6.687 0.000 11.600 0.208
## .Hopeful 8.598 1.603 5.362 0.000 8.598 0.153
## .TherComp 12.854 1.848 6.955 0.000 12.854 0.261
## Attention 47.412 6.096 7.778 0.000 1.000 1.000
## ZoomFatigue 44.007 5.893 7.468 0.000 1.000 1.000
## .QOL 16.029 2.259 7.095 0.000 0.364 0.364
## Mobility 47.765 5.638 8.472 0.000 1.000 1.000
## .TherAlli 35.697 4.250 8.399 0.000 0.810 0.810
##
## Defined Parameters:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## ind_Attntn_ThA 0.158 0.064 2.462 0.014 0.164 0.164
## ind_ZmFtg_ThrA -0.278 0.060 -4.629 0.000 -0.278 -0.278
## ind_Mblty_ThrA -0.152 0.057 -2.656 0.008 -0.158 -0.158
## total_ind -0.272 0.096 -2.822 0.005 -0.272 -0.272
parameterEstimates(base.out, se = TRUE, zstat = TRUE, pvalue = TRUE,
ci = TRUE, level = 0.95, boot.ci.type = "bca.simple",
standardized = TRUE, cov.std = FALSE, fmi = FALSE,
remove.system.eq = TRUE, remove.eq = TRUE,
remove.ineq = TRUE, remove.def = FALSE,
rsquare = TRUE, add.attributes = FALSE, header = TRUE)
## lhs op
## 1 Attention =~
## 2 Attention =~
## 3 Attention =~
## 4 ZoomFatigue =~
## 5 ZoomFatigue =~
## 6 ZoomFatigue =~
## 7 QOL =~
## 8 QOL =~
## 9 QOL =~
## 10 Mobility =~
## 11 Mobility =~
## 12 Mobility =~
## 13 TherAlli =~
## 14 TherAlli =~
## 15 TherAlli =~
## 16 TherAlli ~
## 17 TherAlli ~
## 18 TherAlli ~
## 19 QOL ~
## 20 QOL ~
## 21 QOL ~
## 22 QOL ~
## 23 Attention ~~
## 24 Attention ~~
## 25 ZoomFatigue ~~
## 26 Fixadur ~~
## 27 Fixanum ~~
## 28 Fixlocs ~~
## 29 Energy ~~
## 30 FAS ~~
## 31 FSS ~~
## 32 HRQOL ~~
## 33 BDI2 ~~
## 34 GAD7 ~~
## 35 SixMinWalk ~~
## 36 TenMWalk ~~
## 37 BalConf ~~
## 38 Disclose ~~
## 39 Hopeful ~~
## 40 TherComp ~~
## 41 Attention ~~
## 42 ZoomFatigue ~~
## 43 QOL ~~
## 44 Mobility ~~
## 45 TherAlli ~~
## 46 Fixadur ~1
## 47 Fixanum ~1
## 48 Fixlocs ~1
## 49 Energy ~1
## 50 FAS ~1
## 51 FSS ~1
## 52 HRQOL ~1
## 53 BDI2 ~1
## 54 GAD7 ~1
## 55 SixMinWalk ~1
## 56 TenMWalk ~1
## 57 BalConf ~1
## 58 Disclose ~1
## 59 Hopeful ~1
## 60 TherComp ~1
## 61 Attention ~1
## 62 ZoomFatigue ~1
## 63 QOL ~1
## 64 Mobility ~1
## 65 TherAlli ~1
## 66 ind_Attention_TherAlli :=
## 67 ind_ZoomFatigue_TherAlli :=
## 68 ind_Mobility_TherAlli :=
## 69 total_ind :=
## 70 Fixadur r2
## 71 Fixanum r2
## 72 Fixlocs r2
## 73 Energy r2
## 74 FAS r2
## 75 FSS r2
## 76 HRQOL r2
## 77 BDI2 r2
## 78 GAD7 r2
## 79 SixMinWalk r2
## 80 TenMWalk r2
## 81 BalConf r2
## 82 Disclose r2
## 83 Hopeful r2
## 84 TherComp r2
## 85 QOL r2
## 86 TherAlli r2
## rhs
## 1 Fixadur
## 2 Fixanum
## 3 Fixlocs
## 4 Energy
## 5 FAS
## 6 FSS
## 7 HRQOL
## 8 BDI2
## 9 GAD7
## 10 SixMinWalk
## 11 TenMWalk
## 12 BalConf
## 13 Disclose
## 14 Hopeful
## 15 TherComp
## 16 Attention
## 17 ZoomFatigue
## 18 Mobility
## 19 TherAlli
## 20 Attention
## 21 ZoomFatigue
## 22 Mobility
## 23 ZoomFatigue
## 24 Mobility
## 25 Mobility
## 26 Fixadur
## 27 Fixanum
## 28 Fixlocs
## 29 Energy
## 30 FAS
## 31 FSS
## 32 HRQOL
## 33 BDI2
## 34 GAD7
## 35 SixMinWalk
## 36 TenMWalk
## 37 BalConf
## 38 Disclose
## 39 Hopeful
## 40 TherComp
## 41 Attention
## 42 ZoomFatigue
## 43 QOL
## 44 Mobility
## 45 TherAlli
## 46
## 47
## 48
## 49
## 50
## 51
## 52
## 53
## 54
## 55
## 56
## 57
## 58
## 59
## 60
## 61
## 62
## 63
## 64
## 65
## 66 a1*b1
## 67 a2*b1
## 68 a3*b1
## 69 ind_Attention_TherAlli+ind_ZoomFatigue_TherAlli+ind_Mobility_TherAlli
## 70 Fixadur
## 71 Fixanum
## 72 Fixlocs
## 73 Energy
## 74 FAS
## 75 FSS
## 76 HRQOL
## 77 BDI2
## 78 GAD7
## 79 SixMinWalk
## 80 TenMWalk
## 81 BalConf
## 82 Disclose
## 83 Hopeful
## 84 TherComp
## 85 QOL
## 86 TherAlli
## label est se z pvalue ci.lower ci.upper
## 1 1.000 0.000 NA NA 1.000 1.000
## 2 1.044 0.063 16.490 0.000 0.933 1.171
## 3 0.880 0.058 15.263 0.000 0.773 1.000
## 4 1.000 0.000 NA NA 1.000 1.000
## 5 1.054 0.061 17.358 0.000 0.939 1.177
## 6 0.911 0.058 15.818 0.000 0.807 1.043
## 7 1.000 0.000 NA NA 1.000 1.000
## 8 1.083 0.055 19.854 0.000 0.977 1.189
## 9 0.960 0.058 16.645 0.000 0.850 1.074
## 10 1.000 0.000 NA NA 1.000 1.000
## 11 1.031 0.050 20.588 0.000 0.933 1.128
## 12 0.976 0.055 17.823 0.000 0.868 1.086
## 13 1.000 0.000 NA NA 1.000 1.000
## 14 1.041 0.057 18.335 0.000 0.936 1.153
## 15 0.908 0.056 16.292 0.000 0.800 1.012
## 16 a1 0.223 0.086 2.591 0.010 0.054 0.389
## 17 a2 -0.392 0.076 -5.138 0.000 -0.558 -0.259
## 18 a3 -0.215 0.076 -2.827 0.005 -0.368 -0.065
## 19 b1 0.707 0.069 10.205 0.000 0.570 0.844
## 20 -0.245 0.064 -3.857 0.000 -0.369 -0.120
## 21 0.193 0.067 2.865 0.004 0.066 0.322
## 22 0.605 0.076 7.932 0.000 0.455 0.752
## 23 -10.121 3.728 -2.715 0.007 -18.015 -3.571
## 24 19.499 3.798 5.134 0.000 13.027 28.198
## 25 -16.936 3.695 -4.584 0.000 -24.379 -9.727
## 26 12.097 2.036 5.941 0.000 8.358 16.177
## 27 7.792 2.051 3.799 0.000 4.010 11.986
## 28 16.616 2.043 8.134 0.000 12.906 21.150
## 29 10.294 1.908 5.395 0.000 6.920 14.536
## 30 10.227 1.966 5.202 0.000 6.901 14.157
## 31 17.332 2.339 7.412 0.000 13.356 22.806
## 32 10.806 1.553 6.956 0.000 7.911 13.844
## 33 8.313 1.669 4.982 0.000 5.359 11.742
## 34 12.565 1.712 7.340 0.000 9.394 16.115
## 35 12.607 2.023 6.230 0.000 8.884 16.582
## 36 8.104 1.204 6.733 0.000 5.941 10.546
## 37 12.659 1.681 7.532 0.000 9.672 16.372
## 38 11.600 1.735 6.687 0.000 8.344 15.219
## 39 8.598 1.603 5.362 0.000 5.761 12.129
## 40 12.854 1.848 6.955 0.000 9.634 17.007
## 41 47.412 6.096 7.778 0.000 36.195 60.174
## 42 44.007 5.893 7.468 0.000 34.011 56.991
## 43 16.029 2.259 7.095 0.000 12.158 21.649
## 44 47.765 5.638 8.472 0.000 37.555 60.350
## 45 35.697 4.250 8.399 0.000 28.206 46.122
## 46 50.168 0.537 93.397 0.000 49.177 51.273
## 47 50.118 0.537 93.350 0.000 49.081 51.198
## 48 49.333 0.507 97.237 0.000 48.431 50.368
## 49 49.967 0.530 94.317 0.000 48.874 50.956
## 50 49.845 0.563 88.505 0.000 48.708 50.921
## 51 50.401 0.536 94.104 0.000 49.333 51.501
## 52 49.803 0.519 95.939 0.000 48.822 50.865
## 53 50.091 0.532 94.120 0.000 49.036 51.076
## 54 49.828 0.507 98.238 0.000 48.805 50.803
## 55 50.149 0.556 90.222 0.000 49.033 51.201
## 56 50.024 0.541 92.456 0.000 48.873 51.029
## 57 50.298 0.538 93.488 0.000 49.253 51.407
## 58 49.948 0.538 92.892 0.000 48.893 51.009
## 59 50.159 0.528 94.999 0.000 49.184 51.268
## 60 49.568 0.491 100.918 0.000 48.633 50.523
## 61 0.000 0.000 NA NA 0.000 0.000
## 62 0.000 0.000 NA NA 0.000 0.000
## 63 0.000 0.000 NA NA 0.000 0.000
## 64 0.000 0.000 NA NA 0.000 0.000
## 65 0.000 0.000 NA NA 0.000 0.000
## 66 ind_Attention_TherAlli 0.158 0.064 2.462 0.014 0.037 0.285
## 67 ind_ZoomFatigue_TherAlli -0.278 0.060 -4.629 0.000 -0.425 -0.181
## 68 ind_Mobility_TherAlli -0.152 0.057 -2.656 0.008 -0.268 -0.045
## 69 total_ind -0.272 0.096 -2.822 0.005 -0.487 -0.099
## 70 0.797 NA NA NA NA NA
## 71 0.869 NA NA NA NA NA
## 72 0.689 NA NA NA NA NA
## 73 0.810 NA NA NA NA NA
## 74 0.827 NA NA NA NA NA
## 75 0.678 NA NA NA NA NA
## 76 0.803 NA NA NA NA NA
## 77 0.861 NA NA NA NA NA
## 78 0.764 NA NA NA NA NA
## 79 0.791 NA NA NA NA NA
## 80 0.862 NA NA NA NA NA
## 81 0.782 NA NA NA NA NA
## 82 0.792 NA NA NA NA NA
## 83 0.847 NA NA NA NA NA
## 84 0.739 NA NA NA NA NA
## 85 0.636 NA NA NA NA NA
## 86 0.190 NA NA NA NA NA
## std.lv std.all std.nox
## 1 6.886 0.893 0.893
## 2 7.189 0.932 0.932
## 3 6.062 0.830 0.830
## 4 6.634 0.900 0.900
## 5 6.993 0.909 0.909
## 6 6.046 0.824 0.824
## 7 6.632 0.896 0.896
## 8 7.182 0.928 0.928
## 9 6.370 0.874 0.874
## 10 6.911 0.889 0.889
## 11 7.127 0.929 0.929
## 12 6.748 0.885 0.885
## 13 6.639 0.890 0.890
## 14 6.909 0.921 0.921
## 15 6.032 0.860 0.860
## 16 0.231 0.231 0.231
## 17 -0.392 -0.392 -0.392
## 18 -0.224 -0.224 -0.224
## 19 0.708 0.708 0.708
## 20 -0.255 -0.255 -0.255
## 21 0.193 0.193 0.193
## 22 0.631 0.631 0.631
## 23 -0.222 -0.222 -0.222
## 24 0.410 0.410 0.410
## 25 -0.369 -0.369 -0.369
## 26 12.097 0.203 0.203
## 27 7.792 0.131 0.131
## 28 16.616 0.311 0.311
## 29 10.294 0.190 0.190
## 30 10.227 0.173 0.173
## 31 17.332 0.322 0.322
## 32 10.806 0.197 0.197
## 33 8.313 0.139 0.139
## 34 12.565 0.236 0.236
## 35 12.607 0.209 0.209
## 36 8.104 0.138 0.138
## 37 12.659 0.218 0.218
## 38 11.600 0.208 0.208
## 39 8.598 0.153 0.153
## 40 12.854 0.261 0.261
## 41 1.000 1.000 1.000
## 42 1.000 1.000 1.000
## 43 0.364 0.364 0.364
## 44 1.000 1.000 1.000
## 45 0.810 0.810 0.810
## 46 50.168 6.503 6.503
## 47 50.118 6.499 6.499
## 48 49.333 6.753 6.753
## 49 49.967 6.781 6.781
## 50 49.845 6.482 6.482
## 51 50.401 6.866 6.866
## 52 49.803 6.728 6.728
## 53 50.091 6.473 6.473
## 54 49.828 6.835 6.835
## 55 50.149 6.454 6.454
## 56 50.024 6.518 6.518
## 57 50.298 6.594 6.594
## 58 49.948 6.694 6.694
## 59 50.159 6.683 6.683
## 60 49.568 7.064 7.064
## 61 0.000 0.000 0.000
## 62 0.000 0.000 0.000
## 63 0.000 0.000 0.000
## 64 0.000 0.000 0.000
## 65 0.000 0.000 0.000
## 66 0.164 0.164 0.164
## 67 -0.278 -0.278 -0.278
## 68 -0.158 -0.158 -0.158
## 69 -0.272 -0.272 -0.272
## 70 NA NA NA
## 71 NA NA NA
## 72 NA NA NA
## 73 NA NA NA
## 74 NA NA NA
## 75 NA NA NA
## 76 NA NA NA
## 77 NA NA NA
## 78 NA NA NA
## 79 NA NA NA
## 80 NA NA NA
## 81 NA NA NA
## 82 NA NA NA
## 83 NA NA NA
## 84 NA NA NA
## 85 NA NA NA
## 86 NA NA NA
fitmeasures(base.out)
## npar fmin chisq df
## 55.000 0.190 76.026 80.000
## pvalue baseline.chisq baseline.df baseline.pvalue
## 0.605 2556.997 105.000 0.000
## cfi tli nnfi rfi
## 1.000 1.002 1.002 0.961
## nfi pnfi ifi rni
## 0.970 0.739 1.002 1.002
## logl unrestricted.logl aic bic
## -9062.990 -9024.977 18235.980 18417.387
## ntotal bic2 rmsea rmsea.ci.lower
## 200.000 18243.141 0.000 0.000
## rmsea.ci.upper rmsea.pvalue rmr rmr_nomean
## 0.035 0.997 1.485 1.575
## srmr srmr_bentler srmr_bentler_nomean crmr
## 0.026 0.026 0.028 0.028
## crmr_nomean srmr_mplus srmr_mplus_nomean cn_05
## 0.030 0.026 0.028 269.013
## cn_01 gfi agfi pgfi
## 296.502 0.999 0.998 0.592
## mfi ecvi
## 1.010 0.930