library(ggplot2)
library(dplyr)
library(tidyr)
library(tidyverse)
#install.packages("rlang")
#install.packages("ggpattern")
library(ggpattern)
 setwd("C:/Users/PAULA/Desktop/Datos doctorado/inercia/Archivos para paper")

my.data <- read.csv("Boxplot.csv", check.names = FALSE)


names(my.data) <- gsub(" ", "", names(my.data))
names(my.data) <- trimws(names(my.data))
names(my.data)[1] <- "Fecha"


my.data <- my.data[, colSums(!is.na(my.data)) > 0]

data_all <- my.data %>%
  pivot_longer(
    cols = -Fecha,
    names_to = c("Type","Species","Depth"),
    names_pattern = "(Tmax|Tmin)_(M|N)([0-9]+)",
    values_to = "Temp"
  ) %>%
  mutate(
    Depth = factor(Depth, levels = c("1","5","12")),
    Species = recode(Species,
                     "M" = "Modiola caroliniana",
                     "N" = "Nassella tenuis")
  )%>%
  filter(!is.na(Temp), is.finite(Temp))


labels_species <- c(
  "Modiola caroliniana" = expression(italic("Modiola caroliniana")),
  "Nassella tenuis"     = expression(italic("Nassella tenuis"))
)

# ---- gráfico ----
ggplot(data_all, aes(x = Type, y = Temp, fill = Species)) +
  geom_boxplot(alpha = 0.8) +
  facet_wrap(~Depth, labeller = label_both) +
  scale_fill_manual(
    values = c("Modiola caroliniana" = "grey40",
               "Nassella tenuis"     = "white"),
    labels = labels_species
  ) +
  labs(x = "", y = "Temperature (°C)", fill = "Species") +
  theme_minimal() +
  theme(legend.position = "bottom")