# Figure showing number of vertebrate taxa detected for 3 methods (visual survey, camera trap, eDNA)
# for three vegetation types

library(ggplot2)
library(ggsci)

taxa_detected_data <- read.csv("vertebrate_detections_vegtype.csv", header = TRUE)

p_taxa <- ggplot(taxa_detected_data, aes(fill=Categories, y=Number, x=Vegetation)) + 
  geom_bar(position="stack", stat="identity") + 
  xlab("Vegetation type") +
  ylab("Number of taxa") +
  scale_y_continuous(expand = c(0,0),
                     limits = c(0,25)) +
  scale_x_discrete(limits=c("Willow", "Maple", "Conifer")) +
  theme_bw() + theme(panel.border = element_blank(), panel.grid.major = element_blank(),
                     panel.grid.minor = element_blank(), axis.line = element_line(colour = "black")) +
  theme(axis.text.x = element_text(colour = "black")) +
  theme(axis.text.y = element_text(colour = "black"))

p_taxa + scale_fill_npg(name = "Detection method")