Рассчитать неизменность измерений для зависимых групп в длинном формате - PullRequest
0 голосов
/ 06 марта 2020

У меня есть набор данных, который выглядит следующим образом:

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)
...