r语言 - 带有嵌套变量的 lme 函数的问题：(any(notIntX <- ！apply(X， 2， const))

Question

我试图拟合一个混合效应模型，其中我有嵌套的协变量(VarX5|VarX6(，并被视为固定效应。

然而，我正在尝试调整数据，出现以下消息：

library(nlme)
library(lme4)
dados$VarCat=as.factor(dados$VarCat)
dados$VarX5=as.factor(dados$VarX5)
dados$VarX6=as.factor(dados$VarX6)
model <- lme(log(Resp)~log(VarX1)+log(VarX2)+(VarX3)+(VarX4)+VarX5|VarX6 ,random = ~1|VarCat, 
dados, method="REML")
Error in if (any(notIntX <- !apply(X, 2, const))) { : 
valor ausente onde TRUE/FALSE necessário

Answer 1

您没有正确指定分类变量。在nlme:lme中，它们进入random选项(错误是试图将因子直接传递到公式中(。在lme4: lmer中，它们被括在公式内的括号中。

model <- nlme::lme(
log(Resp) ~ log(VarX1) + log(VarX2) + VarX3 + VarX4,
random = ~ VarCat | VarX5:VarX6,
data = dados)
model <- lme4::lmer(
log(Resp) ~ log(VarX1) + log(VarX2) + VarX3 + VarX4 + (VarCat | VarX5:VarX6),
data = dados)

我不知道你对这个模型想要什么，所以:运算符可能不是最合适的。在lmer帮助页面中查看formula参数的详细信息。