R中线性模型的自相关

Question

我正在尝试为我的BSc荣誉项目解决数据的自动关联问题，但这是我的第一次，我真的可以使用一些建议。

我试图在12年的时间内模拟月平均物种丰度和平均月海表温度（SST）;

    modelw<-lm(SPECIES~SST,
           data=dat, na.action=na.exclude)
summary(modelw)

我有残差和ACF

Modelw residuals

Modelw ACF

我可以看到存在自相关，所以我尝试使用以下方法进行修复;

> library(nlme)
    modelw2 <- gls(SPECIES ~SST, data=dat, 
                  correlation = corAR1(form=~SPECIES),
                  na.action=na.omit)
    summary(modelw2)

这没有用，因为（我认为）如果其中一个变量是日期，它应该被使用？有人可以指点一下我可以尝试的脚本的正确方向，它将与平均丰度和SST一起使用吗？

非常感谢！

Answer 1

我认为这是经过纠正的剧本;

library(nlme)
modelw2 <- gls(SPECIES ~SST, data=dat, 
               correlation = corAR1(form=~MONTH),
               na.action=na.omit)
summary(modelw2)

form应该是包含日期的列。