similar to: what is the function in R for MARS and PRIM?

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

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

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

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

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,

Hi, I was wondering whether a package that can perform dynamic linear models on times series data was available for R? Many Thanks, Gavin Simpson %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. & ECRC [E] gavin.simpson at

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

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, I noticed the following problem with xy.coords() in R 2.2.1-patched (version info at the foot of this email) and R 2.3.0 unstable (subversion no: r37123): > xy.coords(x = matrix(1:20, ncol = 2)) Error in xy.coords(x = matrix(1:20, ncol = 2)) : argument "y" is missing, with no default > xy.coords(x = matrix(1:20, ncol = 2), y = NULL) $x [1] 1 2 3 4 5 6 7 8 9

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

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

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>

I'm using R 2.1.1 and just successfully installed packages tseries, fseries. I try to run example http://www.maths.lth.se/help/R/.R/library/tseries/html/garch.html But it shows > x.arch <- garch(x, order = c(0,2)) # Fit ARCH(2) Error: couldn't find function "garch" Then I run command > help.search("garch") it shows the R information.