Несколько общих замечаний:
- Вы можете установить модели SARIMA с функцией R *
arima
; на самом деле нет необходимости загружать дополнительный пакет для установки моделей Basi c SARIMA. - Набор данных
AirPassengers
является частью базы R, а не частью astsa
.
В arima
несезонные и сезонные параметры указываются через аргументы функций order
и seasonal
соответственно. Вы можете найти подробную информацию о том, какие аргументы принимает функция, если вы наберете ?arima
в терминале R.
Например, чтобы указать SARIMA (0,1,1) (0,1,1) 12 модель, которую вы делаете
arima(log(AirPassengers), order = c(0,1,1), seasonal = list(order = c(0,1,1), period = 12))
Обратите внимание, что аргумент seasonal
представляет собой list
с двумя элементами order
(обозначает порядок сезонных AR, порядок сезонных разностей, и порядок сезонных терминов MA) и period
.
Ошибка, с которой вы сталкиваетесь с astsa::sarima
, кажется, вызвана слишком большим количеством аргументов; убедитесь, что у вас действительно есть 3 несезонных аргумента, 3 (необязательных) сезонных аргумента и один (необязательный) период.
Над моделью будет
sarima(log(AirPassengers), 0, 1, 1, 0, 1, 1, 12)
, которая дает такие же результаты подбора как arima
(на самом деле sarima
просто звонит arima
"под капот").