Filtrando e ajustando a periodicidade
Seu próximo passo é mesclar seus dados de temperatura com os dados de voos do capítulo anterior.
Relembre que, no capítulo anterior, seus dados de voos vão de 2010 a 2015 em períodos mensais. Já os dados de temperatura vão de 2007 a 2015 em períodos diários. Antes de mesclar, você deve filtrar seus dados e ajustar a periodicidade para mensal.
Para converter a periodicidade de objetos xts, você pode usar to.period(), que permite converter rapidamente seus dados para uma frequência mais baixa. Por padrão, esse comando produz valores específicos ao longo de todo o período (ou seja, Open-High-Low-Close, ou OHLC), úteis em análises financeiras, mas que podem não ser relevantes em todos os contextos.
Neste caso, você deve definir o argumento OHLC como FALSE. Em vez de gerar colunas OHLC no seu objeto xts mensal, essa configuração simplesmente seleciona uma linha de cada período como representativa de todo o período. Você pode especificar qual linha usando o argumento indexAt.
Tanto os dados temps_xts quanto os dados flights_xts (do capítulo anterior) estão disponíveis no seu workspace.
Este exercício faz parte do curso
Estudo de Caso: Analisando Séries Temporais de Cidades em R
Instruções do exercício
- Filtre seu objeto
temps_xtspara incluir apenas observações de 2010 a 2015. Salve comotemps_xts_2. - Use
to.period()para converter seus dados diários de temperatura para periodicidade mensal. Lembre-se de especificar o período para o qual deseja converter ("months"). Você também precisa definirOHLCcomoFALSEpara evitar a criação de novas colunas OHLC. Por fim, defina o argumentoindexAtcomo"firstof"para selecionar a primeira observação de cada mês. - Use duas chamadas de
periodicity()para comparar a periodicidade e a duração dos seus novos dados mensais de temperatura com os dadosflights_xtsdo capítulo anterior.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# 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(___)