請教版上各位高手
我試著用R跑Mixed-Effects Model
資料有108個受試者,每個受試者測了8次,總共有864筆資料
使用Mixed-Effects Model ,套件是nlme
M1 <- lme(V16 ~ 1 + V5 + V3 + V4 + V3*V4, random = ~1 | V11, data = data)
但要使用augPred 函數
aug.Pred1 <- augPred(M1,primary = ~V4, level = 0:1, length.out = 2 )
卻一值出現錯誤訊息
Error in [[<-.data.frame(*tmp*, nm, value = c(1L, 1L, 1L, 1L, 1L, :
替換資料裡有 864 列,但資料有 108
請問該如何調整呢
謝謝大家