殘差（）函數錯誤：替換有x行並且數據有y行

Question

我有一個數據集，它有多個人閱讀的每個單詞的閱讀時間。

我正在計算我的數據中每個人的閱讀時間殘差。單詞長度和表達順序（特定單詞的順序）是計算每個人的迴歸的因素。

讀取時間進行了對數轉換（logRT），字長由nchar（）計算。演示順序也是對數轉換的。

model1<-lmer(logRT~wlen+log(order)+(1|subject), data=mydata) 

然後，我嘗試通過執行以下操作以獲取每個數據點的剩餘列，

mydata$logResid<-residuals(model1) 

然後，我得到這個錯誤。

Error in `$<-.data.frame`(`*tmp*`, "LogResid", value = c(0.145113408056189, : 
replacement has 30509 rows, data has 30800 

有沒有人有任何建議？我完全困惑。由於這是我每天都在做的分析，迄今爲止沒有這樣的錯誤。它更令人困惑。

Answer 1

我會說你應該嘗試

model1 <- lmer(logRT~wlen+log(order)+(1|subject), data=mydata, 
       na.action=na.exclude) 

，並看看是否有幫助;它應該在適當的地方填寫NA值。

從?na.exclude：

... when ‘na.exclude’ is used the residuals and predictions are padded to the correct length by inserting ‘NA’s for cases omitted by ‘na.exclude’.