similar to: best ordination method for binary variables

Hi, I'm trying to do a non-metric multidimensional scaling using isoMDS. However, I have some '0' distances in my data, and I'm not sure how to deal with them. I'd rather not drop rows from the original data, as I am comparing several datasets (morphology and molecular data) for the same individuals, and it's interesting to see how much morphological variation can be

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,

My question is with respect to mCLUST and the values of BIC and log likelihood. The relevant part of my R script is: ######################### BEGIN MDS ANALYSIS ######################### #load data data <- read.table("Ecoli33_Barry.dis", header = TRUE, row.names = 1) #perform MDS Scaling mds <- metaMDS(data, k = Dimensions, trymax = 20, autotransform =TRUE, noshare = 0.1,

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,

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.

Hi the list, I need to include some C code in R, but the behavior of the C code is strange : Here is my code : --- 8< --- Rprintf("\n XXXX mTraj=%f mClus=%f",mTraj[i+nbId*c],mClustersCenter[j+nbClusters*c]); Rprintf("\nDistA=%d Tmp=%d",dist,tmp); tmp = mTraj[i+nbId* c] - mClustersCenter [j+nbClusters* c]; Rprintf("\nDistB=%d

Hi, I could use a hand solving a fairly straightforward ordination plot problem: I am conducing NMDS on some community data for roughly 300 localities and 650 species. I have a community matrix, a species attribute matrix, and an environmental attributes matrix. After running metaMDS {vegan} on the community matrix, I can successfully use either of the other two matrices to draw hulls around

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

I am severely rusty re. multivariate / ordination analysis, having done my last work 40 years ago (in plant ecology). I am interested in exploring applications of multivariate analytic approaches to data from the history of motion picture films. I'd very much appreciate any pointers as to possibly appropriate proceedures. The individual "samples" may be individual films and the

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

Dear List, I'm trying to interpret the results of the Kruskal's Non-metric Multidimensional Scaling algorithm (isoMDS, MASS package). The 'goodness of fit' is reported as "The final stress achieved (in percent)". What does this mean exactly? I've tried to google for an answer but I've not come up with a definitive answer. Regards, Dieter -- Dieter Vanderelst

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

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

Hi, I am working with a species-by-trait .csv file (columns=traits, rows=species) and get the following warning message when trying to plot results of both metaMDS and pcoa:  "Warning message: In ordiplot(x, choices = choices, type = type, display = display,  :   Species scores not available" I am using a Gower's transformation in both procedures within the metaMDS or pcoa

Dear all, I'm completely new to R and at first I must say that it is a great program! But I have a problem with the function isoMDS from the MASS package. I have this code which I load with source() from a file: x <- c(163.59514923926784, 150.01448475257115, ...... {here are some more values}) x.sort <- sort(x) x.dist <- dist(x.sort) library(MASS) x.mds <- isoMDS(x.dist)

Hi, I'm trying to use isoMDS to project a directed graph to 2-dim vectors, but I got an error. #here is the code to create the graph using igraph package and run isoMDS on it. library(igraph) library(MASS) g<-make_graph(c(1,2, 2,3, 2,4, 3,4, 4,5, 5,6, 3,6, 1,6, 2,5),directed=TRUE) dist<-distances(g, mode="out") loc<-isoMDS(dist) # below is content of the dist matrix

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

I am starting with a matrix in which rows are vegetation plots and columns are various characteristics including ID# and elevation. I removed elevation and ID columns to avoid having those characteristics influence the distances between points which I calculated using the "dist" command. The resulting distance file was then used in isoMDS. What I want to know is whether I can

If tbl is an object of class 'dist', you can do this: a <- sammon(tbl, k=3) But you can't do this: b <- isoMDS(tbl, k=3) Wouldn't it be sensible to have identical interfaces to sammon() and isoMDS() ? I think all that would be needed is to change this: isoMDS <- function(d, y=cmdscale(d, 2), maxit=50, trace=TRUE) { ...into this: 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,