Výběr podmnožiny dat a úprava periodicity

Dalším krokem je sloučení dat o teplotách s daty o letech z předchozí kapitoly.

Z předchozí kapitoly víš, že data o letech pokrývají období 2010–2015 v měsíčních intervalech. Data o teplotách naproti tomu sahají od roku 2007 do roku 2015 v denních intervalech. Před sloučením je proto potřeba vybrat správné období a upravit periodicitu na měsíční.

Ke změně periodicity xts objektů slouží funkce to.period(), která umožňuje rychlý převod dat na nižší frekvenci. Ve výchozím nastavení funkce vypočítá pro každé období specifické hodnoty (konkrétně Open-High-Low-Close, neboli OHLC), které jsou užitečné ve finanční analýze, ale ne vždy relevantní v jiných kontextech.

V tomto případě nastav argument OHLC na FALSE. Místo generování OHLC sloupců v měsíčním xts objektu toto nastavení jednoduše vybere jeden řádek z každého období jako jeho reprezentativní hodnotu. Pomocí argumentu indexAt pak určíš, který řádek se použije.

V pracovním prostředí máš k dispozici data temps_xts i flights_xts z předchozí kapitoly.

Toto cvičení je součástí kurzu

Případová studie: Analýza časových řad městských dat v R

Zobrazit kurz

Pokyny k cvičení

Vyber z objektu temps_xts pouze pozorování z let 2010 až 2015. Výsledek ulož jako temps_xts_2.
Pomocí to.period() převeď denní data o teplotách na měsíční periodicitu. Nezapomeň zadat požadované období ("months"), nastavit OHLC na FALSE, aby se negenerovaly nové OHLC sloupce, a argument indexAt nastavit na "firstof", aby se z každého měsíce vybralo první pozorování.
Dvěma voláními funkce periodicity() porovnej periodicitu a rozsah nových měsíčních dat o teplotách s daty flights_xts z předchozí kapitoly.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

# 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(___)

Upravit a spustit kód