用gls预测

Question

我正在学习如何使用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”

有人可以建议为什么会出现这些错误吗？另外，在这里使用什么方法进行预测的正确方法呢？

Answer 1

在进一步研究该主题后，我意识到predict.gls本身不是一个函数，而是在帮助文件中用于识别和解释将gls对象传递给预测函数所需的参数。