If you only have seven species, why not draw the label for the species (the names() component of your data) at the species score, rather than a colour?
plot(OBJ, display = c("sites","species"), scaling = 3, type = "n") points(OBJ, display = "sites", scaling = 3, type = "p") points(OBJ, display = "species", scaling = 3, type = "t") Not tested, am away at the moment but heading home. If you are still having trouble (and I see now why you wanted matplot) then email back and I'll given a longer example illustrating who to build up plots by hand. Gavin > Dear List, > > I am running constrained correspondence analysis for abundance data of 7 > birds. > > However, I would like to check which bird prefers which environment > gradient > by showing the species with different colors of the dots in cca plot > (package vegan). > > Please kindly help and thank you > > Elaine > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.