Я успешно смоделировал продольные данные из исследования с использованием функции роста() в пакете Lavaan для R. Я не могу найти его где-либо документально зафиксированным, как извлечь предсказанные траектории для каждого участника. Я могу найти только предсказанную траекторию для группы в целом (данный в разделе «Перехватывает» итоговый вывод).Извлечение отдельных траекторий в модели кривой роста Lavaan
Использование lavPredict
model.syntax <- '
# intercept and slope with fixed coefficients
i =~ 1*t1 + 1*t2 + 1*t3 + 1*t4
s =~ 0*t1 + 1*t2 + 2*t3 + 3*t4
# regressions
i ~ x1 + x2
s ~ x1 + x2
# time-varying covariates
t1 ~ c1
t2 ~ c2
t3 ~ c3
t4 ~ c4
'
fit <- growth(model.syntax, data=Demo.growth)
summary(fit)
head(lavPredict(fit))
Который будет производить предсказанные оценки для каждого отдельного
i s
[1,] 1.1378809 0.676301228
[2,] -2.5421940 -1.425974525
[3,] -0.1279434 0.966734762
[4,] 1.1682777 1.477200679
[5,] -0.5141435 0.006995809
[6,] -1.2646865 0.524024630
Спасибо, я считаю, что это то, что я ищу. Однако я попробовал это со многими моделями, и все они дают аналогичную ошибку. "Ошибка в' [<-' ('* TMP *',, ov.x.dummy.lv.idx, значение = с (2.14028787487844,: подстрочные вне границ" я, возможно, следует упомянуть две вещи, может быть связано с этой ошибкой. Во-первых, Lavaan дает мне предупреждающее сообщение в моих моделях о том, что некоторые случаи пустые и поэтому удалены, а два, а не простой термин перехвата и склона, некоторые из моих моделей имеют более высокие порядковые члены, представляющие квадратичные а также кубическое изменение с течением времени. –
Я не использовал lavPredict в производстве, однако, играя с ним, я знаю, что ему не нравятся отсутствующие данные. Я думаю, что для такого рода вещей лучше всего задать вопрос о лаване Страница групп google. –
Можете ли вы опубликовать воспроизводимый пример модели вы фитинга пожалуйста: Id думаю, функция lavaan файлы справки будут иметь один. – user20650