Eu tenho uma sugestão de gráfico logo abaixo. Vou postar meu código, o resultado obtido e a seguir fazer alguns comentários:
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))
Agora vou dar as explicações e justificativas para este gráfico:
ggplot(dados, aes(x=warfare, y=Abund, group=warfare:habitat_F_S, fill=habitat_F_S))
: determina qual conjunto de dados será plotado, bem como os eixos x e y e as cores. Além disso, criei a interação warfare:habitat_F_S,
para que os dados fossem agrupados por warfare
e habitat_F_S
simultaneamente
geom_boxplot()
: defini que queria um boxplot
stat_summary(fun.y=median, geom="line", lwd=1, aes(group=habitat_F_S, colour=habitat_F_S))
: uni as medianas dos boxplots, não as médias. Particularmente, acho que fica melhor assim, pois a linha horizontal no meio do boxplot é a mediana. Caso tu faça questão que seja a média, como pedido na pergunta, altere fun.y=median
para fun.y=mean
.
facet_wrap(~ Binomial, nrow=5)
: este é o comando que cria o grid com 5 linhas para cada espécie
labs(x="Guerra", y="Abundância", colour="Habitat", fill="Habitat")
: alterei os nomes dos eixos e legenda para ficarem mais explicativos
theme(axis.text.x = element_text(angle = 45, hjust = 1))
: rotacionei os labels do eixo x em 45 graus por questões estéticas