2

R の Lavaan パッケージの growth() 関数を使用して、調査から縦方向のデータをうまくモデル化しました。グループ全体の予測軌道のみを見つけることができます (要約出力の「Intercepts」部分の下に表示されます)。

4

1 に答える 1

2

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
于 2016-07-19T23:33:36.937 に答える