#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, : «Сигма» не является положительно определенной
Что это значит? И могу ли я что-нибудь с этим поделать? Я не слишком хорошо знаком со слишком техническим языком и был бы признателен, если бы вы объяснили все просто.
Спасибо.