Hola Mundo
El FPR mide la intensidad del incendio. Pero cómo se puede identificar un FPR es alto? La mayoría ocurre entre 10 y 40. Pero entonces los de más de 40 sonMUY graves?
Comenta una 🍍 en el tweet si aprendiste algo nuevo!
El código para hacer la gráfica!
# Hecho por @GonzalezGouveia con gusto para
# Twitter #DatosDeMiercoles por @R4DS_es
# cargar paquetes
library(tidyverse)
# leer los datos de internet
incendios <- read_csv("https://bit.ly/2lRiAiB")
# filtrando para intensidad mayor que cero
incendios_filtered <- incendios %>%
filter(intensidad > 0)
# calculando cuantiles
quant_lim <- quantile(incendios_filtered$intensidad)
# mirando los cuantiles se puede saber
# la distribución de los datos
# etiquetas
bajo <- 'FRP mayor a 10'
medio <- 'FRP entre 10 y 40'
muy_alto <- 'FRP menor a 40'
# creando las etiquetas basados en los cuantiles
incendios_filtered <- incendios_filtered %>%
mutate(intensidad_lab = case_when(
intensidad <= 10 ~ bajo,
intensidad <= 40 ~ medio,
TRUE ~ muy_alto
))
# ordenando los factores para el facet
incendios_filtered$intensidad_lab <-
factor(incendios_filtered$intensidad_lab,
ordered = TRUE,
levels = c(bajo, medio, muy_alto),
labels = c(bajo, medio, muy_alto))
# otra prueba
incendios_filtered %>%
group_by(intensidad_lab,
semana) %>%
summarise(n = n()) %>%
ggplot(aes(x = semana,
y = n,
color = intensidad_lab)) +
geom_point() +
geom_smooth() +
facet_wrap(~ intensidad_lab) +
labs(title = 'Incendios en el Amazonas',
subtitle = 'por semana e intensidad',
y = 'cantidad',
x = 'semanas desde el 1ero enero 2019',
caption = 'FRP (fire radiate power, en megawatts) medida de intensidad del incendio') +
theme(legend.position = 'none')
ggsave('ddm_amazonas.png',width = 6, height = 3)
# fin del código
# Sigueme @GonzalezGouveia --> https://bit.ly/2WkLmFb