spe.pres<- colSums(spe>0) # Somme des sites où chaque espèce est présente. 
hist(spe.pres, main=”Cooccurrence des espèces”, las=1, xlab=”Fréquence”, breaks=seq(0,30, by=5), col=”grey”)