site.pres<- rowSums(spe>0) # Nombre d'espèces présentes dans chaque site
hist(site.pres, main=”Richesse en espèces”, las=1, xlab=”Fréquence des sites”, ylab=”Nombre d'espèces”, breaks=seq(0,30, by=5), col=”grey”)