2016-07-19 28 views
1

我已成功使用Lavaan軟件包中的grow()函數對R進行研究,從中成功建立了縱向數據的模型。我無法在任何地方找到它記錄如何爲每個參與者提取預測軌跡。我只能找到整個組的預測軌跡(根據摘要輸出的「截距」部分給出)。提取Lavaan增長曲線模型中的單個軌跡

+1

你可以發佈你的擬合模型的可重複的例子,請:身份證猜lavaan功能幫助文件將有一。 – user20650

回答

0

使用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 
+0

謝謝,我相信這是我正在尋找的。然而,我嘗試了很多模型,所有這些都給出了類似的錯誤。 「[<-'('* tmp *',,ov.x.dummy.lv.idx,value = c(2.14028787487844,: 下標越界)」錯誤「 我應該提到兩件事可能與這個錯誤有關:一,Lavaan在我的模型中給我一個關於某些情況是空的並因此被刪除的警告信息,以及兩個,而不是有一個簡單的截距和斜率項,我的一些模型的高階項代表二次方或立方隨時間變化 –

+0

我在製作中沒有使用過lavPredict,但是,在玩它的時候,我知道它不喜歡缺少數據,我認爲對於這樣的事情,你最好打賭會問關於lavaan的問題谷歌組合頁面。 –