Panel graphics in ggplot2

Question

I have 25 species and would like to

1. write only a script in which I could have a separate chart for each species (in total would be 25 small charts, in a grid 5x5).

2. I would also like to make two different Plots for savanna and forest;

3. put a trend line linking the average of Plots (see an example in figure 2).

Could someone help me with my script? figure 1 is the graph I made in R with the script that follows in the link below and figure 2 is an example of what I would like to do.

My data and script can be accessed by link: https://drive.google.com/open?id=1UBHOY6KdIu3SmxJRnK374rTsmYQCd5GK

Answer 1

I have a graphic suggestion just below. I will post my code, the result obtained and then make some comments:

library(ggplot2)

dados <- read.table(file="data_fran.csv", sep=";", header=TRUE)

dados$warfare <- factor(dados$warfare, levels=c("War1", "War2", "aPost"))

ggplot(dados, aes(x=warfare, y=Abund, group=warfare:habitat_F_S, fill=habitat_F_S)) +
  geom_boxplot() +
  stat_summary(fun.y=median, geom="line", lwd=1, 
    aes(group=habitat_F_S, colour=habitat_F_S)) +
  facet_wrap(~ Binomial, nrow=5) +
  labs(x="Guerra", y="Abundância", colour="Habitat", fill="Habitat") +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

Now I will give the explanations and justifications for this graph:

• ggplot(dados, aes(x=warfare, y=Abund, group=warfare:habitat_F_S, fill=habitat_F_S)): determines which data set will be plotted, as well as the x and y axes and colors. Also, I created the interaction warfare:habitat_F_S, for the data to be grouped by warfare and habitat_F_S simultaneously

• geom_boxplot(): I defined that I wanted a boxplot

• stat_summary(fun.y=median, geom="line", lwd=1, aes(group=habitat_F_S, colour=habitat_F_S)): uni the medians of the boxplots, not the averages. Particularly, I think it looks better this way, because the horizontal line in the middle of the boxplot is the median. If you insist on the average, as asked in the question, change fun.y=median for fun.y=mean.

• facet_wrap(~ Binomial, nrow=5): this is the command that creates the grid with 5 lines for each species

• labs(x="Guerra", y="Abundância", colour="Habitat", fill="Habitat"): I changed the axes and caption names to make them more explanatory

• theme(axis.text.x = element_text(angle = 45, hjust = 1)): I rotated the x-axis Abels at 45 degrees for aesthetic reasons