#r #forecast
Вопрос:
Привет, я хочу вычислить STLM
модель из forecast
пакета, используя elm
функцию из nnfor
пакета для серии с учетом сезонных колебаний.
Я могу сделать это, используя modelfunction
аргумент then.
Давайте воспользуемся этим временем.:
y<- ts(rnorm(300,mean=0,sd=1), start=c(2012,2), frequency=12)
> head(y,10)
Feb Mar Apr May Jun
2012 1.83672399 0.09163453 1.43897758 -0.06213232 0.21279124
Jul Aug Sep Oct Nov
2012 1.43642897 1.36854342 -0.38119914 -0.94666416 1.41982674
Он отлично работает с функциями из forecast
пакета (например: nnetar
)
model1 =stlm(y,modelfunction=forecast::nnetar)
model$model
Series: x.sa
Model: NNAR(24,1,12)[12]
Call: modelfunction(y = x.sa)
Average of 20 networks, each of which is
a 24-12-1 network with 313 weights
options were - linear output units
sigma^2 estimated as 0.0004526
Но не из моделей прогнозирования из других пакетов.
model2 =stlm(y,modelfunction=nnfor::elm)
> Error in ` .default`(fitted(fit), allseas) :
time-series/vector length mismatch
Какой формат мы должны использовать для аргумента modelfunction ?
Есть какие-нибудь идеи, чтобы помочь мне ? Большое спасибо.