rm(list = ls())
library(knitr)
library(foreign)
library(lavaan)
## This is lavaan 0.5-22
## lavaan is BETA software! Please report any bugs.
setwd(dir = "/Users/ivanropovik/OneDrive/Projects/APVV IK/DATA")
data <- read.spss("APVV_ANALOGIES_ARTICLE_imputed.sav", to.data.frame=TRUE)
## Warning in read.spss("APVV_ANALOGIES_ARTICLE_imputed.sav", to.data.frame
## = TRUE): APVV_ANALOGIES_ARTICLE_imputed.sav: Unrecognized record type 7,
## subtype 18 encountered in system file
## Warning in read.spss("APVV_ANALOGIES_ARTICLE_imputed.sav", to.data.frame
## = TRUE): APVV_ANALOGIES_ARTICLE_imputed.sav: Unrecognized record type 7,
## subtype 24 encountered in system file
#attach(data)
data$P_STROOP_INH_CAS_RS <- with(data, ((max(P_STROOP_INH_CAS) + 1) - P_STROOP_INH_CAS))
data$P_TMT_stried_RS <- with(data, ((max(P_TMT_stried) + 1) - P_TMT_stried))
data$AL_HINT_SUM_RS <- with(data, ((max(AL_HINT_SUM) + 1) - AL_HINT_SUM))
respc.model <- '
Att =~ a*P_STROOP_INH_CAS + b*WJ_VIZ_POROV
Fluency =~ c*P_FLUE_PISM_SUM + d*P_DF_spravne_spolu
Shifting =~ e*P_TMT_stried + f*P_FLUE_STR_POC
WM =~ g*WJ_CIS_RADY + h*TOH_SCORE
gF =~ i*WJ_PRIEST_VZTAH + j*WJ_KV_VYVODZ
AtL =~ k*Learn + l*AL_HINT_SUM
Fluency ~ r*Att
Shifting ~ s*Att
WM ~ t*Att
gF ~ u*WM
AtL ~ v*WM
WJ_VIZ_POROV ~~ P_STROOP_INH_CAS
gF ~~ AtL
'
fitted.model <- sem(model = respc.model, data = data, meanstructure = FALSE,
std.lv = FALSE, estimator = "ML", test = "standard",
orthogonal = TRUE, std.ov = TRUE, likelihood = "wishart", bootstrap = 2000)
summary(fitted.model, standardized = TRUE, rsquare = TRUE)
## lavaan (0.5-22) converged normally after 39 iterations
##
## Number of observations 210
##
## Estimator ML
## Minimum Function Test Statistic 61.485
## Degrees of freedom 47
## P-value (Chi-square) 0.076
##
## Parameter Estimates:
##
## Information Expected
## Standard Errors Standard
##
## Latent Variables:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## Att =~
## P_STROOP_I (a) 1.000 0.567 0.567
## WJ_VIZ_POR (b) -0.968 0.146 -6.632 0.000 -0.549 -0.549
## Fluency =~
## P_FLUE_PIS (c) 1.000 0.481 0.481
## P_DF_sprv_ (d) 1.034 0.235 4.393 0.000 0.497 0.497
## Shifting =~
## P_TMT_strd (e) 1.000 0.750 0.750
## P_FLUE_STR (f) -0.624 0.121 -5.168 0.000 -0.468 -0.468
## WM =~
## WJ_CIS_RAD (g) 1.000 0.665 0.665
## TOH_SCORE (h) 0.424 0.131 3.235 0.001 0.282 0.282
## gF =~
## WJ_PRIEST_ (i) 1.000 0.575 0.575
## WJ_KV_VYVO (j) 1.365 0.211 6.461 0.000 0.785 0.785
## AtL =~
## Learn (k) 1.000 0.732 0.732
## AL_HINT_SU (l) -1.037 0.171 -6.080 0.000 -0.759 -0.759
##
## Regressions:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## Fluency ~
## Att (r) -0.778 0.167 -4.662 0.000 -0.917 -0.917
## Shifting ~
## Att (s) 1.264 0.215 5.878 0.000 0.955 0.955
## WM ~
## Att (t) -0.752 0.161 -4.658 0.000 -0.640 -0.640
## gF ~
## WM (u) 0.730 0.159 4.603 0.000 0.844 0.844
## AtL ~
## WM (v) 0.538 0.145 3.696 0.000 0.489 0.489
##
## Covariances:
## Estimate Std.Err z-value P(>|z|) Std.lv
## .P_STROOP_INH_CAS ~~
## .WJ_VIZ_POROV -0.189 0.065 -2.934 0.003 -0.189
## .gF ~~
## .AtL 0.093 0.046 2.001 0.045 0.470
## .Fluency ~~
## .Shifting 0.000 0.000
## .gF 0.000 0.000
## .AtL 0.000 0.000
## .Shifting ~~
## .gF 0.000 0.000
## .AtL 0.000 0.000
## Std.all
##
## -0.275
##
## 0.470
##
## 0.000
## 0.000
## 0.000
##
## 0.000
## 0.000
##
## Variances:
## Estimate Std.Err z-value P(>|z|) Std.lv Std.all
## .P_STROOP_INH_C 0.679 0.083 8.169 0.000 0.679 0.679
## .WJ_VIZ_POROV 0.699 0.084 8.332 0.000 0.699 0.699
## .P_FLUE_PISM_SU 0.769 0.097 7.934 0.000 0.769 0.769
## .P_DF_sprvn_spl 0.753 0.098 7.648 0.000 0.753 0.753
## .P_TMT_stried 0.437 0.106 4.108 0.000 0.437 0.437
## .P_FLUE_STR_POC 0.781 0.085 9.158 0.000 0.781 0.781
## .WJ_CIS_RADY 0.557 0.092 6.070 0.000 0.557 0.557
## .TOH_SCORE 0.920 0.094 9.830 0.000 0.920 0.920
## .WJ_PRIEST_VZTA 0.669 0.079 8.512 0.000 0.669 0.669
## .WJ_KV_VYVODZ 0.383 0.089 4.286 0.000 0.383 0.383
## .Learn 0.464 0.092 5.037 0.000 0.464 0.464
## .AL_HINT_SUM 0.424 0.096 4.431 0.000 0.424 0.424
## Att 0.321 0.087 3.671 0.000 1.000 1.000
## .Fluency 0.037 0.061 0.607 0.544 0.160 0.160
## .Shifting 0.050 0.106 0.471 0.638 0.089 0.089
## .WM 0.261 0.077 3.376 0.001 0.590 0.590
## .gF 0.095 0.054 1.769 0.077 0.287 0.287
## .AtL 0.408 0.100 4.064 0.000 0.761 0.761
##
## R-Square:
## Estimate
## P_STROOP_INH_C 0.321
## WJ_VIZ_POROV 0.301
## P_FLUE_PISM_SU 0.231
## P_DF_sprvn_spl 0.247
## P_TMT_stried 0.563
## P_FLUE_STR_POC 0.219
## WJ_CIS_RADY 0.443
## TOH_SCORE 0.080
## WJ_PRIEST_VZTA 0.331
## WJ_KV_VYVODZ 0.617
## Learn 0.536
## AL_HINT_SUM 0.576
## Fluency 0.840
## Shifting 0.911
## WM 0.410
## gF 0.713
## AtL 0.239
fitMeasures(fitted.model)
## npar fmin chisq
## 31.000 0.147 61.485
## df pvalue baseline.chisq
## 47.000 0.076 543.661
## baseline.df baseline.pvalue cfi
## 66.000 0.000 0.970
## tli nnfi rfi
## 0.957 0.957 0.841
## nfi pnfi ifi
## 0.887 0.632 0.971
## rni logl unrestricted.logl
## 0.970 -3333.484 -3302.594
## aic bic ntotal
## 6728.967 6832.580 210.000
## bic2 rmsea rmsea.ci.lower
## 6734.355 0.038 0.000
## rmsea.ci.upper rmsea.pvalue rmr
## 0.063 0.758 0.048
## rmr_nomean srmr srmr_bentler
## 0.048 0.048 0.048
## srmr_bentler_nomean srmr_bollen srmr_bollen_nomean
## 0.048 0.048 0.048
## srmr_mplus srmr_mplus_nomean cn_05
## 0.048 0.048 218.551
## cn_01 gfi agfi
## 247.248 0.953 0.922
## pgfi mfi ecvi
## 0.574 0.966 0.591
residuals(fitted.model, type = "cor")$cor
## P_STRO WJ_VIZ P_FLUE_P P_DF__ P_TMT_ P_FLUE_S WJ_CIS
## P_STROOP_INH_CAS 0.000
## WJ_VIZ_POROV 0.000 0.000
## P_FLUE_PISM_SUM 0.029 -0.001 0.000
## P_DF_spravne_spolu -0.030 -0.020 0.000 0.000
## P_TMT_stried 0.009 -0.005 0.036 -0.008 0.000
## P_FLUE_STR_POC -0.062 -0.043 0.033 -0.014 0.000 0.000
## WJ_CIS_RADY -0.001 0.066 0.075 -0.071 0.037 -0.026 0.000
## TOH_SCORE -0.092 0.039 -0.111 -0.015 -0.152 0.096 -0.041
## WJ_PRIEST_VZTAH 0.067 0.004 0.035 0.041 0.042 -0.075 0.036
## WJ_KV_VYVODZ 0.152 -0.005 0.087 0.022 -0.021 -0.051 0.008
## Learn -0.014 0.060 0.087 0.049 -0.056 -0.022 0.034
## AL_HINT_SUM -0.005 0.077 -0.043 -0.012 0.010 -0.040 0.052
## TOH_SC WJ_PRI WJ_KV_ Learn AL_HIN
## P_STROOP_INH_CAS
## WJ_VIZ_POROV
## P_FLUE_PISM_SUM
## P_DF_spravne_spolu
## P_TMT_stried
## P_FLUE_STR_POC
## WJ_CIS_RADY
## TOH_SCORE 0.000
## WJ_PRIEST_VZTAH -0.087 0.000
## WJ_KV_VYVODZ 0.022 0.000 0.000
## Learn -0.030 -0.017 -0.010 0.000
## AL_HINT_SUM 0.040 -0.029 -0.003 0.000 0.000