similar to: Displaying a PLS biplot under a 3d surface plot

Hi list, I'm struggling with plotting a 3d data set. I have a data frame with x,y and z values. They occur in no particular order (ie x and y are not ordered). However, I want to visualize the data using a contour plot. Now, the x,y and z data going into contour() are supposed to be ordered. My question is: Is their a quick way to re-arrange my data such that it can be fed to contour()?

Hello, I'd like to reproduce the standard biplot with pc.biplot=TRUE using xyplot in lattice in order to assign different symbols to groupings of observations (similar to the example in fig. 11.2 on page 285 in DAAG). In order to reproduce the biplots I need to know how to scale the observations and variables. In ?biplot.princomp I found that the variables are scaled by 'lambda ^

Hi, I think that the documentation for the biplot function `biplot.princomp' is inconsistent with what it actually does. Here is what the documentation states: pc.biplot If true, use what Gabriel (1971) refers to as a "principal component biplot", with lambda = 1 and observations scaled up by sqrt(n) and variables scaled down by sqrt(n). Then inner products between

> biplot standardGeneric for "biplot" defined from package "stats" function (x, ...) standardGeneric("biplot") <environment: 0x0d4444d8> Methods may be defined for arguments: x Use showMethods("biplot") for currently available ones. > > > showMethods("biplot") Function: biplot (package stats) x="ANY"

Hi, I am trying to understand how the biplot.prcomp is constructed so I can manipulate it to emphasise particular observations and reduce the number of variables shown. The prcomp model I have ran has cor=TRUE and scale=TRUE I have worked out from looking at str(prcomp.model) that... prcomp.model$x = the observations ploted in the biplot prcomp.model$rotation = the variables that form the

I'm attempting to modify how biplot draws its red vectors (among other things). This is how I've started: Biplot <- function(xx, comps = c(1, 2), cex = c(.6, .4)) { ## Purpose: Makes a biplot with princomp() object to not show arrows ## ---------------------------------------------------------------------- ## Arguments: xx is an object made using princomp() ##

[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

Dear R Community, We am pleased to announce the release of the 'bcpa' package. The source code and binaries is now on CRAN (http://cran.r-project.org/web/packages/bpca/index.html). The bpca package implements biplot (2d and 3d) and diagnostic tools of the quality of the reduction. It depends R (>= 2.6.0), scatterplot3d, rgl and MASS packages. Regards, --

Hello all, I just started in on R today. I am at the point where I am trying to distinguish different groups of data (hydrochemical data) in biplot. I search and found a post about converting the sample number in biplot to a symbol like x or + "> temp <- matrix(runif(50), nrow=10) > temp.pca <- princomp(temp) > biplot(temp.pca, xlabs=c("A", "A",

I have 2 questions about the Biplot function: Both of the questions refer to the following type of graph, which is a biplot of a principal component analysis: biplot(prcomp(dataset)) 1. Does anyone know how to change the appearance of data points on the biplot? As it is currently, for this type of graph, the label for each data point shows up on the graph - does anyone know how to get rid

Dear all, I've been learning biplot (Gabriel, 1971) and some days ago I sent for this list a procedural function with invitation for a collaborative package. Jari Oksanen made some suggestions and I agree with all. So, I reworked the function under the object-oriented programming (OOP/S3). I think it is now a good frame for more resources. Below it is the function and a small script to

Hi! Is there a way to get biplot and identify to work togheter. Having the output of prcomp I would like to draw a biplot that. Instead of plotting the sample (row-names) names plots some "pch" symbols. (thats easy with xlabs) But now I would like to add using identify the names to only some of the points. I have noticed that both biplot.prcomp and biplot.default does a lot of

Hi everybody, Excuse me for this silly question, but after searching the help archives I'm still unable to find my way to read the source code of the "biplot" function in R. I have installed the mvbutils package, and tried: > fixr(biplot) which only gives me: function (x, ...) UseMethod("biplot") with no further details.... I'd like to read the code to see if

I am doing a principle component analysis on a dataset with a lot of different variables and have constructed a biplot of the data. Unfortunately, as can be seen on the attached image, the biplot is very messy, cluttered, and hard to read. I have performed a few modifications including outlier removal from a few of the variables, which has made the plot better, however it still is too cluttered

Dear R-helpers, The vectors in my biplot are completely obscured by the ~1400 labels R is printing on my biplot. I don't really care about the labels. How can I make the biplot without the annoying labels? See attached, if that helps you see my problem. Many thanks, Mark Na

Dear list, I want to perform a biplot, using customized titels for the x and y axis. Setting xlab="" and ylab="" resulted in an error, e.g.: > data(USArrests) > biplot(princomp(USArrests),xlab="",ylab="") Error in biplot.default(t(t(scores[, choices])/lam), t(t(x$loadings[, : length of dimnames[1] not equal to array extent > How do I

Hi everyone. I performing a simple PCA using the princomp function. Then, I use the biplot function to show it. However, the function use line number to represent samples. I would like to know if there's a way to use a dot (point) instead of the line number when using the biplot function. With regards, Phil -- View this message in context:

Dear list, I'm trying to make a biplot, but instead of plotting the row number for each observation, plot a group factor. Example: prcomp(iris[,1:4]) -> PCA biplot(PCA) #this makes a nice biplot but with row names Instead of row numbers I want to plot iris[,5], which is a factor. I can do this: plot(PCA$x[,1], PCA$x[,2], type="n") text(PCA$x[,1], PCA$x[,2], labels=iris[,5],

I am trying to do a surface profile plot. data is X                  Y(1)             Z(1) 1-jan-02       2002        number 2-jan-02       2002        number . . . 1-jan-03       2003 (Y2)     number Z(2) 2-jan-03       2003 (Y2)     number Z(2) . . . until dec 31 2007.   I used the plot3d funtions to build a scatter point plot. Call rinterface.rrun("library(rgl)") Call

Thanks Andris, Michael and Petr for your prompt and kind feedbacks. I will try generating my own biplot from low-level graph commands... I hope it will work. Best regards, Marco -- Marco Manca, MD University of Maastricht Faculty of Health, Medicine and Life Sciences (FHML) Cardiovascular Research Institute (CARIM) PO Box 616 6200 MD Maastricht E-mail: m.manca at path.unimaas.nl Office