#Species distribution
(ab <- table(unlist(spe))) #note that when you put an entire line of code in brackets like this, the output for that operation is displayed right away in the R console

barplot(ab, las=1, xlab="Abundance class", ylab="Frequency", col=grey(5:0/5))