similar to: Reading biplot function's source code

Hi, I wonder if it is possible with the graphics package to overlay one or several plots (barplots, for example) over a map. Data for the map is in a data frame with the latitude and longitude coordinates, and then: > plot(map$long, map$lat, type ="l") produces the map. I want to put each barplot in specific locations on the map, namely at the center of "statistical

Hi, I'm sure this one is very easy.... I am trying to write a function where one of its arguments has two posible (strings) values, defaulting to one of them if none is specified. My problem is that when evaluating the function the following warning is produced: "the condition has length > 1 and only the first element will be used" I've read the help page of if() but I

Dear list, I would like to plot arrows with different colors according to arrow length, and also (if possible) with head size proportional to arrow length. The idea is to make a quiver-like plot of matlab with wind speed data. So far, I´ve been able to use different colors, but I need to find a more efficient way to recode arrow length intervals into colors. On the contrary, I can't define

Hello, In doing a series of non-linear estimations of a function which is a sum of a varying number of sinusoids, I would like to "autogenerate" the arguments needed by nls() depending on that number. For example, when there are two sinusoids: > nls( y ~ mu + A1 * cos(2*pi*f1*x - P1) + A2 * cos(2*pi*f2*x - P2), data = some.xy.data, start = list( mu=some.value0,

Hi, I'm having difficulty to replicate in R a nested (hierarchical) anova example found in p. 308 of Zar, J.H. 1996. Bostatistical Analysis. Prentice Hall. 3rd ed. The example (15.1) is as follows: The variable is blood cholesterol concentration in women (in mg/100 ml of plasma). This variable was measured after the administration of one of three different drugs, each drug having been

*This message was transferred with a trial version of CommuniGate(tm) Pro* A simple question, perhaps to simple, but I am new to R: For a data frame with k variables and n1+n2 observations (two groups) how can i substract the group means for every variable, so both group variables become centered to mean zero? Thanks in advance for your help H?ctor

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

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() ##

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

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

[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

Hello everybody, I am afraid this might be a stupid question, but after scanning a little bit through the help archive and the internet, I didn't come up with an answer myself. I have done a principle component analysis on a data set, which yields the mean response of 14 olfactory neurons to 16 different odors. When plotting the result in a biplot, it nicely shows me the odors in

Which is the parameter that is used to decrease the size of ylabs plotted in biplot? I tried playing with cex and cex.lab I am not getting it right pc <- princomp(USArrests) biplot(pc, xlabs = rep("", nrow(USArrests)),ylabs=(colnames(USArrests))) Thanks../Murli [[alternative HTML version deleted]]

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

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

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