# Représentation du nombre de sites fouillés par décennie ###

ensembles <- as_tibble(read.csv("ensembles_général.csv", sep = ";",
header = TRUE))

summary(ensembles)

ensembles.sans.dbl <- distinct(ensembles, ensemble_id, .keep_all = TRUE)

ensembles.sans.dbl$décennie_fouille <-
recode(ensembles.sans.dbl$décennie_fouille,
        "1970?" = "non rens.",
        "1980?" = "non rens.",
        "1990?" = "non rens.")

ensembles.sans.dbl <-
ensembles.sans.dbl[ensembles.sans.dbl$décennie_fouille != "non
rens.",]

ensembles.sans.dbl$décennie_fouille <-
as.factor(ensembles.sans.dbl$décennie_fouille)
summary(ensembles.sans.dbl$décennie_fouille)

dec.fouille <- ggplot(data = ensembles.sans.dbl, mapping = aes(x = as.Date(décennie_fouille, format = "%Y"))) +
  geom_histogram bins = 30, fill = "#46ACC8", color = "black") +
  scale_fill_manual(values = wes palette("FantasticFox1", n = 1)) +
  scale_x_date(breaks = "10 years", labels = date_format("%Y"),
  limits = c(as.Date("1868-05-20"), as.Date("2026-05-20"))) +
  ylab("Fréquence") +
  xlab("Année") +
  theme_classic() +
  theme (axis.text = element_text(size = 12),
  axis.title = element_text(size = 14, face = "bold"))
dec.fouille

ggsave(file = "histogramme_dates_fouilles.svg", width = 8, height = 6)