Filtrado y ajuste de periodicidad
Tu siguiente paso es combinar los datos de temperatura con los datos de vuelos del capítulo anterior.
Recuerda que tus datos de vuelos abarcan de 2010 a 2015 en periodos mensuales. En cambio, tus datos de temperatura cubren de 2007 a 2015 en periodos diarios. Antes de fusionarlos, deberías filtrar los datos y ajustar la periodicidad a mensual.
Para convertir la periodicidad de objetos xts, puedes usar to.period(), que te permite convertir rápidamente tus datos a una frecuencia inferior. De forma predeterminada, este comando produce valores específicos a lo largo de todo el periodo (concretamente, Open-High-Low-Close u OHLC), útiles en análisis financiero pero no siempre relevantes.
En este caso, debes establecer el argumento OHLC en FALSE. En lugar de generar columnas OHLC en tu objeto xts mensual, este ajuste tomará simplemente una fila de cada periodo como representativa de todo el periodo. Puedes indicar qué fila usar con el argumento indexAt.
Tanto los datos temps_xts como los datos flights_xts (del capítulo anterior) están disponibles en tu espacio de trabajo.
Este ejercicio forma parte del curso
Estudio de caso: Análisis de series temporales de una ciudad en R
Instrucciones del ejercicio
- Filtra tu objeto
temps_xtspara incluir solo observaciones de 2010 a 2015. Guarda el resultado comotemps_xts_2. - Usa
to.period()para convertir tus datos diarios de temperatura a periodicidad mensual. Asegúrate de especificar el periodo al que quieres convertir ("months"). También necesitas establecerOHLCenFALSEpara evitar generar nuevas columnas OHLC. Por último, establece el argumentoindexAten"firstof"para seleccionar la primera observación de cada mes. - Usa dos llamadas a
periodicity()para comparar la periodicidad y la duración de tus nuevos datos mensuales de temperatura con los datosflights_xtsdel capítulo anterior.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# 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(___)