[R] Stress with MDS

Montignies Francois francois.montignies at insa-lyon.fr
Wed Feb 20 17:27:45 CET 2008


Hi,

I am looking for the best multidimensional configuration for my data (47*47
distance matrix).
I ve tried classical metric (cmdscale) and non metric MDS (isoMDS, nmds)
but it is now difficult to choose the best solution because of the
uncertainties in the definitions of the "stress" function.

So, same problem, several questions :

1. Statistical consideration : With "cmdscale" we get eigen values. What is
the best way to choose optimal dimensionality? With the eigen values and
corresponding GOF like we do with PCA ? If I compute stress1, does it make
any sense (I saw it in some publications)?

2. With isoMDS and nmds we get the final stress but i can't find the source
code so i don't know what is in the box. Obviously, I got different values
from isoMDS and  nmds . I started from the same initial configuration
(cmdscale) and the same parameters (maxit,tol)to compare results.
I tried to compute stress1 by myself and find values closed to nmds with
ndms config, but far away from isoMDS with isoMDS config (plus a strange
increasing value between k=4 and k=5). Could you help me please?  I lost
myself...

k	isoMDS$stress  stress1(isoMDS)  nmds$stress	stress1(nmds)
2	0,18830413	0.2912164	0.2758062	0.2658789
3	0,11521339	0.1866746	0.1754007	0.1727632
4	0,08733106	0.1638274	0.1281730	0.1271329
5	0,06942862	0.1991569	0.09756043	0.0970992
6	0,05751437	0.1563326	0.07846889	0.07822841

Here is my stress1 function

stress1<-function(datadist,fitteddist)
{sqrt(sum((datadist-fitteddist)^2)/sum(datadist^2))}

Best regards



More information about the R-help mailing list