我尝试使用predict函数为具有splines::bs项的glmer模型获取模型拟合值。predict函数抛出以下错误：

Error in `[.data.frame`(fr, vars) : undefined columns selected

字符串
以下是使用模拟数据的可再现示例：

# Reproducible example of glmer spline predict problem

library(tibble)
library(dplyr)
library(splines)
library(lme4)

## Simulate data
set.seed(1234L)

n_groups <- 100L

num_per_group <- rpois(n_groups, 50L) + 1L

group_mean <- rnorm(n_groups, sd = 0.2)

x_vals <- c(-20:19)

simulated_data <- tibble(x = sample(x_vals,
                                    sum(num_per_group),
                                    replace = TRUE),
                         group_no = rep(c(1:n_groups),
                                        num_per_group),
                         group_mean = rep(group_mean,
                                          num_per_group)) %>%
  arrange(x, group_no) %>%
  mutate(lin_pred = 3 - 0.05*abs(x) + group_mean,
         y = exp(lin_pred + rnorm(sum(num_per_group),
                                  mean = 0, sd = 0.1)),
         y = round(y))

# Fit model

glmer_spline <- glmer(y ~ bs(x,
                             knots = c(0),
                             degree = 1L) +
                        (1 | group_no),
                      data = simulated_data,
                      family = poisson(link = "log"),
                      control = glmerControl(optimizer="bobyqa")
)

summary(glmer_spline)

# This throws an error:
# Error in `[.data.frame`(fr, vars) : undefined columns selected

prediction_output_splines <- predict(glmer_spline,
                                     newdata = simulated_data,
                                     type = "response",
                                     allow.new.levels = TRUE,
                                     re.form=~0)

型
我正在使用R版本4.2.3（2023-03-15），tibble 3.2.1，dplyr 1.1.2，lme 4 1.1-34，在RStudio“Mountain Hydrangea”Release（547 dcf 86，2023-07-06）中用于运行Ventura 13.4.1的2019 MacBook Pro上的macOS。
我想知道我的代码中是否有错误，我是否应该使用不同的函数来获得预测，我试图做的事情是否存在某种缺陷，或者是否在其中一个包中存在错误。任何帮助都非常感谢。