CommencezCommencez gratuitement

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

La fonction tsCV() calcule les erreurs de validation croisée pour séries temporelles. Elle exige de préciser 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 utiliserez tsCV() pour calculer et tracer les valeurs de MSE jusqu'à 8 pas d'avance, en utilisant la méthode naive() appliquée aux données goog. L'exercice utilise les graphiques de ggplot2, que vous ne connaissez peut-être pas, mais nous avons fourni suffisamment de code pour que vous puissiez compléter le reste.

Assurez-vous de consulter les diapositives sur tsCV() dans le cours. Les données goog ont été chargées dans votre espace de travail.

Cette activité fait partie du cours

Prévision en R

Voir le cours

Instructions de l’exercice

  • En utilisant les données goog et la fonction de prévision naive(), calculez les erreurs en validation croisée jusqu'à 8 pas d'avance. Assignez le résultat à e.
  • Calculez les valeurs de MSE pour chaque horizon de prévision et retirez les valeurs manquantes dans e en précisant le deuxième 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). Réfléchissez à votre connaissance des fonctions. Si MSE = mse est fourni dans la liste des arguments d'une fonction, alors mse doit renvoyer à un objet qui existe dans votre espace de travail à l'extérieur de la fonction, tandis que MSE est la variable qui fait référence à cet objet à l'intérieur de votre fonction.

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

# 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