在使用mgcv包运行gam模型时，我遇到了一个奇怪的错误消息，我无法理解：
“model.frame.default（formula = death ~ pm10 + Lag（resid1，1）+：变量长度不同（在“Lag（resid 1，1）”中找到）"。
模型1中使用的观测数与偏差残差的长度完全相同，因此我认为这个误差与数据大小或长度的差异无关。
我在网页here上发现了一个相当相关的错误信息，但该帖子没有得到足够的答案，所以它对我的问题没有帮助。
可再现的示例和数据如下：

library(quantmod)
library(mgcv) 
require(dlnm)

df <- chicagoNMMAPS
df1 <- df[,c("date","dow","death","temp","pm10")] 
df1$trend<-seq(dim(df1)[1]) ### Create a time trend

运行模型

model1<-gam(death ~ pm10 + s(trend,k=14*7)+ s(temp,k=5),
data=df1, na.action=na.omit, family=poisson)

获取偏差残差

resid1 <- residuals(model1,type="deviance")

在模型1中添加一天滞后偏差

model1_1 <- update(model1,.~.+ Lag(resid1,1),  na.action=na.omit)

model1_2<-gam(death ~ pm10 + s(trend,k=14*7)+ s(temp,k=5) + Lag(resid1,1), data=df1, 
na.action=na.omit, family=poisson)

这两个模型都产生了相同的错误消息。