Utilisation de la prévision de prévision (auto.arima())

Question

J'ai récemment mis à jour vers la version 4.03() et la dernière ligne (ligne 4) de l'exemple ci-dessous donne maintenant un message d'erreur (comme indiqué en bas). Notez que forecast() dans la ligne 4 est alimenté par la sortie de auto.arima() à partir de la ligne 3 (qui fonctionne sans erreurs). Est-ce que quelque chose a changé dans le paquet de prévisions?

En outre, le message d'erreur disparaît lorsque le terme de zoo dans la ligne 3 est remplacé par un terme ts en utilisant le code suivant:

autarimod <- auto.arima(log(as.ts(zooinpdat))) ##New line 3 

Alors, est la combinaison prévision (auto.arima()) pas accepter plus longtemps les objets du zoo? Si c'est le cas, existe-t-il une meilleure façon de gérer cela que la méthode as.ts()?

library(zoo) 
library(forecast) 

inpdat <- c(353.03, 383.06, 407.9, 420.58, 345.96, 299.73, 286.42, 291.03, 
    297.71, 300.92, 272.13, 283.58, 331.72, 372.95, 404.78, 403.04, 
    374.57, 332.94, 284.37, 311.78, 307.27, 302.42, 283.52, 288.64, 
    337.19, 416.35, 418.65, 431.51, 407.74, 319.28, 297.33, 314.83, 
    290.49, 309.38, 294.5, 330.63, 371.2, 418.76, 440.05, 467.23, 
    384.32, 329.81, 300.4, 318.9, 355.06, 329.93, 293.43, 297.76, 
    340.42, 393.09, 395.2, 443.13, 396.45, 341.96, 307.95, 322, 339.63, 
    312.12, 304.31, 310.95) 

zooinpdat <- zooreg(inpdat, frequency=12, start=as.yearmon("May 1965")) 

autarimod <- auto.arima(log(zooinpdat)) ##Line 3 

for_arima <- forecast(autarimod, level=0.98, h=48) ##Line 4 


Error in .cbind.ts(list(e1, e2), c(deparse(substitute(e1))[1L], deparse(substitute(e2))[1L]), : 
not all series have the same frequency 

Answer 1

Le paquet prévu est destiné à ts objets, objets non zoo. Certaines fonctions fonctionneront bien avec les objets zoo, mais il n'y a aucune garantie. En particulier, quand je fais des changements au paquet, je ne vérifie jamais si les changements causeront des problèmes si zoo objets sont utilisés.

Vous pouvez corriger votre erreur en utilisant

zooinpdat <- as.ts(zooinpdat)