Dear All,
I am using a Procrustes analysis to compare two NMDS ordinations for the
same set of sites. One ordination is based on fish data, the other is based
on invertebrate data. Ordinations were derived using metaMDS() from the
{vegan} library as follows:
fish.mds<-metaMDS(fish.data, distance="bray", k=3, trymax=100,
wascores=TRUE, trace=TRUE, zero="add")
invert.mds<-metaMDS(bugcal.a, distance="bray", k=3, trymax=100,
wascores=TRUE, trace=TRUE, zero="add"),
where fish.data and invert.data are site-by-abundance matrices for fish and
invertebrates, respectively.
I have then used protest() to test the significance between the two
ordinations:
protest.results<-protest(fish.mds, invert.mds, scale=TRUE, symmetric=TRUE,
permutations=999)
The problem is, the results of the protest analysis seem to be affected by
the original ordering of sites in the data input files. That is, if I
re-sort one of the input files based on some criteria, I sometimes see a
change in the strength and significance of the concordance results. In an
attempt to figure out what is happening, I have re-ordered each dataset in a
number of ways, and plotted the resulting ordinations. I have seen that
while the ordering of sites in the input file does not change the spatial
relationship between them (i.e., does not change the fundamental ordination
solution), it can change how the sites are oriented with respect to the
different NMDS axes. I was of the belief that a difference in orientation
with respect to the NMDS axes should not affect the Procrustes comparison of
two ordinations, as the protest function should be rotating and reflecting
one matrix with respect to the other to find the closest match between them
(i.e., it should be accounting for differences in how the two solutions are
oriented in space). However, I appear to see different results depending on
how sites are oriented with respect to each NMDS axis.
When comparing two ordination solutions with Protest, is it necessary to
ensure that the original data input files are ordered in the same way?
Thanks in advance for any insight.
Sincerely,
--
Christy Dolph
Ph.D. Candidate
Water Resources Science
University of Minnesota
200 Hodson Hall, 1980 Folwell Ave
St. Paul, MN 55108
dolph...@umn.edu
612-868-1565
[[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.
