PEC 2 Raincloud Plot

if (!require('lubridate')) install.packages('lubridate');
library(lubridate)

path = 'estat_sdg_02_10_en.CSV'
df_obesity<- read.csv(path, row.names=NULL)

summary(df_obesity)

   DATAFLOW         LAST.UPDATE            freq               unit               bmi                geo             TIME_PERIOD     OBS_VALUE    
 Length:468         Length:468         Length:468         Length:468         Length:468         Length:468         Min.   :2008   Min.   : 5.90  
 Class :character   Class :character   Class :character   Class :character   Class :character   Class :character   1st Qu.:2014   1st Qu.:19.23  
 Mode  :character   Mode  :character   Mode  :character   Mode  :character   Mode  :character   Mode  :character   Median :2017   Median :36.50  
                                                                                                                   Mean   :2017   Mean   :35.87  
                                                                                                                   3rd Qu.:2019   3rd Qu.:50.38  
                                                                                                                   Max.   :2022   Max.   :64.80  
                                                                                                                                  NA's   :6      
   OBS_FLAG        
 Length:468        
 Class :character  
 Mode  :character  
                   
                   
                   
                   

df_obesity$bmi_fact <- as.factor(df_obesity$bmi)

# Crear el gráfico y asignarlo a 'P1'
P1 <- ggplot(df_obesity, aes(x = bmi_fact, y = OBS_VALUE, fill = bmi)) +
  geom_flat_violin(position = position_nudge(x = .25, y = 0), adjust = 2) +
  geom_point(position = position_jitter(width = .15), size = .25) +

  geom_boxplot(aes(x = as.numeric(bmi_fact) + 0.25, y = OBS_VALUE), 
               outlier.shape = NA, alpha = 0.3, width = .1, colour = "black") +
  
  ylab('OBS_VALUE') + xlab('bmi') +
  coord_flip() + theme_cowplot() +
  guides(fill = FALSE, colour = FALSE) +
  ggtitle("Raincloud Plot obesity")

# Mostrar el gráfico
P1