[問題] augPred function 出現錯誤訊息?

作者: f12508 (心寬念存)   2024-07-27 09:51:05
請教版上各位高手
我試著用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
請問該如何調整呢
謝謝大家
作者: Wush978 (拒看低質媒體)   2024-07-28 08:31:00
看起來是資料的問題,是嗎?
作者: f12508 (心寬念存)   2024-07-28 19:50:00
但資料跑lme 是沒問題的
作者: locka (locka)   2024-07-29 12:22:00
請問能提供測試資料嗎?

Links booklink

Contact Us: admin [ a t ] ucptt.com