similar to: Problem in NMDS

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,

Hi, I've been using R for a while now but I've got a problem with metaMDS (in the vegan package) that I can't quite figure out. I have a set of proportion data (from 0-1, rows sum to 1) that I apply metaMDS to using the command: nMDS.set=metaMDS(sqrt(test.set),distance="euclidean",k=3,zerodist="add",autotransform=FALSE) I am using a squared-chord distance

Dear List, This is Elaine, a postgraduate studying in bird distributions in East Asia. I want to calculate Simpson dissimilarity index, based on a presence/absence matrix of bird species in islands in East Asia. (matrix row: 36 islands/matrix column: species ID) (R package vegan to make NMDS and R package betapart) In most papers using vegan for NMDS and betapart for dissimilarity

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

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

Hi R-experts, I am using the envfit function over an ordination of floristic data. The problem is that every time that I run it changes the results. Sometimes dramatically, selecting variables that the first time were not significant. I do not get what could be the problem or if is normal given the permutations are different. # the NMDS ordination gap_flor_NMDS_chord <- metaMDS(gaps_flor,

Hello R-user community! I am running R 2.7.0 on a Power Book (Tiger). (I am still R and statistics beginner) Presently I try to run the function metaMDS (vegan) using an existing dissimilarity-matrix. As I would like to start with this matrix I thought I could just give the matrix using the x= -argument Test<-metaMDS(x=Dist.Gower) Fehler in inherits(comm, "dist") :

Hi! Thanks for providing great help in R-related statistics. Now, however I'm stuck. I'm not a statistics person but I was recommended to use R to perform a nmds plot and PerMANOVA of my dataset. Sample(treatment) in the columns and species (OTU) in the rows. I have 4 treatments (Ambient Temperature, Ambient temperature+Low pH, High temperature, High temperature+low pH), and I have 16

Dear r-helpers! How can I integrate other distances (in the form of a dist object) into function metaMDS? The problem: metaMDS needs the original data.frame for the calculation and only the default distances of function vegdist are allowed. Any suggestions are greatly appreciated! Thank you, Kim -- Jetzt kostenlos herunterladen: Internet Explorer 8 und Mozilla Firefox 3.5 - sicherer, schneller

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,

Dr. Stevens, Hi, my name is Trey Scott, and I'm a grad student of Brian McCarthy's. He referred me to you because of your expertise in handling complex R problems. We were hoping you could help us solve a nagging problem that is prohibiting me from producing graphicl output. Here is a simple mock-up of the matrix I'm using a b c d e f 1i 1 4

Hi, I am trying to create 3 mds plots: one with axis 1 vs axis 2, one with axis 2 vs axis 3, and one with axis 1 vs axis 3. When inputting my code, I only end up with one diagram and when inputting mdsg.mds$dims, the program returns 2 for 2 dimensions. How can I create the other two plots? Any help would be greatly appreciated, Calla Carbone The table I use is number of each taxa on each

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

Hello, We are trying to find a function to compute "stress" in our multidimensional scaling analysis of a dissimilarity matrix. We've used "dist()" to create the matrix and "cmdscale()" for the scaling. In order to determine the number of dimensions we would like to plot stress vs. dimensions. However, we cannot find a pre-made command. It seems that other

Hi, I'm trying to use metaMDS with a dissimilarity matrix of angles, not Bray-Curtis, and I wanted to know if there is an in-built function to produce a plot of stress values against dimensions, that could be used to determine the 'true' dimension of the solution. The number of objects is only a lowly 8 so any solution higher than 2-dimensional is not likely to be interpretable,

H ello R-experts, I want to do ordination plots using vegan metaMDS. I have a where many cells have zero values. Data structure: X[1:10,1:14] Height.1 Height.2 Height.3 Height.4 Height.5 Height.6 Height.7 Height.8 Height.9 Height.10 Height.11 Height.12 Height.13 D30I1A 46 0 0 0 0 0 0 0 0 0 39 0 98 D30I1B

Dear R-community, dear Jari Oksanen! I use metaMDS (package vegan) to calculate NMDS. In a lot of papers I read that it is recommended to use previous best solutions as a new starting configuration to get better results and to avoid local minima. On the help page I found that a previous.best-command is already implemented in metaMDS: metaMDS(comm, distance = "bray", ...,plot = FALSE,

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.,

Dear r-Community, Step1: I would like to calculate a NMDS (package vegan, function metaMDS) with species data. Step2: Then I want to plot environmental variables over it, using function envfit. The Problem: One of these environmental variables is cos(EXPOSURE). But for flat releves there is no exposure. The value is missing and I can't call it 0 as 0 stands for east and west. Therefore I

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")