similar to: StructTS and arima and missing values

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 have a bit of a problem with a legend in png, eps and pdf plots produced from a custom plotting function. I was writing a little function to produce some stock plots for a routine analysis conducted in our lab. I have a wrapper function figures() (see below) that produces eps, png and pdf versions of the lab output. When I draw the plot on the X11() device the correct amount of space is

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

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

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

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.

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

Dear list, The help page for library/require contains the following paragraph in the section "Packages that require other packages": The source code for a package that requires one or more other packages should have a call to 'require', preferably near the beginning of the source, and of course before any code that uses functions, classes or methods from the

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, 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, 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'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 `.'

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, Having reset-up a print server from installation disks today, I did something stupid and blindly copied the contents of /etc/samba from the previous installation into the new installation. Now smbd refuses to start (nmbd starts OK). I am currently using a blank smb.conf file and the smb.conf file from the old installation. testparm doesn't report any errors in the smb.conf from the old

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