similar to: plotting other axes for PCA

Is it possible to make disappear the id numbers from scatter.dudi (mc analysis) ? a <- as.factor(c(1, 2, 3, 2, 1)) b <- as.factor(c(3, 2, 3, 1, 1)) x <- as.factor(c(1, 2, 2, 1, 3)) y <- as.factor(c(2, 2, 3, 1, 1)) dat <- data.frame(a=a, b=b,x=x,y=y) summary(dat) dat require(ade4) dat.acm <- dudi.acm(dat, scann = FALSE, nf = 2) scatter.dudi(dat.acm) Thank you very much !

Hi everyone, I am doing PCA with labdsv package. I was trying to create a biplot graphs in order to observe arrows related to my variables. However when I run the script for this graph, the console just keep saying: *Error in nrow(y) : element 1 is empty; the part of the args list of 'dim' being evaluated was: (x)* could please someone tell me what this means? what i am doing

Hi everybody,   I am doing principal component analysis (PCA) using "prcomp' function. When i did "Biplot", i did not found interesting result and it is based on Principal component (PC) 1 vs PC2. Now, i want to see"Biplot" in combination of either PC1 vs PC3 or PC2 vs PC 3. I did not get the ideas. Does any one have ideas ? I am optimistic on getting some idea.

Dears, I've been learning biplot (Gabriel, 1971) and I found the function 'biplot', inside of the package 'stats', useful but, a bit limited. So, I'm thinking to start a colaborative package to enhance this methods to other multivariate methods. In this way, I would like to start it, making public a new function (biplot.pca, still in development, but running) that make

Is there an option to draw concentration ellipses in biplots ? It seems really nice to summarize large number of points of each group. Cheers../ Murli

He everybody, Well I have a biplot CCA-like origined from plot.cca (vegan package). I need to rotate on y axis the lines and symbols of constrained and sites representation. If I do that on an image editor, I rotate everything, including titles, axes labels and positions. I just need to rotate the inner par and keep the variables names (constrained) and symbols in the new positions but with the

Hi all, I've been trying to figure out how to superimpose a histogram on a biplot that shows the relative contribution of each axis. I have been using the NIPALS function (http://biomserv.univ-lyon1.fr/~dray/ files/softwares/nipals.R) to run principal component analyses. Here is a toy example. source("http://biomserv.univ-lyon1.fr/~dray/files/softwares/nipals.R")

Hello I am trying to write a function that would plot timeseries easily... I aim at plotting two time-series on 2 different y axis sharing the same x-axis I succeded in doing so: plotTimeSerie(x,y,y2,[a lot of other args]){ ... plot() axis.POSIXct(side=1) #here I build the x-axis points() #here I plot the first time serie related to y-axis ... axis(side=2,[some args])

Hi all I have produced a biplot for a PCA (see attached pdf) that I ran however the names of the variables which are placed at the end of the arrows overlap and are thus unreadable. Similarly some of the numbered points overlap. I was wondering if there was a way to edit the biplot to move the label names and if not what the best alternative is. Thanks Anna pca<-biodata[,3:10]

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

Hello, Is there a simple way of stacking/merging two dataframes in R? I want to stack them piece-wise, not simply add one whole dataframe to the bottom of the other. I want to create as follows: x.frame: aX1 bX1 cX1 ... zX1 aX2 bX2 cX2 ... zX2 ... ... ... ... ... aX99 bX99 cX99 ... zX99 y.frame: aY1 bY1 cY1 ... zY1 aY2 bY2 cY2 ... zY2 ... ... ... ... ... aY99 bY99 cY99 ...

I've tried a few things both with prcomp(), and rda() and its friends in vegan (including biplot.rda and ordiplot), but can't find a solution. I'd like to associate subsets of the points in a resulting biplot ("sites" in the rda object) with different plotting colors/text styles to emphasize certain sets of points. I can't figure out how to keep the arrows (for

[was sent a wrong R-help address; manually resent by MM] Hello I'am using the 'biplot' and 'biplot.pincomp' functions of the 'mva' package for my studies. The biplot represents both the observations and the variables of a matrix of multivariate data on the same plot. The observations are represented by their numbers (the line of the data matrix), but I would need to

x=rmvnorm(2000, rep(0, 6), diag(c(5, rep(1,5)))) x=scale(x, center=T, scale=F) pc <- princomp(x) biplot(pc) There are a bunch of red arrows plotted, what do they mean? I knew that the first arrow labelled with "Var1" should be pointing the most varying direction of the data-set (if we think them as 2000 data points, each being a vector of size 6). I also read from

Hi, I have a large data set that I input as a matrix, where I have 1:x rows with names AX, then x+1: y rows named AY, etc. The idea is that I have to count how many rows exactly I have with name AX and how many I have with name AY (or find which row numbers have names AX). Is there any way in R to count a number of rows with a name matching a required pattern? Thanks in advance, Katie --

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

The script below successfully produces a biplot of the data but the 'site names' (rows) and the names of the 'response variables' (columns) are shown as simple numerals (rather than the column and row names). How might I 'enforce' the use of the row/column names used in the datafile (section of datafile shown below)? Can anyone help, please? Section of datafile sample a b

Hi, I am doing PCA on several columns of data in a data.frame. I am interested in particular rows of data which may have a particular combination of 'types' of column values (without any pre-conception of what they may be). I do the following... # My data table. allDat <- read.table("big_select_thresh_5", header=1) # Where some rows look like this... # PDB SUNID1

Hi- I'm trying to make my PCA (princomp) colored. In my csv excel sheet, I have the first column numbered according to the groupings I want to assign to the PCA. I've played around with trying to set this first column as the color vector, but haven't had any luck. Any suggestions? Thanks, Hillary [[alternative HTML version deleted]]

I have a decent sized matrix (36 x 11,000) that I have preformed a PCA on with prcomp(), but due to the large number of variables I can't plot the result with biplot(). How else can I plot the PCA output? I tried posting this before, but got no responses so I'm trying again. Surely this is a common problem, but I can't find a solution with google? The University of Dundee is a