Ozono y una introducción a ggplot2
El ozono atmosférico se encuentra más que nada en la estratósfera, y su concentración varía a lo largo del año. Su presencia es clave porque absorbe radiación ultravioleta dañina. En los años '80 se hizo evidente que las concentraciones de ozono estaban en baja, y que la causa eran los CFC producidos por la industria. En consecuencia, se acordó terminar con su uso y reemplazarlos con productos más amigables. Se cree que los niveles de ozono se van a recuperar en 2030-2040. Este hecho es uno de los hitos más importantes en la historia ambiental moderna, y representa un antecedente clave para otro tipo de cambios que tenemos que hacer ya, vinculados al cambio climático y pérdida masiva de biodiversidad.
El paquete de R {ggplot2} es enorme y versátil. A continuación, apenas algunas de sus virtudes con un ejemplo sencillo de series temporales. Los datos corresponden a mediciones de ozono atmosférico tomadas en Perth, Australia, entre los años 1984 y 2016. Los obtuve de un laboratorio del NOAA. Una novedad es que los bajé y pre-procesé por línea de comando en linux/bash, lo que van a encontrar como primer script más abajo. Luego, por medio del paquete {plyr} en R, calculé promedio y desvío estándar por mes y por año (también detallado en scripts comentados al final).
En {ggplot2}, si uno usa los parámetros estándares, el gráfico se ve así...
Todo los detalles pueden ser ajustados desde 'theme()', pero también existe una opción global, 'theme_set()', con la que uno puede cambiar el estilo general del gráfico y el tamaño de la letra. El que sigue es de estilo 'classic'...
Hace poco empecé a armar arreglos de gráficos directamente desde R, algo muy útil para publicaciones. En el ejemplo que sigue, hecho con {ggpubr}, muestro el mismo gráfico con diferentes opciones para los colores...
Por último, presento los promedios anuales de ozono entre 1984 y 2016 con el estilo de gráfico 'bw', y con rótulos del eje X ajustados para mejor visualización. Todavía no se vislumbra la recuperación en la concentración de ozono...
Script de pre-procesamiento
Script para graficar
Fuente
Earth System Research Laboratory, NOAA
El paquete de R {ggplot2} es enorme y versátil. A continuación, apenas algunas de sus virtudes con un ejemplo sencillo de series temporales. Los datos corresponden a mediciones de ozono atmosférico tomadas en Perth, Australia, entre los años 1984 y 2016. Los obtuve de un laboratorio del NOAA. Una novedad es que los bajé y pre-procesé por línea de comando en linux/bash, lo que van a encontrar como primer script más abajo. Luego, por medio del paquete {plyr} en R, calculé promedio y desvío estándar por mes y por año (también detallado en scripts comentados al final).
En {ggplot2}, si uno usa los parámetros estándares, el gráfico se ve así...
Todo los detalles pueden ser ajustados desde 'theme()', pero también existe una opción global, 'theme_set()', con la que uno puede cambiar el estilo general del gráfico y el tamaño de la letra. El que sigue es de estilo 'classic'...
Hace poco empecé a armar arreglos de gráficos directamente desde R, algo muy útil para publicaciones. En el ejemplo que sigue, hecho con {ggpubr}, muestro el mismo gráfico con diferentes opciones para los colores...
Por último, presento los promedios anuales de ozono entre 1984 y 2016 con el estilo de gráfico 'bw', y con rótulos del eje X ajustados para mejor visualización. Todavía no se vislumbra la recuperación en la concentración de ozono...
Script de pre-procesamiento
Script para graficar
Fuente
Earth System Research Laboratory, NOAA
Comentarios