【发布时间】:2016-11-23 00:27:57
【问题描述】:
我已经使用 Ravaan 包中的 growth() 函数成功地对一项研究的纵向数据进行了建模。我在任何地方都找不到关于如何为每个参与者提取预测轨迹的文档。我只能找到整个组的预测轨迹(在摘要输出的“拦截”部分下给出)。
【问题讨论】:
-
您能否发布一个您正在拟合的模型的可重现示例:我猜 lavaan 函数帮助文件会有一个。
标签: r factor-analysis r-lavaan
我已经使用 Ravaan 包中的 growth() 函数成功地对一项研究的纵向数据进行了建模。我在任何地方都找不到关于如何为每个参与者提取预测轨迹的文档。我只能找到整个组的预测轨迹(在摘要输出的“拦截”部分下给出)。
【问题讨论】:
标签: r factor-analysis 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, value = c(2.14028787487844, : subscript out of bounds) 中的错误”我也许应该提到两件事可能与此错误有关。一, Lavaan 在我的模型中给我一个警告信息,关于某些情况是空的,因此被删除,另外两个,而不是简单的截距和斜率项,我的一些模型具有更高阶项,表示随时间的二次和/或三次变化。跨度>