I have 4 sampling locations and 23 species abundance data. Using the same concept and formula of Simpson index of diversity, instead of calculating the probability of individuals of different species present in a community (sampled site), I calculated the probability of occurrence of individuals of a single species among the communities (sampling sites) in R (vegan package). My outcomes are very satisfactory with the total abundance results.