CommencerCommencer gratuitement

Utiliser tsCV() pour la validation croisée de séries temporelles

La fonction tsCV() calcule les erreurs de validation croisée pour séries temporelles. Elle nécessite de spécifier la série, la méthode de prévision et l’horizon de prévision. Voici l’exemple utilisé dans la vidéo :

> e = tsCV(oil, forecastfunction = naive, h = 1)

Ici, vous allez utiliser tsCV() pour calculer et tracer les valeurs de MSE jusqu’à 8 pas à l’avance, avec la méthode naive() appliquée aux données goog. L’exercice utilise les graphiques ggplot2, que vous ne maîtrisez peut-être pas encore, mais nous avons fourni suffisamment de code pour que vous puissiez compléter le reste.

N’oubliez pas de vous référer aux diaporamas sur tsCV() dans le cours. Les données goog ont été chargées dans votre environnement de travail.

Cet exercice fait partie du cours

Prévision en R

Afficher le cours

Instructions

  • En utilisant les données goog et la fonction de prévision naive(), calculez les erreurs en validation croisée jusqu’à 8 pas à l’avance. Assignez le résultat à e.
  • Calculez les valeurs de MSE pour chaque horizon de prévision et éliminez les valeurs manquantes dans e en spécifiant le second argument. L’expression pour calculer le MSE a été fournie.
  • Tracez les valeurs de MSE résultantes (y) en fonction de l’horizon de prévision (x). Mobilisez vos connaissances sur les fonctions. Si MSE = mse est fourni dans la liste des arguments d’une fonction, alors mse doit faire référence à un objet qui existe dans votre espace de travail en dehors de la fonction, tandis que MSE est la variable qui désigne cet objet à l’intérieur de votre fonction.

Exercice interactif pratique

Essayez cet exercice en complétant cet exemple de code.

# Compute cross-validated errors for up to 8 steps ahead
e <- tsCV(___, forecastfunction = ___, h = ___)

# Compute the MSE values and remove missing values
mse <- colMeans(e^2, na.rm = ___)

# Plot the MSE values against the forecast horizon
data.frame(h = 1:8, MSE = mse) %>%
  ggplot(aes(x = h, y = ___)) + geom_point()
Modifier et exécuter le code