[R] Estimate "between-axes" vs "within-axes heterogeneity of multivariate matrices
Nikos Alexandris
nikos.alexandris at felis.uni-freiburg.de
Thu Dec 23 06:13:27 CET 2010
This time with a more-R oriented question:
Is the mrpp {vegan} package [1] useful in trying to check, or get a clue about
the differences between- and within-axes (or variables or dimensions or
columns) of a multivariate matrix?
The description explains:
" ...(MRPP) provides a test of whether there is a significant difference
between two "or more groups of sampling units. ..."
"... difference may be one of location (differences in mean) or one of spread
(differences in within-group distance) ..."
and
"... Function mrpp operates on a data.frame matrix where rows are observations
and responses data matrix. The response(s) may be uni- or multivariate. ..."
Question: what about the observations being actually the columns? Is a simple
transposing of the matrix enough? Any other alternatives or hints?
Thak you, Nikos
---
[1] <http://cc.oulu.fi/~jarioksa/softhelp/vegan/html/mrpp.html>
--%<---
Nikos:
> My question(s) in the end might be silly but I am no expert on this, so
> here it goes:
>
> Noy-Meir (1973), Pielou (1984) and a few others have pointed to
> non-centered PCA being in some cases useful. They clearly explain that "it
> is the case" when multi-dimensional data display distinct clusters (which
> have zero, or near-zero, projections in some subset of the axes) and the
> task is (exactly) to separate this clusters among the principal
> components.
>
> I have done my complete work using prcomp() and tested combinations of
> center=FALSE/TRUE and scale=FALSE/TRUE. I would like to now check this
> "between-axes" vs "within-axes" heterogeneity of my data and cross-check
> results with the various tested PCA-versions.
>
> Is there any (official or custom) function available in R that could answer
> this question? Some relative/comparative (preferrable simple and intuitive)
> measure(s)? Something that would graphically perhaps give an indication
> without time-consuming clustering, sampling or whatsoever processing?
>
> Even though the above mentoined authors mention some measure for the
> assymetry of the yielded compoenents ( uncentered -> unipolar, centered ->
> bipolar) I find the concept a bit hard to understand.
>
> Isn't there a quick way (function) to just say (with numbers of plots of
> course) "well, it seems that the data are heterogenous looking at between-
> axes" or the other way around "it looks like the variables differ within,
> more than between"?
>
> Apologies for repeating the same question (trying to understand the problem
> myself). Thank you, Nikos
More information about the R-help
mailing list