similar to: Clustering analysis with ordination plots

Hi all, I have site data with plant species cover and am looking for trends. I'm kind of new to this, but have done lots of reading and can't find an answer. I tried decorana (I know it's been replaced by ca.) and see a trend, but I'm not sure what it means. Is there a way to get the loadings/eigenvectors of the axes (like in PCA)? Is there a way to do this with rda() too? How

Hello, I am trying to learn how to do PCA. I found a tutorial online, but what I'm seeing in my installation does not match what is in the tutorial. Specifically, if I select Statistics:Dimensional Analysis: Principal-components analysis, the dialog I am presented with does not include an option to create an ordination plot. My question is why not? Is this a problem with my data set,

hello all R experts, how do I calculate the reliability between the two groups using the ICCs? I'll appreciate your reply, Thanks Sincerely, Supreet kaur, Biomedical research engineer, Nationwide Childrens Hospital, Columbus, OH (614)355-3509 [[alternative HTML version deleted]]

Hi all, I am using the vegan package to run a prcincipal components analysis on forest structural variables (tree density, basal area, average height, regeneration density) in R. However, I could not find out how to extract factor loadings (correlations of each variable with each pca axis), as is straightforwar in princomp. Do anyone know how to do that? Moreover, do anyone knows

Hi Listers, I am trying to reflect a PCA biplot in the x-axis (i.e. PC1) but am not having much success. In theory I believe all I need to do is multiply the site and species scores for the PC1 by -1, which would effectively flip the biplot. I am creating a blank plot using the plot command and accessing the results from a call to rda. I then use the calls to scores to obtain separate site and

Dear friends, I am beginning to use R software in my academic research and I'm having some problems regarding the use of PCA. I have a table with 24445 rows and 9 columns, and I used the function prcomp() to do the analysis. Working with an example?: x<-read.table("test.txt", header=T) row.names(x)<-x[,1] x<-x[,-1] require(stats) pca<-prcomp(x, scale=T) names(pca) ##

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

G'Day R users! Following an ordination using prcomp, I'd like to test which variables singnificantly contribute to a principal component. There is a method suggested by Peres-Neto and al. 2003. Ecology 84:2347-2363 called "bootstrapped eigenvector". It was asked for that in this forum in January 2005 by J?r?me Lema?tre: "1) Resample 1000 times with replacement entire

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 All, Can anyone point me a hint (package) how to draw a 3D plot using the first 3 components from PCA? Thanks a lot, FD [[alternative HTML version deleted]]

Hi Im currently trying to plot my NMDS data together with fitted variables (envfit funct) on an ordination plot. The plot function shows two displays="sites" and "sp". I was wondering how to plot it so that the sites come up as different points for different sites but the species come up as actual names? It looks a little busy at the moment with everything in. Sya -- View

Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: 7bit X-Virus-Scanned: by amavisd-new X-Spam-Checker-Version: SpamAssassin 2.63 (2004-01-11) on hypatia.math.ethz.ch X-Spam-Level: **** X-Spam-Status: No, hits=4.2 required=5.0 tests=MSGID_FROM_MTA_HEADER,RCVD_IN_BL_SPAMCOP_NET autolearn=no version=2.63 Dear r-help list, my question is about ordination technics:

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070401/33921c2a/attachment.pl

Hello, R-Help - I am trying to plot the results of an ordination from package vegan. The tricky part for me right now is getting the colors of the ellipses denoting the 95% confidence intervals of the group centroids to match the colors of the points for those same groups. >From and earlier post, I saw the code to make a plot of the ordination using different colors for my different groups.

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

I have a general question about using thin plate splines in the surf and ordisurf routines. My rudimentary knowledge of a gam is that with each predictive variable there is a different smooth for each one and then they are added together with no real interaction term (because they don't handle this well?). Now, If I have two variables that have a high D^2 score and a low GCV score (I am

Hi, I have a problem with making PCA plots that are readable. I would like to set different sympols instead of the numbers of my samples or their names, that I get plotted (xlabs). How is this possible? With points, i don´t seem to get the right data plotted onto the PCA plot, as I do not quite understand from where it is taken. I dont know how to plot the correct columns of the prcomp

Hello, supposing that I have two or three clear categories for my data, lets say pet preferece across fish, cat, dog. Lets say most people rate their preference as being mostly one of the categories. I want to do pca on the data to see three 'groups' of people, one group for fish, one for cat and one for dog. I would like to see the odd person who likes both or all three in the

?for doesn't return anything help.search("for") doesn't return anything- Is the for loop so prevelant in computer programing that the documentation is implicit or is R paradigm to discourage the use of the for loop. I will post data probably tonight, but here is my problem. I have preformed an MDS on a set of data. I have the scores of the four axes that are the optimal

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