我正在调查环境/海洋因素如何影响渔业上岸量。我在 R 中有以下 GAMM，使用该gam函数：
GAMM1 <- gam(log(weight.landed + 1) ~ s(day.in.series) + s(month) + s(wind.speed) + s(wind.direction,bs="cc") + s(SST) + s(Chl.a) + s(fisherID,bs="re") + offset(log(effort)), data=mydata)
fisherID是识别个体渔民的变量，是 GAMM 的随机部分。
effort用作抵消，以使渔获物的重量标准化为使用的渔具数量。
day.in.series是时间序列 (1 - 6000) 中的有序天数，month是显示季节性影响的日历月。所有其他变量都是各种环境数据。
ACF:
和部分 ACF:
图表明数据中存在自相关。在 ACF 图中，在延迟 30 左右之前有明显的滞后，我不完全确定部分 ACF 图告诉我什么。
根据我的研究，我需要将 corARMA 合并到我的模型中。可以使用该gam功能来完成，还是gamm需要使用该功能？我也不确定随机变量 ( fisherID) 和偏移量是否会以任何方式影响这一点。