Subsetare și ajustarea periodicității
Următorul pas este să combini datele despre temperatură cu datele despre zboruri din capitolul anterior.
Reține din capitolul anterior că datele despre zboruri acoperă perioada 2010–2015, cu periodicitate lunară. Prin contrast, datele despre temperatură acoperă perioada 2007–2015, cu periodicitate zilnică. Înainte de a le combina, ar trebui să faci o subsetare a datelor și să ajustezi periodicitatea la lunar.
Pentru a converti periodicitatea obiectelor xts, poți folosi to.period(), care îți permite să treci rapid la o frecvență mai mică. Implicit, această funcție produce valori specifice pentru întreaga perioadă (și anume, Open-High-Low-Close sau OHLC), utile în analiza financiară, dar care nu sunt relevante în orice context.
În cazul de față, setează argumentul OHLC la FALSE. În loc să genereze coloane OHLC în obiectul xts lunar, această setare va prelua pur și simplu un rând din fiecare perioadă ca reprezentativ pentru întreaga perioadă. Poți specifica ce rând să fie preluat folosind argumentul indexAt.
Atât datele temps_xts, cât și datele flights_xts (din capitolul anterior) sunt disponibile în spațiul tău de lucru.
Acest exercițiu face parte din cursul
Studiu de caz: Analiza seriilor de timp pentru un oraș în R
Instrucțiuni pentru exercițiu
- Fă un subset al obiectului
temps_xtspentru a include doar observațiile din perioada 2010–2015. Salvează rezultatul catemps_xts_2. - Folosește
to.period()pentru a converti datele zilnice despre temperatură la periodicitate lunară. Asigură-te că specifici perioada dorită ("months"). De asemenea, seteazăOHLClaFALSEpentru a evita generarea de coloane OHLC noi. În final, setează argumentulindexAtla"firstof"pentru a selecta prima observație din fiecare lună. - Folosește două apeluri ale funcției
periodicity()pentru a compara periodicitatea și durata noilor date lunare despre temperatură cu dateleflights_xtsdin capitolul anterior.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
# 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(___)