我正在学习如何使用 gls 进行预测，但是当我尝试这样做时，我不断收到错误消息。这是我正在使用的代码。
data("airquality")
airquality$Date=as.Date(paste(1973,airquality$Month,airquality$Day,sep="-"))
air1<-subset(airquality,complete.cases(Ozone))
library(nlme)
model1=gls(Ozone~Wind*Date,air1,correlation=corAR1())
到目前为止，这工作正常。现在，我想预测特定日期特定风速的臭氧。所以我用。
newdata<-data.frame(Wind=10,Date="1973-05-04")
newdata$Date<-as.Date(newdata$Date)
predict(model1,newdata)
最后一行代码给了我以下错误：
eval 中的错误（expr、envir、enclos）：找不到对象“臭氧”
然后我尝试了以下代码来预测：
predict.gls(model1,newdata)
这是我得到的错误，尽管该nlme库已包含在内。
错误：找不到函数“predict.gls”
谁能告诉为什么会出现这些错误？另外，在这里使用什么代码进行预测的正确方法是什么？