如何在 r 中对 lmList 对象执行 Breusch 异教测试？

我正在使用 lme4 包中的 lmList 对一个数据框进行多重回归。在每个回归上，我都想做一个异方差测试。尝试了 Breusch 异教测试但没有成功。这是示例数据：

library(lme4)
set.seed(123)
dat <- data.frame(group=c(rep("a",10), rep("b",10), rep("c",10)), 
                   x=rnorm(30), y1=rnorm(30), y2=rnorm(30), y3=rnorm(30))

regs <- lmList(cbind(y1, y2, y3) ~ x | group, dat)

导致九次回归。到目前为止我的尝试（使用“lmtest”包）是：

library(lmtest)

lapply(regs, bptest)

然而，这仅对 a、b 和 c 进行测试，而不是对这些中的每个 y1、y2 和 y3 进行测试。所以三个测试而不是九个。还尝试在 a 组内“到达”以对该组进行三个测试：

asum <- summary(regs$a)
bptest(asum$`Response y1`)

导致错误消息“model.matrix.default(terms(formula), model.frame(formula)) 中的错误：model.matrix() 中的模型框架和公式不匹配”​​。

关于如何使这项工作的任何想法？