Когда я запускаю модель, используя функцию sem из semTools, используя набор данных с пропущенными значениями, все работает достаточно хорошо.И когда я использую пакет мышей для вменения пропущенных значений, импутация работает, я получаю 5 наборов данных без пропущенных данных.
Используемая функция sem выглядит следующим образом:
sem_full <- sem(model=full_model, data=my_data, ordered=c("r01","r02","r03","r04","r05","r06","r07","r08","r09","r10","r11","r12","r13","r14","r17","r15","r16","r19","r20","r21","r22","r23","r24","q01","q02","q03","q04","q05","q09","q11","q13","q15","q22","q23","q24","q17","q18","q20","q10","q12","q14","q19","q21","q25","q26","q27"))
и модельуказано так:
full_model <- '
first_one =~ NA*r01 + r02 + r03 + r04 + r05 + r06 + r07 + + r08 + r09 + r10 + r11 + r12 + r13 + r14 + r17
first_two =~ NA*r13 + r14 + r15 + r16 + r19 + r20 + r21 + r22 + r23 + r24
f1 =~ NA*q01 + q02 + q03 + q04 + q05
f2 =~ NA*q09 + q11 + q13 + q15
f3 =~ NA*q22 + q23 + q24
f4 =~ NA*q17 + q18 + q20
f5 =~ NA*q10 + q12 + q14
f6 =~ NA*q19 + q21
f7 =~ NA*q25 + q26 + q27
first_one + first_two ~ f1
first_one + first_two ~ f2
first_one + first_two ~ f3
first_one + first_two ~ f4
first_one + first_two ~ f5
first_one + first_two ~ f6
first_one + first_two ~ f7
first_one + first_two ~ sex
first_one + first_two ~ age
first_one + first_two ~ freq
first_one + first_two ~ sat
f1 ~~ 1*f1
f2 ~~ 1*f2
f3 ~~ 1*f3
f4 ~~ 1*f4
f5 ~~ 1*f5
f6 ~~ 1*f6
f7 ~~ 1*f7
first_one ~~ 1*first_one
first_two ~~ 1*first_two'
Теперь, когда я использую
sem_mice <- runMI(full_model, data=mice_sem.imp,fun="sem", ordered=c("r01","r02","r03","r04","r05","r06","r07","r08","r09","r10","r11","r12","r13","r14","r17","r15","r16","r19","r20","r21","r22","r23","r24","q01","q02","q03","q04","q05","q09","q11","q13","q15","q22","q23","q24","q17","q18","q20","q10","q12","q14","q19","q21","q25","q26","q27"), meanstructure = TRUE)
, чтобы получить результаты для всех вмененных наборов данных, я получаю сообщение об ошибке:
Ошибкаin eigen (VarCov, симметричный = TRUE, only.values = TRUE): бесконечные или отсутствующие значения в 'x'
Что я делаю неправильно?Отсутствуют пропущенные значения ни в одном из вмененных наборов данных, и когда я запускаю функцию sem для каждого из вмененных наборов данных в отдельности, она работает нормально (но я не получаю объединенный результат).
Все переменныетипа integer, за исключением пола (числовой) и возраста (числовой).И кроме пола, возраста и частоты все они содержат упорядоченные значения 1 2 3 4.
Кажется, что предиктор пола (закодированный как фиктивный 0, 1) вызвал некоторые проблемы.Но я не понимаю, почему, тогда, он прекрасно работает при вычислении одной модели, а не объединенной модели через функцию runMI.