Sottoinsiemi e modifica della periodicità
Il prossimo passo è unire i tuoi dati di temperatura con i dati sui voli del capitolo precedente.
Ricorda che, nel capitolo precedente, i dati sui voli coprono il periodo dal 2010 al 2015 con periodicità mensile. Al contrario, i tuoi dati di temperatura vanno dal 2007 al 2015 con periodicità giornaliera. Prima di unirli, dovresti creare un sottoinsieme dei dati e portare la periodicità a mensile.
Per convertire la periodicità degli oggetti xts, puoi usare to.period(), che ti permette di convertire rapidamente i dati a una frequenza inferiore. Per impostazione predefinita, questo comando produce valori specifici sull'intero periodo (ossia Open-High-Low-Close, o OHLC), utili nell'analisi finanziaria ma non sempre rilevanti.
In questo caso, imposta l'argomento OHLC su FALSE. Invece di generare colonne OHLC nel tuo oggetto xts mensile, questa impostazione prenderà semplicemente una riga per ciascun periodo, rappresentativa dell'intero mese. Puoi specificare quale riga usare con l'argomento indexAt.
Sia i dati temps_xts sia i dati flights_xts (dal capitolo precedente) sono disponibili nel tuo workspace.
Questo esercizio fa parte del corso
Caso di studio: analisi di serie temporali cittadine in R
Istruzioni dell'esercizio
- Crea un sottoinsieme dell'oggetto
temps_xtsincludendo solo le osservazioni dal 2010 al 2015. Salvalo cometemps_xts_2. - Usa
to.period()per convertire i dati giornalieri di temperatura a periodicità mensile. Assicurati di specificare il periodo di destinazione ("months"). Devi anche impostareOHLCsuFALSEper evitare di generare nuove colonne OHLC. Infine, imposta l'argomentoindexAtsu"firstof"per selezionare la prima osservazione di ogni mese. - Usa due chiamate a
periodicity()per confrontare la periodicità e la durata del nuovo insieme di dati mensili di temperatura con i datiflights_xtsdel capitolo precedente.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# 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(___)