D<-read.table("SCHOOLPLEZIER.dat",header=TRUE)
install.packages("lavaan")
install.packages("stringr")
install.packages("recoder")
library(lavaan)
library(stringr)
library(recoder)
D$item1<-recoder(D$item1,"1:4;2:3;3:2;4:1")
D$item2<-recoder(D$item2,"1:4;2:3;3:2;4:1")
D$item3<-recoder(D$item3,"1:4;2:3;3:2;4:1")
D$item4<-recoder(D$item4,"1:4;2:3;3:2;4:1")
D$item5<-recoder(D$item5,"1:4;2:3;3:2;4:1")
D$item6<-recoder(D$item6,"1:4;2:3;3:2;4:1")
model1 <- "
f1 =~ item1 + item2 + item4 + item7
f2 =~ item3 + item5 + item6 + item8"
fit <- cfa(model = model1, data = D, std.lv = TRUE)
summary(fit, fit.measures = TRUE)
Theta<-inspect(fit,"est")$theta
Sigma<-fitted(fit)$cov
sum(diag(Sigma-Theta))/sum(diag(Sigma))
L <- inspect(fit,"est")$lambda # dit is de geschatte factorladingenmatrix
P <- inspect(fit,"est")$psi # dit is de geschatte covariantiematrix van de gemeenschappelijke factoren
install.packages("lavaan")
install.packages("stringr")
install.packages("recoder")
library(lavaan)
library(stringr)
library(recoder)
D$item1<-recoder(D$item1,"1:4;2:3;3:2;4:1")
D$item2<-recoder(D$item2,"1:4;2:3;3:2;4:1")
D$item3<-recoder(D$item3,"1:4;2:3;3:2;4:1")
D$item4<-recoder(D$item4,"1:4;2:3;3:2;4:1")
D$item5<-recoder(D$item5,"1:4;2:3;3:2;4:1")
D$item6<-recoder(D$item6,"1:4;2:3;3:2;4:1")
model1 <- "
f1 =~ item1 + item2 + item4 + item7
f2 =~ item3 + item5 + item6 + item8"
fit <- cfa(model = model1, data = D, std.lv = TRUE)
summary(fit, fit.measures = TRUE)
Theta<-inspect(fit,"est")$theta
Sigma<-fitted(fit)$cov
sum(diag(Sigma-Theta))/sum(diag(Sigma))
L <- inspect(fit,"est")$lambda # dit is de geschatte factorladingenmatrix
P <- inspect(fit,"est")$psi # dit is de geschatte covariantiematrix van de gemeenschappelijke factoren