lmer模型残差的提取——处理缺失值

lmer模型残差的提取——处理缺失值,r,lme4,R,Lme4,我正试图从一个模型中提取残差,我已经使用lmer()估计了残差,并将其存储为一个新变量。有一些缺少的值,所以我尝试使用子集来跳过这些行,以便让残差在数据帧中的正确点结束 这就是我所做的: df$newvariable[!is.na(df$na)]<-resid(model) df$newvariable[!is.na(df$na)]要满足MCVE在R标签中发布此错误的要求,请包括在问题本身中完全重现此错误所需的所有数据和代码。话虽如此,我想您需要在训练模型之前输入缺失值。尝试在lmer(

我正试图从一个模型中提取残差,我已经使用lmer()估计了残差,并将其存储为一个新变量。有一些缺少的值,所以我尝试使用子集来跳过这些行,以便让残差在数据帧中的正确点结束

这就是我所做的:

df$newvariable[!is.na(df$na)]<-resid(model)

df$newvariable[!is.na(df$na)]要满足MCVE在R标签中发布此错误的要求,请包括在问题本身中完全重现此错误所需的所有数据和代码。话虽如此,我想您需要在训练模型之前输入缺失值。尝试在
lmer()
函数中设置
na.action=na.exclude
,然后使用
resist
提取残差。对于模型中一个或多个变量缺少数据的每个观测值,得出的残差长度应等于NA值(而不是像
NA.omit
那样忽略这些观测值)。为了满足MCVE在R标签中发布此错误的要求,请在问题本身中包含完全重现此错误所需的所有数据和代码。话虽如此,我想您需要在训练模型之前输入缺失值。尝试在
lmer()
函数中设置
na.action=na.exclude
,然后使用
resist
提取残差。对于模型中一个或多个变量缺少数据的每个观测值,得出的残差长度应等于NA值(而不是像
NA.ommit
那样忽略这些观测值)。
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