forecast包调用报错：参数'M'缺失无默认值，求解决方案

问题重现 #

用户运行以下R代码构建ARIMA和ETS模型并执行预测：

ARIMA_Modelb <- Arima(bankura_ts, order = c(0,0,1), 
                      seasonal = c(2,0,2))
ETS_Modelb <- ets(bankura_ts, damped = TRUE, model = "AZN")

# ARIMA模型预测绘图
autoplot(forecast(ARIMA_Modelb, h=72)) + theme_bw() + 
  ylab("Rainfall") + xlab("Datetime") + 
  scale_x_date(date_labels = '%b - %Y', 
               breaks = '5 year', minor_breaks = '6 month') +
  ggtitle("Bankura Rainfall Forecast 2021-2026 
  ARIMA Model")
# ETS模型预测绘图
autoplot(forecast(ETS_Modelb, h=72)) + theme_bw() + 
  ylab("Rainfall") + xlab("Datetime") + 
  scale_x_date(date_labels = '%b - %Y', breaks = '5 year', 
               minor_breaks = '6 month') +
  ggtitle("Bankura Rainfall Forecast 2021-2026 
  ETS Model")

执行后触发报错：

Error in forecast(fit1, h = 12) :  argument "M" is missing, with no default

用户已尝试重新安装forecast包，但问题仍存在，且其他调用forecast函数的代码也出现相同错误。

排查与解决步骤 #

1. 解决函数命名冲突 #

报错中的fit1并非用户代码中的对象，说明存在其他包的forecast函数覆盖了forecast包的同名函数。先确认当前forecast函数的来源：

find("forecast")

如果输出包含非package:forecast的结果，需明确指定命名空间调用forecast包的函数，修改预测代码：

# 用forecast::指定命名空间
autoplot(forecast::forecast(ARIMA_Modelb, h=72)) + theme_bw() + 
  ylab("Rainfall") + xlab("Datetime") + 
  scale_x_date(date_labels = '%b - %Y', 
               breaks = '5 year', minor_breaks = '6 month') +
  ggtitle("Bankura Rainfall Forecast 2021-2026 ARIMA Model")

autoplot(forecast::forecast(ETS_Modelb, h=72)) + theme_bw() + 
  ylab("Rainfall") + xlab("Datetime") + 
  scale_x_date(date_labels = '%b - %Y', breaks = '5 year', 
               minor_breaks = '6 month') +
  ggtitle("Bankura Rainfall Forecast 2021-2026 ETS Model")

2. 更新forecast包与依赖 #

旧版本forecast可能存在兼容性问题，执行以下代码安装最新版本并检查依赖：

install.packages("forecast", dependencies = TRUE)
library(forecast)
packageVersion("forecast")

确保版本为最新稳定版（截至2024年推荐v8.23及以上），若依赖包缺失，根据提示补装。

3. 验证时间序列对象格式 #

确认bankura_ts是合法的时间序列对象：

class(bankura_ts)
frequency(bankura_ts)

若对象类型不是ts/msts或频率信息缺失，重新转换为标准时间序列：

# 示例：假设数据框包含日期和降雨量列，月度数据设置frequency=12
bankura_ts <- ts(bankura_ts$Rainfall, frequency = 12)

4. 重启R会话 #

包加载冲突可能在会话中残留，重启R后重新加载相关包再运行代码：

library(forecast)
library(ggplot2)
# 重新构建模型并执行预测
ARIMA_Modelb <- Arima(bankura_ts, order = c(0,0,1), seasonal = c(2,0,2))
autoplot(forecast::forecast(ARIMA_Modelb, h=72)) + theme_bw() + ... # 后续绘图代码不变

内容的提问来源于stack exchange，提问作者Vivek sreekumar menon