spe.pres <- colSums(spe>0) #compute the number of sites where each species is present. 
hist(spe.pres, main="Species occurrence", las=1, xlab="Frequency of occurrences", breaks=seq(0,30, by=5), col="grey")