Функция модели в STLM

#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 ?

Есть какие-нибудь идеи, чтобы помочь мне ? Большое спасибо.