I can't seem to prevent overlap, even after realigning the site labels. Is there any way I can prevent this overlap?
WD <- "C:/Users/..."
setwd(WD)
require(vegan)
phylum.dat <- read.csv("Condensed Phylum data for R.csv", stringsAsFactors = T, header = T, row.names = 1)
phylum.t1 <- decostand(phylum.dat, "total")
phylum.t2 <- decostand(phylum.t1, "max")
phylum.t2.dca <- decorana(phylum.t2)
phylum.t2.dca.DW <- decorana(phylum.t2, iweigh = 1)
names(phylum.t2.dca)
summary(phylum.t2.dca)
phylum.t2.dca.taxonscores <- scores(phylum.t2.dca,display=c("species"), choices=c(1,2))
plot(phylum.t2.dca, display = "none", type = "text", main = "Phylum")
point_colors <- c(rep("blue", 4), rep("red", 4), rep("green", 4))
points(phylum.t2.dca, display = "site", choices = 1:2, pch = c(15, 16, 17, 18), col = point_colors)
text(phylum.t2.dca, display = "species", choices = 1:2, cex = 0.8, col = "darkcyan")
text(phylum.t2.dca, display = "sites", cex = 0.8, pos = 4)
treat=c(rep("Treatment1",4), rep("Treatment2",4), rep("Treatment3",4))
ordihull(phylum.t2.dca,groups=treat, draw = "lines", col= c("blue", "red", "green"), label=F)