Search code examples
rconstraintspcaveganvegdist

how to preform a PCoA with constraind factors

I have performed the following PCoA, using the following code:

    spe.dist <- vegdist(spe, method = "bray")
spe.braysq.pcoa <- pcoa(sqrt(spe.dist))
biplot.pcoa(spe.braysq.pcoa, sqrt(spe))

However, I now want to perform a constrained PCoA. I want to constrain it for the factors field and moment_datacollection and have the other factors still unconstrained. I have searched on the internet to find a way to do this. However, I could not find it. Do you guys know how to do this? Or did I set up the code the wrong way?

Solution

  • Use vegan function dbrda (which is named after acronym for distance-based redundancy analysis).