我正在尝试从我估计的模型中提取残差，使用lmer（）来存储为新变量。有一些缺失值，所以我试图使用子集来跳过这些行，以便残差最终在数据帧中的正确点。
这就是我所做的：
df$newvariable[!is.na(df$na)]<-resid(model)
使用lmer函数估计模型。
我得到的错误是：
Error in `$<-.data.frame`(`*tmp*`, "newvariable", value = numeric(0)) : replacement has 0 rows, data has 3332
In addition: Warning message:
In is.na(df$na) : is.na() applied to non-(list or vector) of type 'NULL'
我尝试过其他一些事情，包括：
df$newvariable<-resid(model,na.action=na.omit)
这至少产生数据，但我知道它不完整，所以不适合它：
Error in `$<-.data.frame`(`*tmp*`, "e.bi", value = c(6.75180624992464, :
replacement has 3083 rows, data has 3332
我对R很新，对我正在做的事情不太了解。作为参考，我实际上试图复制在这篇期刊文章http://journal.frontiersin.org/article/10.3389/fpsyg.2016.00891/full中完成的分析，其中包括代码 - 请参阅第5页的特定部分。我不确定为什么我的不行！任何帮助将不胜感激。