similar to: axis.Date problem, bug(?)

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

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

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)

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

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, My colleague has been having a problem with the following data and plotting commands. The example below is part of a larger set of plots, but I've isolated the problem to this example using this small dataset (below), which kills rgui consistently. My version info > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status

Hi, Following the Writing R Extensions manual, I created a method for the cor function. As cor is not a generic, I followed the advice of section 6.1 of the same manual and did the following: cor <- function(x, ...) UseMethod("cor") cor.default <- stats::cor cor.symcoca <- function{ some code } I used package.skeleton to create the basic set-up of my package, containing the

Dear list, I just hit an error that stopped my make && make check-devel operation on my linux box (FC3, i686 P4 2GB RAM). Just to note that I've been building the development branch(?) for some time and this is the first hint of a problem. 1) updated the src tree using svn update 2) ran ../configure --with-recommended-package=no from my build directory 3) got: R is now configured

Hi I am having a particular problem with some glm models I am running. I have been adapting code from Bill Venables 'Programmers niche' in RNews Vol 2/2 to fit ca. 1000 glm models to a combination of species 0/1 data (as Y) and related physicochemical data (X), to automate the process of fitting this many models. I have successfully managed to fit all the models and have stored the

Hi, consider the following example: I have a matrix like this: spp.mds$points [,1] [,2] CLAP0 1.79703164 -11.66716182 CLAP30 3.87034797 -7.48168377 YBI0000 10.27538316 -3.32226184 YBI0100 0.58463806 -1.25748701 hir10000 5.82907036 -4.09695960 hir10016 -2.22113472 -4.86563557 gyn30000 5.14643671 -3.58706541 gyn30032 2.99552434 -6.21815132 ...

Dear List, Consider the following example: dat <- data.frame(var1 = rnorm(100), var2 = rnorm(100), var3 = rnorm(100), var4 = rnorm(100)) oldpar <- par(mfrow = c(2,2), no.readonly = TRUE) invisible(lapply(dat, function(x) { plot(density(x), main = deparse(substitute(x))) } ) )

Dear List, Can I build a binary package (.zip) for Windows on my Linux machine from my package sources? There is no C, C++, Fortran code involved, just plain ol' R. I read through the article by Jun and Rossini, but (on first reading) this seems more targeted at building a Windows version of R and Windows package binaries that contain C, C++, Fortran code that needs to be compiled. Thanks,

Dear list, In using the scatter.smooth() function (modreg) on a small data set (100 obs) the following error was produced: > scatter.smooth(Na, S) Warning message: k-d tree limited by memory. ncmax= 200 I haven't used scatter.smooth much but when I have, I haven't seen this message before. gc() returns > gc() used (Mb) gc trigger (Mb) Ncells 417693 11.2 667722

Dear List, I am thinking about ways in which I might analyse some stratigraphic data. The nature of the data series I have generates a number of issues: 1) The data I have in mind come from a sediment core sequence taken from the bottom of a lake. The sequence is sliced into a priori defined slices, in this case 0.2cm per slice. in this way a sequence of 0.2cm slices is produced for the

Hi, I've had a look the bug list and searched though the R documentation, email lists etc. but didn't see anything on this: when I do: summary(species.glm1, correlation = TRUE) I get a correlation matrix like this: Correlation of Coefficients: ( p I(H C pH * 1 I(pH^2) * B 1 Ca . . 1 I(Ca^2) . . B attr(,"legend") [1] 0 ` ' 0.3 `.'

Hi, On a new Dell laptop, with a fresh FC3 installation (with latest updates applied) I am experiencing make hanging consistently after/during building grDevices. This happens when using the configure flag --with-lapack (I have the LAPACK rpm distributed with FC3 installed). The last two lines of output are: make[5]: Leaving directory

1) As long as the developer version of tcltk is installed, I now have no problem compiling R with tcltk support on Fedora Core 3. After extracting the gunzipped directory, I switch to the new directory and execute the following two commands: ./configure ----enable-R-shlib make 2) After compiling JGR as specified on their web site, the JGR program is initialized with the 'run' file.