'How to I plot the names of each sample region on top of a Canonical analysis of Principal Coordinates (CAP)

I am trying to include the region labels of the samples instead of plotting them just as numbers in the CAP plot. Essentially, how do I do I add labels to the site scores?

This is the code I've used.

'''
env<- read_excel("Environmental.xlsx")
bio<-read_excel("Biotic 3.xlsx")

bio2= wisconsin(sqrt(bio))
cap = capscale(formula = bio2 ~ ., data= env, dist = "bray")

plot(cap,xlim= c(-3,3), ylim=c(-3,3), display = "species", type ="text")
points(cap, display = "sites", col = 'grey')
text(cap3, display ='bp', col = 'blue')
'''

Thank you

r


Sources

This article follows the attribution requirements of Stack Overflow and is licensed under CC BY-SA 3.0.

Source: Stack Overflow

Solution Source