similar to: Subsetted 1-D arrays (PR#4110)

Hello, R-devel! I posted to R-help, and (inappropriately) to R-bugs, this R/S-Plus difference: > LETTERS[c(NA,2)] in S is c("","B"), but in R is c("NA","B") Kurt Hornik <Kurt.Hornik@ci.tuwien.ac.at> wrote: > I think we do not want to change this. ... > R> is.na(LETTERS[c(NA,2)]) [1] TRUE FALSE > so we really have NA but it is

If all rows are in the same "group", rowsum() returns a vector instead of a (1xN) matrix, contrary to documentation: R> print(z <- rowsum(matrix(1:12, 3,4), rep("x",3))) [1] 6 15 24 33 R> dim(z) NULL It worked correctly in version 1.4.0 but was broken by version 1.6.1. I'm currently using 1.7.1 under Solaris 2.8. --please do not edit the information

I would like to bring to your attention the following error message which didn't appear on previous versions (long time ago?) Thanks for all your effort Karl Version 2.2.1 Patched (2006-01-21 r37153) > f <- rep(c(1,2),each=5) > x <- tapply(f,f,sum) > colnames(x) Error in dn[[2]] : subscript out of bounds ------------------------------- Karl Thomaseth, Ph.D. Research

Hi, I have what I think is an easy question. I have a data frame, called stockdata, of stock prices that looks like this: date ticker close 1 01/02/1998 GE 24.667 2 01/05/1998 GE 25.104 3 01/06/1998 GE 24.771 4 01/07/1998 GE 24.979 5 01/08/1998 GE 24.750 6 01/02/1998 HIT 71.125 7 01/05/1998 HIT 72.313

Any doubts about R's big-league status should be put to rest, now that we have a Sudoku Puzzle Solver. Take that, SAS! See package "sudoku" on CRAN. The package could really use a puzzle generator -- contributors are welcome! -- David Brahm (brahm at alum.mit.edu) [[alternative HTML version deleted]] _______________________________________________ R-packages mailing list

Any doubts about R's big-league status should be put to rest, now that we have a Sudoku Puzzle Solver. Take that, SAS! See package "sudoku" on CRAN. The package could really use a puzzle generator -- contributors are welcome! -- David Brahm (brahm at alum.mit.edu) [[alternative HTML version deleted]] _______________________________________________ R-packages mailing list

as.character.factor contains this line (where cx=levels(x)[x]): if ("NA" %in% levels(x)) cx[is.na(x)] <- "<NA>" Is it possible that this is no longer the desired behavior? These two results don't seem very consistent: > as.character(as.factor(c("AB", "CD", NA))) [1] "AB" "CD" NA > is.na(.Last.value)[3] [1] TRUE

The following code produces 6 plots on a page, but the first is distorted and different from the others: par(mfrow=c(3,2), las=2) for (i in 1:6) { frame() par(mar=c(7, 7, 1, 1)) axis(2); box(); abline(h=seq(0,1,.5), col=2:4) } The first plot's axes are mis-aligned with the plotting area implied by the box. It seems to be a result of calling par(mar) after frame(). Is this expected

Thanks to David Meyer [david.meyer@ci.tuwien.ac.at] and David Brahm [brahm@alum.mit.edu] who suggested: m[ cbind(index.i, index.j) ] <- vals This works fine for the example I gave. Unfortunately, this approach doesn't extend to using the row and column names to make assignments: > m <- matrix("",ncol=3,nrow=3) > dimnames(m) <-

R> x <- 0. + 1:8000 R> y <- outer(-x, x, pmin) Error: cannot allocate vector of size 1000000 Kb Why does R need to allocate a gigabyte to create an 8000 x 8000 matrix? It doesn't have any trouble with outer(-x, x, "+"). Thanks. -- David Brahm (brahm at alum.mit.edu) Version: platform = i686-pc-linux-gnu arch = i686 os = linux-gnu system = i686, linux-gnu status =

After upgrading to R-2.3.1 on Linux Redhat, I was suprised by this: R> x <- c(721.077, 592.291, 372.208, 381.182) R> sum(x) - 2066.758 [1] 4.547474e-13 Now I understand that floating point arithmetic is not precise, but 1) the result is exactly 0 in R-2.2.1 (patched) on the same machine, 2) .Machine$double.eps = 2.2e-16, so the error seems quite large. Also note I get the same

I sent a message to <r-announce at stat.math.ethz.ch> last Thursday ("New package: colSums"), and still haven't seen it echoed on r-help or on the web archive (in fact there is no r-announce web archive for 2002). Is something broken? Did I need to use <r-announce at lists.R-project.org> instead? -- -- David Brahm (brahm at alum.mit.edu)

Here's another curious problem with empty subsets: R> x <- as.factor(1:10) R> x[] <- NA Error in "[<-.factor"(*tmp*, , value = NA) : Argument "i" is missing, with no default Bug, or merely different from S-Plus? -- David Brahm (brahm at alum.mit.edu) Peter Dalgaard <p.dalgaard at biostat.ku.dk> in thread [Trouble with 3-dots] said:

In a clean environment under R-2.1.0 on Linux: > x <- 1:5 > x[3] <<- 9 Error: Object "x" not found Isn't that odd? (Note x <<- 9 works just fine.) Why am I doing this? Because I'm stepping through code that normally lives inside a function, where "<<-" is appropriate. -- David Brahm (brahm at alum.mit.edu)

I'd like to play with ARIMA models of stock prices, but I am a complete novice. Could some kind soul explain the relationship among packages "ts", "tseries", "dse", "dse2", and "fracdiff"? Are they 'competing' products or does one depend on another? Where would be the best place for a novice to begin? Thanks for any advice. PS. I

Under R-2.2.1, a POSIXlt date created with "strptime" has an unknown Daylight Savings Time flag: > strptime(20051208, "%Y%m%d")$isdst [1] -1 This is true on both Linux (details below) and Windows. It did not occur under R-2.1.0. Any ideas? TIA! > Sys.getenv("TZ") TZ "" Version: platform = i686-pc-linux-gnu arch = i686 os = linux-gnu

help(library) says: "`.Last.lib' is called when a package is loaded." Probably should be "...is detached." -- David Brahm (brahm@alum.mit.edu) --please do not edit the information below-- Version: platform = sparc-sun-solaris2.6 arch = sparc os = solaris2.6 system = sparc, solaris2.6 status = major = 1 minor = 3.1 year = 2001

I'm playing with predict.Arima in the 3/19/02 development snapshot of R-devel. The following produces an error message because x is not of class "ts": R> x <- rnorm(20) R> obj <- arima(x, c(2,0,0)) R> predict(obj) Error in round(x, digits) : Non-numeric argument to mathematical function Granted the documentation for arima says x should be a time-series, but

In prior versions, as recent as R-1.5.0pre (2002-04-08), Tcl/Tk support worked just fine, with "configure" (no flags) finding /usr/local/lib/tclConfig.sh and /usr/local/lib/tkConfig.sh. In Monday's official release of R-1.5.0, Tcl/Tk support now fails for me (under Solaris 2.6): ... checking for tclConfig.sh... no checking for tclConfig.sh in library (sub)directories... no

I view plots on my screen with X11(width=.455*11, height=.455*8.5), which creates a small window with the American standard aspect ratio. Under R-2.0.1 and earlier, the default fonts were a reasonable size, but under R-2.1.0 they are too big and fat. I now have to either set pointsize=10 in X11(), or par(cex=.7) afterwards. The NEWS file has this to say about X11 fonts: The changes to font