similar to: problem using evaluating a formula

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

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

Just wondering if this is the expected behaviour. I was wanting to produce quoted text within \code{}, without manually entering the '"'. \dQuote{} seems advisable after reading the Writing R Extensions manual, so I tried \code{\dQuote{mytext}} expecting it to produce "mytext" in monospace font (with ' ' round it in the R help files) but it appears that

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

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

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

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

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, 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, A few weeks ago a discussion took place regarding Fedora Core 4 and compiling R on that platform using the new version 4 of gcc and its gfortran compiler. gcc was recently updated to 4.0.1 on FC4 (4.0.1-5 in Red Hat Land) so I thought I'd give compiling R a go on my laptop which needed updating anyway. I've had trouble using the optimisation flags I used to use under FC3 with

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.

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