similar to: message posting

Hi, In a recent reply I sent to the list, my initial response was to suggest, politely that by reading the posting guide and learning to use the supplied tools one can often help oneself. I then went on to say that reading the posting guide would have led you to RSiteSearch("insert query"). Luckily I went to the posting guide to check this was there before posting and much to my

Dear List, I need to add a subscript and a superscript to some of the ions in the labels on some plots. I have got to here but now I'm stuck: plot(1:10, xlab = expression(paste("nm SO"[4], " ", mu, "eq cm"^{-2}, " yr"^{-1}))) Which gives almost what I require. No matter what I tried, however, I could not get bot a sub script *and* a superscript

Hi I'm having a problem with model.frame, encapsulated in this example: y1 <- matrix(c(3,1,0,1,0,1,1,0,0,0,1,0,0,0,1,1,0,1,1,1), nrow = 5, byrow = TRUE) y1 <- as.data.frame(y1) rownames(y1) <- paste("site", 1:5, sep = "") colnames(y1) <- paste("spp", 1:4, sep = "") y1 model.frame(~ y1) Error in model.frame(formula, rownames,

Dear List, I often have to merge two or more data frames containing unique row names but with some columns (names) common to the two data frames and some columns not common. This toy example will explain the kind of setup I am talking about: mat1 <- as.data.frame(matrix(rnorm(20), nrow = 5)) mat2 <- as.data.frame(matrix(rnorm(20), nrow = 4)) rownames(mat1) <- paste("site",

Sorry to disturb you, but can you help me to solve one little problem? I want to draw a graphic after another with R but I cannot find the first one after that. Do you know the command to hold the graphic with R? I remember with Matlab you may use "hold on". Thanks. -- Baohua Yang Email:yangbaohua@gmail.com [[alternative HTML version deleted]]

Dear List, I have the following function: coca <- function(x, ...) { if(is.null(class(x))) class(x) <- data.class(x) UseMethod("coca", x) } and a default method coca.default <- function(x, y, method = c("predictive", "symmetric"), reg.method = c("simpls", "eigen"), weights = NULL,

##data y1 <- matrix(c(3,1,0,1,0,1,1,0,0,0,1,0,0,0,1,1,0,1,1,1), nrow = 5, byrow = TRUE) y2 <- matrix(c (3,0,10,3,3,0,0,1,1,0,0,0,0,0,1,0,1,0,0,2,1,0,1,1,0,2,1,1,4,1), nrow = 5, byrow = TRUE) y1 <- as.data.frame(y1) y2 <- as.data.frame(y2) rownames(y1) <- rownames(y2) <- paste("site", 1:5, sep = "") colnames(y1) <-

Dear List, I have the following problem with axis.Date, here is an artificial example: dates <- scan(what = "character") "25/03/2000" "26/03/2000" "27/03/2000" "28/03/2000" "29/03/2000" "30/03/2000" "31/03/2000" "01/04/2000" "02/04/2000" "03/04/2000" "04/04/2000"

Dear list, given this formula: > fmla <- formula(y1 ~ spp1 + spp2 + spp3 + spp5) > fmla[[3]] spp1 + spp2 + spp3 + spp5 is this the intended behaviour of as.character: > as.character(fmla[[3]]) [1] "+" "spp1 + spp2 + spp3" "spp5" ? Where does the extra "+" come from? > as.character(fmla) [1] "~"

something like "sprintf" in C? so I can do: print(sprintf("the correct result is %3.4f\n", myresult)); ------- Also, I am desperately looking for a "clear console screen" function in R... thanks a lot! [[alternative HTML version deleted]]

Hi In trying to fit a linear model , I use the leaps() function to determine wich predictors I should include in my model. I would like to plot the Mallow's Cp criteria against p with the indexes of selected model variates as points labels Is there already such a function? (I could not find it) Thanks Anne [[alternative HTML version deleted]]

hello I read "Practical Time Series" (Gareth Janacek; 2001) and they presented e.g the smoothing functions msmooth(x,k) or the bivariate function crosscorr(x,y,k), but both didn't work on my machine. I only load the ts library, is another library necessary or did this function change since 2001? Is there a more recent and detailed manual for ts? thanks, cheers Martin -- Martin

where is the cook-distance in R? i can't find it!please help me!

I am attempting to use R through WinEdt. I installed the SWinRegistry package and it seems to install fine. But when I type: > library(RWinEdt) I get: Error in library(WinEdt) : there is no package called 'RWinEdt' I am running Windows XP and R 1.8.1 The SWinRegistry is correctly named as well that is it is named: SWinRegistry_0.3-3.zip, which I know caused

Hello, I suppose this is a basic question but couldn't find a solution.: I have a large matrix with let say 3 columns: V1 V2 V3 a x 2 a x 4 a y 8 b z 16 and I want to compute some statistics based on the levels resulting form the combination of the two first columns e.g.: SUM-> V1 V2 V3 a x 6 a y 8 b z 16 Thanks for your hints . Marc

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

how can i draw a 3d scatter plot? thank you.

Dear list, I am using Sweave to build a small report. I want to produce a series of figures, each figure containing a number of plots and then have them included in the Sweave file. An example would be to : postscript(file = "ANCbwplot%03d.eps", onefile = FALSE, other options...) oldpar <- par(mfrow = c(2,2)) .... do lots of plots to produce a number of eps files .... par(oldpar)

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 have seen, several times, dots (like this: "y ~." ) in formula descriptions, noticeably in R help. I am unable to see what it does correspond to. Any ideas ? -- --~~ Toulouse, Grenoble, Auch, Arcachon, B??ziers, Paris, Saragosse, L??vignac Sur Save, habitats naturel du Valdo. ~~-- < http://www.le-valdo.com>