Thanks Mikko for the interest in the post.
This is my code:
library(xlsx)
library(lavaan)
library(semPlot)
library(coefficientalpha)
library(psy)
library(semTools)
library(SemPlot)
library(pastecs)
library(corrplot)
library(ggplot2)
data1<-read.xlsx('R_ins2.xlsx', sheetName='Foglio3', header=T)
model1 <- '
# measurement model
var1 =~ I2 +I5 + I6
var2=~ Y1 + Y2 + Y3 + Y4 + Y5 + Y6
var3=~ C1 + C2 + C3 + C4
var4=~ T1 + T2 + T3 + T4
var5=~ H1 + H2 + H3 + H6 + H7+ H9 +I4
var6=~ AW1 + AW2+ AW3+ X6 + X7 + X8
'
m1_fit <- cfa(model1, data=data1)
summary(m1_fit, standardized=TRUE, rsq=T, fit.measures = TRUE)
reliability(m1_fit)
Thanks,
Elena