У меня есть фрейм данных (длиной 3575 наблюдений) с двумя столбцами: дата и концентрация кислорода.Дата указана в формате ГГГГ-ММ-ДД в течение 25 лет, а значениями кислорода являются числа в диапазоне от 0,05 до 7,5.Я пытаюсь найти выбросы в моем наборе данных с помощью пакета tsoutliers, и моя проблема заключается в синтаксисе, а не в значениях самих данных.
Впервые я использовал функцию auto.arima с этим кодом:
auto.arima(Site7$Oxygen, seasonal = TRUE, stepwise = FALSE)
Мне очень жаль, что я не смог привести минимальный воспроизводимый пример, но с учетом масштаба данных и тестовбыло бы довольно сложно поместить все это здесь.
Для результатов, которые я получил:
Series: Date$Oxygen
ARIMA(3,1,2)
Coefficients:
ar1 ar2 ar3 ma1 ma2
0.2988 0.2439 0.1431 -0.4348 -0.5259
s.e. 0.0947 0.0686 0.0316 0.0948 0.0772
sigma^2 estimated as 0.5999: log likelihood=-4145.53
AIC=8303.06 AICc=8303.09 BIC=8340.15
Теперь, когда я пытаюсь визуализировать аддитивные выбросы, я использовал tso () команда с этим кодом:
AOoutliers <- tso(Date$Oxygen,types = "AO", maxit = 10, tsmethod = c("arima"), args.tsmethod = list(order =c(3,1,2)))
Теперь у меня есть эта ошибка, которую я не могу понять:
Ошибка в tso0 (x = y, xreg = xreg, cval = cval, delta = delta, n.start = n.start: попытка получить слот "y" из объекта базового класса ("числовой") без слотов
AnyПомощь будет принята с благодарностью. Спасибо.