par(mfrow=c(1,3), mar=c(5,4,2,2))
plot(1:6, resR[1,],type="b", col="blue", ylim = range(resR), xlab = "Vegetation stade", ylab = "Richness")
 lines(1:6, resR[2,],type="b", lty=2, col="red")
 lines(1:6, resR[3,],type="b", lty=3, col="orange")
 lines(1:6, resR[4,],type="b", lty=4, col="darkgreen")

plot(1:6, resS[1,],type="b", col="blue", ylim = range(resS), xlab = "Vegetation stade", ylab = "Shannon index")
 lines(1:6, resS[2,],type="b", lty=2, col="red")
 lines(1:6, resS[3,],type="b", lty=3, col="orange")
 lines(1:6, resS[4,],type="b", lty=4, col="darkgreen")

plot(1:6, resGS[1,],type="b", col="blue", ylim = range(resGS), xlab = "Vegetation stade", ylab = "Simpson index")
 lines(1:6, resGS[2,],type="b", lty=2, col="red")
 lines(1:6, resGS[3,],type="b", lty=3, col="orange")
 lines(1:6, resGS[4,],type="b", lty=4, col="darkgreen")

legend(2.8, 0.6,  c("Pologne","Bourgogne","Provence","Corse"), lty = 1:4, col=c("blue", "red", "orange", "darkgreen"), bty="n")
par(mfrow=c(1,1))