setwd("C:/Users/PAULA/Desktop/Datos doctorado/inercia/Archivos para paper")

require(ggplot2)
library(tidyr)
library(dplyr)
read.csv("coverage.csv")
my.data<-read.csv("coverage.csv")
colnames(my.data)[1] <- "Month"
long.data <- my.data %>%
  pivot_longer(cols = -Month,
               names_to = "Species",
               values_to = "Coverage")
long.data$Species <- recode(long.data$Species,
                            "M..caroliniana" = "Modiola caroliniana",
                            "N..tenuis" = "Nassella tenuis")
long.data$Species <- recode(long.data$Species,
                            "Modiola caroliniana" = "Mitalic",
                            "Nassella tenuis" = "Nitalic")

long.data$Month <- factor(long.data$Month,
                          levels = c("June-2024","July-2024","August-2024","September-2024",
                                     "October-2024","November-2024","December-2024",
                                     "January-2025","February-2025","March-2025","April-2025",
                                     "May-2025","June-2025"))
ggplot(long.data, aes(x = Month, y = Coverage, color = Species, shape = Species, group = Species)) +
  geom_line(linewidth = 1.2) +
  geom_point(size = 3) +
  scale_color_manual(
    name = NULL,  
    values = c("Mitalic" = "grey70", "Nitalic" = "black"),
    labels = c(expression(italic("Modiola caroliniana")),
               expression(italic("Nassella tenuis")))
  ) +
  scale_shape_manual(
    name = NULL,  # 👈 mismo "name", así se fusionan
    values = c("Mitalic" = 16, "Nitalic" = 17),
    labels = c(expression(italic("Modiola caroliniana")),
               expression(italic("Nassella tenuis")))
  ) +
  labs(x = "Month", y = expression("Coverage (m"^2*")")) +
  theme_minimal(base_size = 11) +
  theme(panel.grid = element_blank(),
    axis.text.x = element_text(angle = 45, hjust = 1),
    legend.title = element_blank(),   # 👈 asegura que no aparezca título
    legend.position = "bottom"
  ) +
  scale_x_discrete(labels = c("Jun","Jul","Aug","Sep","Oct","Nov","Dec",
                              "Jan","Feb","Mar","Apr","May","Jun"))