I am comparing different sites based on their floristic composition in R. Therefore, I have created a huge community datamatrix (presence/absence data) from 53 sites including over 1000 species. To Perform PCoA on the data I first transformed the raw community datamatrix into a distance matrix using Bray-Curtis distance and took the square root of the resulting values to avoid negative eigenvalues. Than I performed PCoA using cmdscale()
and calculated the explained variance of the first Coordinates as followed:
spe.b.pcoa$eig[1]/sum(spe.b.pcoa$eig)
spe.b.pcoa$eig[2]/sum(spe.b.pcoa$eig)
As result I get explained variance of the first coordinate 9.3 % and second coordinate 7.3 % which appears to be extremely low. However, when I plot the results along the two axes it looks fairly nice and the plot shows exactly what I have expected. So, does anyone know why I get such low values of explained variance and if this really matters for me? Has this something to do with the huge amount of species which are, please correct me if I am wrong, in that case my variables?
Interestingly enough I found no PCoA plots on the internet which actually indicate the explained variance of the axes....
Thank you for your help!!!