r - na.fail.default 中的错误 : missing values in object - but no missing values

Question

我正在尝试使用这些数据运行 lme 模型:

tot_nochc=runif(10,1,15)
cor_partner=factor(c(1,1,0,1,0,0,0,0,1,0))
age=runif(10,18,75)
agecu=age^3
day=factor(c(1,2,2,3,3,NA,NA,4,4,4))
dt=as.data.frame(cbind(tot_nochc,cor_partner,agecu,day))
attach(dt)

corpart.lme.1=lme(tot_nochc~cor_partner+agecu+cor_partner *agecu, 
                  random = ~cor_partner+agecu+cor_partner *agecu |day, 
                  na.exclude(day))

我收到此错误代码:

Error in na.fail.default(list(cor_partner = c(1L, 1L, 2L, 1L, 1L, 1L, : missing values in object

我知道论坛中有类似的问题。然而，就我而言:

• cor_partner 没有缺失值；
• 整个对象被编码为一个因素(至少从全局环境中显示)。

我可以使用 na.action 排除这些 NA 值，但我更想知道为什么该函数会读取缺失值 - 以准确了解我的数据发生了什么。

Answer 1

tl;dr 您必须立即在整个数据帧上使用 na.exclude() (或其他)，以便剩余的观察结果保持匹配变量...

set.seed(101)
tot_nochc=runif(10,1,15)
cor_partner=factor(c(1,1,0,1,0,0,0,0,1,0))
age=runif(10,18,75)
agecu=age^3
day=factor(c(1,2,2,3,3,NA,NA,4,4,4))
## use data.frame() -- *DON'T* cbind() first
dt=data.frame(tot_nochc,cor_partner,agecu,day)
## DON'T attach(dt) ...

现在尝试:

library(nlme)
corpart.lme.1=lme(tot_nochc~cor_partner+agecu+cor_partner *agecu, 
              random = ~cor_partner+agecu+cor_partner *agecu |day, 
              data=dt,
              na.action=na.exclude)

我们收到收敛错误和警告，但我认为这是因为我们使用的是一个微小的虚构数据集，其中没有足够的信息，而不是因为代码存在任何固有问题。