Модель Лаваана с двоичными переменными и модерируемым посредничеством

#simulation

Вопрос:

В настоящее время я пытаюсь указать свою модель, чтобы использовать параметр для моделирования (я опубликовал это здесь).

Теперь я попытался указать свою модель модели следующим образом:

# declare binary IV a factor

 US$unexp <- factor(US$unexp)
 

# order binary DV

 US[,"sign"] <- lapply(US[, "sign"], ordered)
 

# interaction product term (to keep it simple(r): only consider one moderation fc*cn)

 US$cnfc <- US$cn*US$fc_1_o_agr
 

# model specification

 model <- 'fc_1_o_agr ~ a1*unexp
      joy ~ a2*fc_1_o_agr
      sign ~ b1*fc_1_o_agr   b2*joy   b3*cn   b4*cnfc   c*unexp
      # when cn = 0
      a1b30 := a1*b3
      total0 := a1b30   b3   c
      # when cn = -1sd (sd = 1.32)
      a1b3low := a1   b3*-1.32
      totallow := a1b3low   c
      # when cn =  1sd (sd = 1.32)
      a1b3high := a1   b3*1.32
      totalhigh := a1b3high  c'
 

# fit model

 fit <- cfa(model, data = US, ordered = "sign", std.lv = T)
 

Однако, когда я затем бегу
output <- sim(NULL, n = rep(70:300,2), model = fit, generate = fit)

Я получаю следующую ошибку:

Ошибка в МАССЕ::mvrnorm(n = sample.nobs[g], mu = Mu.hat[[g]], Sigma = COV, : «Сигма» не является положительно определенной

Что это значит? И могу ли я что-нибудь с этим поделать? Я не слишком хорошо знаком со слишком техническим языком и был бы признателен, если бы вы объяснили все просто.

Спасибо.