У меня есть набор данных, который выглядит следующим образом:
id v1 v2 v3 time
s01 4 3 1 0
s02 2 3 1 0
...
s01 4 5 4 1
s02 3 1 4 1
Моя спецификация модели Lavaan выглядит следующим образом:
model <- '
lat =~ v1 + v2 + v3
'
, и я хочу измерить различные инвариантности измерения по всей groups (= моменты времени: 1/2).
Следующий код прекрасно работает, но я не думаю, что он распознает переменную id, поэтому группы рассматриваются как независимые (например, пол), но я хочу, чтобы они оценивались как зависимые. Можете ли вы помочь мне с этим? Есть ли способ определить переменную ID в Lavaan?
MI_conf <- measEq.syntax(model, data = data_long , group="time", return.fit = T, missing = "fiml")
MI_conf_sum <- summary(MI_conf, fit.measures = TRUE)
MI_weak <- measEq.syntax(model, data = data_long , group="time", return.fit = T, missing = "fiml",
group.equal=c("loadings"))
MI_weak_sum <- summary(MI_weak, fit.measures = TRUE)
MI_strong <- measEq.syntax(model, data = data_long , group="time", return.fit = T, missing = "fiml",
group.equal=c("loadings", "intercepts"))
MI_strong_sum <- summary(MI_strong, fit.measures = TRUE)
MI_strict <- measEq.syntax(model, data = data_long , group="time", return.fit = T, missing = "fiml",
group.equal=c("loadings", "intercepts", "residuals"))
MI_strict_sum <- summary(MI_strict, fit.measures = TRUE)