• 1000 все гетеросексуалы) удваиваются и хотят проверить неизменность измерений между полами. Кроме того, я хочу проверить инвариантность для другого фактора группирования (
cohort между диадами). Чтобы учесть предполагаемую несамостоятельность ошибок пар, я рассматриваю пол как повторяющееся измерение в диадах. Да, кстати, все элементы двоичные.
Мой код выглядит так: (Это не мой реальный код, я упростил имена переменных et c. К сожалению, я не знаю, как сгенерировать данные dyadi c в качестве воспроизводимого примера.)
mod.1 <- "facA_m =~ itemA1_m + itemA2_m + itemA3_m + itemA4_m + itemA5_m + ... #12 items
facA_f =~ itemA1_f + itemA2_f + itemA3_f + itemA4_f + itemA5_f + ... #12 items
facB_m =~ itemB1_m + itemB2_m + itemB3_m + itemB4_m + itemB5_m + ... #12 items
facB_f =~ itemB1_f + itemB2_f + itemB3_f + itemB4_f + itemB5_f + ... #12 items"
genderAsLong <- list(facA = c("facA_m", "facA_f"),
facB = c("facB_m", "facB_f"))
gencoh.config <- measEq.syntax(configural.model = mod.1, data = dyadicd,
ordered = c( paste0("itemA", c(1:12), "_m"),
paste0("itemB", c(1:12), "_f") ),
longFacNames = genderAsLong, group = "cohort",
parameterization = "theta", estimator = "WLSMV",
ID.fac = "std.lv", ID.cat = "Wu.Estabrook.2016",
return.fit = TRUE)
Я получаю следующее сообщение об ошибке: Error in if (any(dCOV < 0)) { : missing value where TRUE/FALSE needed
О чем конкретно говорится в этом сообщении, это из-за пропущенных значений (у меня их несколько, не проблема в других аналогичных моделях) или возникла техническая проблема? Было бы очень полезно уже знать, является ли это проблемой кодирования вообще или статистической проблемой. Заранее благодарим за любой совет.