Извлечь остатки из цикла lme4

Question

Я хочу извлечь остатки из линейного смешанного модельного цикла.

results <- matrix(NA, nrow(meth),4,   dimnames=list(rownames(meth),  c("nobs", "coef", "coef.se", "Pval")))

Цикл создает матрицу с результатами для переменной X, но теперь я хочу извлечь из этого цикла остатки, которые бы создали ту же матрицу, что и мет (те же измерения), но затем с остатками. Я старался:

#make matrix
res<- matrix(NA, nrow(meth),ncol(meth), dimnames=list(rownames(meth), colnames(meth)) )
#run loop
for (i in 1:(nrow(meth))) {
  y <- (meth[i,])
  model <- try(lmer(as.formula(f), data=(sample_trans))
               , silent=TRUE)
  if (!inherits(model, "try-error")) {
      res[i,] <- resid(model)
  }
}

но я получаю ошибку:

Error in res[i, ] <- resid(model) : 
number of items to replace is not a multiple of replacement length

В данных есть NA, но это не должно влиять на это, потому что я использовал try. Что-то не так с самой петлей?