用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)

最后一行代码出现以下错误：

Error in eval(expr, envir, enclos) : object 'Ozone' not found

然后我尝试使用以下代码进行预测：

predict.gls(model1,newdata)

这是我遇到的错误，尽管 nlme 库已包含在内。

Error: could not find function "predict.gls"

谁能告诉我为什么会出现这些错误？另外，在这里使用代码进行预测的正确方法是什么？

Answer 1

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

用gls预测

Predict with gls

r

predict

nlme