我正在尝试执行具有随机效应的逐步模型，其中我可以获得 BIC 值。

lmerTest 包说它适用于 lme4，但只有从模型中删除我的一个自变量（这是一个具有两个选项 (TM) 的因素），我才能让它工作

错误代码是：

在错误$<-（*tmp*，式中，值=条款）：用于分配此S4类的子集的方法没有

或者

as_lmerModLmerTest(model) 中的错误：模型不属于“lmerMod”类：无法强制转换为“lmerModLmerTest”类

我在某处读过它可能与 drop1 有关，但我仍然没有弄清楚。我也愿意接受其他包和功能的建议。

之前，在尝试 full.model <- lm ( ... 一切正常。更改为 lmer 后，它不再起作用了。

我现在使用的代码：

full.model <- lme4::lmer(dep ~ TM + ind + (1 | dorp),  data=test)  #lmerTest:: give same outcome

step.model<- lmerTest::step(full.model, direction="both",k=log(16))   # n=16

summary(step.model)

BIC(step.model)

#Example dataset

test <- data.frame(TM = as.factor(c(rep("org", 3), rep("min", 3),rep("org", 3), rep("min", 3),rep("org", 3), rep("min", 3))),
                   dep = runif(18,0,20),
                   ind = runif(18,0,7),
                   dorp = as.factor(c(rep(1,6),rep(2,6),rep(3,6))))