[R] procrustes (vegan) plot of residual differences
Juan Antonio Balbuena
j.a.balbuena at uv.es
Wed May 23 17:19:30 CEST 2012
Hello
This is a simple question but I couldn't google an answer.
In the procrustes function of the vegan package, one uses
plot(procrustes_object, kind=2) to obtain a plot of the residual
differences. For instance:
data(varespec)
vare.dist <- vegdist(wisconsin(varespec))
library(MASS)
mds.null <- isoMDS(vare.dist, tol=1e-7)
mds.alt <- isoMDS(vare.dist, initMDS(vare.dist), maxit=200, tol=1e-7)
vare.proc <- procrustes(mds.alt, mds.null)
plot(vare.proc, kind=2)
In the resulting plot 3 horizontal lines (1 continuous and 2 dashed) are
displayed and I was wondering what they actually represent. (Perhaps mean
residual +/- 95% CI ?).
I would also like to know whether it is possible to customize the plot by,
for instance, removing the lines, changing their appearance, or making them
represent something else, such as mean +/- SD.
Thank you very much for your attention.
Juan Antonio Balbuena
--
Dr. Juan A. Balbuena
Marine Zoology Unit
Cavanilles Institute of Biodiversity and Evolutionary Biology
University of
Valencia
[1]http://www.uv.es/~balbuena
P.O. Box 22085
[2]http://www.uv.es/cavanilles/zoomarin/index.htm
46071 Valencia, Spain
[3]http://cetus.uv.es/mullpardb/index.html
e-mail: [4]j.a.balbuena at uv.es tel. +34 963 543 658 fax +34 963 543 733
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
NOTE! For shipments by EXPRESS COURIER use the following street address:
C/ Catedrático José Beltrán 2, 46980 Paterna (Valencia), Spain.
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
References
1. http://www.uv.es/%7Ebalbuena
2. http://www.uv.es/cavanilles/zoomarin/index.htm
3. http://cetus.uv.es/mullpardb/index.html
4. mailto:j.a.balbuena at uv.es
More information about the R-help
mailing list