Utiliser merge() et tracer dans le temps
Maintenant que vous avez des données de température qui couvrent la même période (2010-2015) et la même fréquence (mensuelle) que vos données de vols, vous êtes prêt à faire la fusion.
Pour fusionner des objets xts par colonne, vous pouvez utiliser merge(). Lorsque deux objets xts partagent la même périodicité, merge() réussit généralement à combiner l'information dans les bonnes lignes. Même lorsque des objets xts n'ont pas la même périodicité, merge() conservera le bon ordre chronologique de ces objets à travers des périodes différentes.
Dans cet exercice, vous fusionnerez vos deux objets xts par colonne et vous générerez de nouveaux graphiques pour explorer le lien entre les retards de vols et la température. temps_monthly et flights_xts sont disponibles dans votre espace de travail.
Cette activité fait partie du cours
Étude de cas : analyser des séries chronologiques urbaines en R
Instructions de l’exercice
- Utilisez
merge()pour combinerflights_xtsettemps_monthly. Comme ces objets xts partagent la même périodicité, votre commande de fusion devrait placer les données de température dans la bonne ligne de votre objetflights_xts. Notez que l'ordre dans lequel vous indiquez les objets à fusionner détermine l'emplacement des colonnes dans l'objet fusionné. Pour rester cohérent, mettezflights_xtsen premier ettemps_monthlyen second. - Examinez les premières lignes de votre objet xts fusionné (
flights_temps) pour confirmer que la fusion a réussi. Vous devriez voir les données de température alignées avec les données de vols. - Utilisez
plot.zoo()pour générer un seul graphique contenant les colonnespct_delayettemps_monthlydeflights_temps. Assurez‑vous de sous‑sélectionner les colonnes pertinentes et de définirplot.typeà"single". Laissez l'argumentltytel quel.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# Use merge to combine your flights and temperature objects
flights_temps <- merge(___, ___)
# Examine the first few rows of your combined xts object
# Use plot.zoo to plot these two columns in a single panel
plot.zoo(___[,c("___", "___")], plot.type = "___", lty = lty)
legend("topright", lty = lty, legend = labels, bg = "white")