similar to: Stress with MDS

Dear useRs, last week I asked you about a problem related to isoMDS. It turned out that in my case isoMDS was trapped. Nonetheless, I still have some problems with other data sets. Therefore I would like to know if anyone here has experience with how well isoMDS performs in comparison to other non-metric MDS routines, like Minissa. I have the feeling that for large data sets with a high

I have been trying to use lavaan (version 0.4-7) for a simple path model, but the program seems to be computing far less degrees of freedom for my model then it should have. I have 7 variables, which should give (7)(8)/2 = 28 covariances, and hence 28 DF. The model seems to only think I have 13 DF. The code to reproduce the problem is below. Have I done something wrong, or is this something I

Hi; I am working with the similarity matrix below and I would like to plot a two-dimensional MDS solution such as each point in the plot has a label. This is what I did: data <- read.table('c:/multivariate/mds/colour.txt',header=FALSE) similarity <- as.dist(data) distance <- 1-similarity result.nmds <- nmds(distance) plot(result.nmds) (nmds and plot.nmds as defined at

Dear all, I'm running a set of nonparametric MDS analyses, using a wrapper for isoMDS, on a 800x800 distance matrix. I noticed that setting the parameter k to larger numbers seriously increases the calculation time. Actually, with k=10 it calculates already longer than for k=2 and k=5 together. It's now calculating for 6 hours, and counting... There is quite a difference between the

Hello, I am using vegan to do an NMDS plot and I would like to suppress the labels for the loading vectors. Is this possible? Alternatively, how can I avoid overlap? Many thanks for the help. Example code: #perform NMDS using metaMDS() function spe.nmds<-metaMDS(data, distance='bray',k=2 , engine = "isoMDS", autotransform=F, trymax=1000) #calculate the loading (i.e.,

hello, subsequently to a NMDS analysis (performed with metaMDS or isoMDS) is it possible to rotate the axis through a varimax-rotation? Thanks in advance. Bernd Panassiti

Dear R-helpers, I am running metaMDS in the vegan package, which uses isoMDS in MASS, to perform Nonmetric Multidimentional Scaling (NMDS). I have seen some authors report a p-value for the NMDS ordination based on randomization of the dataset. As I understand it this is meant to compare the stress in your dataset to multiple runs of randomized data. I do not see a way to perform such a test in

Hi, I'm using nmds command (library vegan) to analyze some fishing data. I'd like to plot not only points, but also the names of species and stations in a specified position. I used the command text(nmds$points[,1], nmds $points[,2],labels=row.names(nmds $points),pos=3,cex=0.5) But the labels are sometimes overlapped. Is there any way to use identify, or a similar command, to plot the

What objective function is maximized by survreg with the default Weibull model? I'm getting finite parameters in a case that has the likelihood maximzed at Infinite, so it can't be a simple maximum likelihood. Consider the following: ############################# > set.seed(3) > Stress <- rep(1:3, each=3) > ch.life <- exp(9-3*Stress) > simLife <- rexp(9,

Dear All, I was wondering whether anyone might be able to provide some advice with an nMDS / R problem. I?m trying to run nMDS on a dataset that contains many missing values and was wondering how I can account for the missing values when running nMDS? It seems as though the data are being grouped depending on where the zero values appear. Any suggestions greatly appreciated. Thank you very much

I've been unable to find a R package that provides the means of performing Clarke & Ainsworth's BIO-ENV procedure or something comparable. Briefly, they describe a method for comparing two separate sample ordinations, one from species data and the second from environmental data. The analysis includes selection of the 'best' subset of environmental variables for explaining

Hi I have just (finally) started to poke around in R and wanted to analyse a stream fish dataset with 28 sites and 18 species. When trying to follow the Vegan manual to run nmds from distance measures calculated by the vegdist function it turns out that I have two sites (streams) with the exactly the same four species (I have used pres-abs data in this case). When I try to run isoMDS I get an

Hello, I can draw a basic stress plot for NMDS with the following code in package Vegan. > stressplot(parth.mds, parth.dis) When I try to specify the line and point types, it gives me error message. > stressplot(parth.mds, parth.dis, pch=1, p.col="gray", lwd=2, l.col="red") Error in plot.xy(xy, type, ...) : invalid plot type In the above code, if I removed line type,

Hi, I am currently working on some data and feel that NMDS would return an excellent result. With my current data set however I have been experiencing some problems and cannot carry out metaMDS. I have tried with a few smaller data sets which I created for practice sake and this has worked fine. I think it is the set up of my data set that is causing me trouble. I have 18 columns and 18 rows,

Hello everyone! metaMDS(cm, distance = "euclidean", k = 2, trymax = 50, autotransform =TRUE, trace = 1, plot = T) (cm is a similarity matrix, in which values are positive integers or 0) I use this command to run NMDS on my matrix "cm". But the stress is very high after analysis. About 14. Actually, there is no improvment comparing with using isoMDS.

Colleagues, I'm developing a library of functions for community ecology analyses and have a couple of questions that I've not been able to answer via faqs or docs. 1) Are there existing functions for: a) Bray-Curtis (polar) ordination? b) non-metric multidimensional scaling (NMDS)? c) canonical correlation analyses (CCA and DCA)? d) TWINSPAN (doubtful...)? 2) If not,

Hello everyone, I have a question on isoMDS. My data set (of vegetation) with 210 samples is in this way: Rotfoehrenau Lavendelweidenau Silberweidenau .... 067_Breg.7 0 2 0 .... 071_Dona.4 0 2 6 .... ... I want to do an isoMDS-analysis with the dissimilarity index

Hello R experts, I've used metaMDS to run NMDS on some fish abundance data, and am also working on correlating environmental data to the NMDS coordinates. I'm fairly new to metaMDS and NMDS in general, so I have what are probably some very basic questions. My fish abundance data consists of 66 sites for which up to 20 species of fish were identified and counted. I ran metaMDS on this data

Dear r-helpers, I just read in an article by Virtanen et al. (2006) where vegetation-environment relationships are studied by fitting smoothed surfaces on an NMDS ordination using GAMs (Wood 2000). The authors describe, that they used R? as goodness-of-fit statistic, which they compare to the R? of fitted vectors. Calculations were carried out using the package vegan (Oksanen). I know that I can

Dear R users, could you please help me figure out why I am getting an error? Initially my data looks like this: > attributes(compl)$names [1] "UserID" "compl_bin" "Sex.x" "PHQ_base" "PHQ_Surv1" "PHQ_Surv2" "PHQ_Surv3" [8] "PHQ_Surv4" "EFE"