Delmängder och justering av periodicitet
Nästa steg är att slå samman temperaturdata med flygdata från föregående kapitel.
Kom ihåg att flygdata sträcker sig från 2010 till 2015 i månatliga perioder. Temperaturdata täcker däremot perioden 2007 till 2015 i dagliga perioder. Innan du slår samman dataseten behöver du skapa en delmängd av dina data och justera periodiciteten till månadsvis.
För att konvertera periodiciteten hos xts-objekt kan du använda to.period(), som gör det enkelt att omvandla dina data till en lägre frekvens. Som standard genererar kommandot specifika värden för hela perioden – nämligen Open-High-Low-Close (OHLC) – vilket är användbart inom finansiell analys men inte alltid relevant i andra sammanhang.
I det här fallet bör du sätta argumentet OHLC till FALSE. I stället för att generera OHLC-kolumner i ditt månatliga xts-objekt väljer den inställningen helt enkelt en rad per period som representativ för hela perioden. Du kan ange vilken rad med argumentet indexAt.
Både temps_xts-data och flights_xts-data (från föregående kapitel) finns tillgängliga i din workspace.
Den här övningen är en del av kursen
Fallstudie: Analys av tidsseridata för städer i R
Övningsinstruktioner
- Skapa en delmängd av ditt
temps_xts-objekt som bara innehåller observationer från 2010 till 2015. Spara den somtemps_xts_2. - Använd
to.period()för att konvertera din dagliga temperaturdata till månadsperiodicitet. Se till att ange vilken period du vill konvertera till ("months"). Du behöver också sättaOHLCtillFALSEför att undvika att nya OHLC-kolumner genereras. Sätt slutligen argumentetindexAttill"firstof"för att välja den första observationen varje månad. - Använd
periodicity()i två anrop för att jämföra periodicitet och varaktighet hos din nya månatliga temperaturdata medflights_xts-data från föregående kapitel.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# Subset your temperature data to include only 2010 through 2015: temps_xts_2
temps_xts_2 <- ___["___/___"]
# Use to.period to convert temps_xts_2 to monthly periodicity
temps_monthly <- to.period(___, period = "___", OHLC = ___, indexAt = "___")
# Compare the periodicity and duration of temps_monthly and flights_xts
periodicity(___)
periodicity(___)