similar to: summary and min max

Hi, The following call: xyplot(incidence ~ year, melanoma, panel=panel.smooth) produces a blank plot region with an error message: Error using packet 1 plot.new has not been called yet > sessionInfo() R version 2.11.1 (2010-05-31) x86_64-pc-linux-gnu locale: [1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C LC_TIME=en_CA.UTF-8 LC_COLLATE=en_CA.UTF-8 LC_MONETARY=C

Hi, I posted in R-devel, thinking it might be a general R problem, but Uwe cannot reproduce it in his GNU/Linux distro, with the same R version. So I'm posting here to ask whether Debian R users are having this problem. I'm getting a segmentation fault as follows: ---<--------------------cut here---------------start------------------->--- R> begt <-

Hello, Suppose we need a function that takes a POSIXct object and need to calculate the time difference between it and GMT time: gmtDiff <- function(time) { time.gmt <- as.POSIXct(format(time, tz="GMT")) time.plt <- as.POSIXlt(time) dlstime <- ifelse(time.plt$isdst > 0, 1, 0) timezone <- as.numeric(difftime(time, time.gmt, units="hours"))

Hi, When using method L-BFGS-B along with a parscale argument, should the lower and upper bounds provided be on the scaled or unscaled values? Thanks. Cheers, -- Seb

Hi, I've run into a ridiculous problem I can't find any solutions for in the archives or help pages: data(barley) cutYield <- with(barley, by(yield, variety, cut, breaks = c(0, 30, 60, 90))) As in this example, I'm using 'by' to return a factor for each level of another factor. The problem is that 'by' returns a list of the factors, and I need all these factors

Hi, I'd like to add, say, the sample size for every group in a bwplot as a parenthetical annotation to the axis. Here's a sketch: --8<---------------cut here---------------start------------->8--- require(Hmisc) age <- sample(1:100, 1000, replace = TRUE) sex <- gl(2, 8, 1000, c("Male", "Female")) grp <- gl(4, 6, 1000, letters[1:4]) bwplot(grp ~ age |

Hello, I'm adding a private library tree in my home directory by adding this to my ~/.Renviron: R_LIBS="~/R/library:${R_LIBS}" so that once in R: R> cat(Sys.getenv("R_LIBS"), "\n") ~/R/library:/usr/local/lib/R/site-library:/usr/lib/R/site-library:/usr/lib/R/library Is this the best way to proceed? Cheers, -- Sebastian P. Luque

Hi, I can't pinpoint exactly since what R version I first noticed this error, but are others seeing this? R> vignette("grid") R> sh: : command not found instead of having the vignette displayed in the proper viewer? This is with: R> sessionInfo() R version 2.9.2 (2009-08-24) x86_64-pc-linux-gnu locale:

Hi, After the latest update of r-cran-* packages in sid, I'm getting the following message on startup: During startup - Warning message: package 'lattice' was built under R version 2.9.2 and help may not work correctly I haven't found anything wrong with the package's functionality so far though. ---<--------------------cut

Hi, I'm getting a segmentation fault as follows: ---<--------------------cut here---------------start------------------->--- R> begt <- as.POSIXct(strptime("10/01/2009 06:00:00", format="%d/%m/%Y %H:%M:%S"), + tz="GMT") R> tser <- seq(begt, by=5, length.out=91000) R> tser.trunc <- format(tser) Error: segfault from C

Hi, I've begun to use the lme4 package, rather than nlme, for more flexibility during modelling, and running the examples in lmer I receive this error message: ---<---------------cut here---------------start-------------->--- R> (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) Error in printMer(object) : no slot of name "status" for this object of class

Hi, Doing 'R CMD check diveMove' is now throwing this message: Data codoc mismatches from documentation object 'sealLocs': Variables in data frame 'sealLocs' Code: id.time.class.lon.lat Docs: class id lat lon time with: R> sessionInfo() R version 2.10.1 (2009-12-14) x86_64-pc-linux-gnu locale: [1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C

Hi, Using the following simple character vector representing a time series with fractional seconds: datetime <- c("20/09/2011 13:00:59.00", "20/09/2011 13:00:59.02", "20/09/2011 13:00:59.04") Conversion to POSIXct runs into problems; the second element is not interpreted correctly: ---<--------------------cut

Fellow R-helpers, Suppose we create a histogram as follows (although it could be any vector with zeroes in it): R> lenh <- hist(iris$Sepal.Length, br=seq(4, 8, 0.05)) R> lenh$counts [1] 0 0 0 0 0 1 0 3 0 1 0 4 0 2 0 5 0 6 0 10 0 9 0 4 0 [26] 1 0 6 0 7 0 6 0 8 0 7 0 3 0 6 0 6 0 4 0 9 0 7 0 5 [51] 0 2 0 8 0 3 0 4 0 1 0 1 0 3

Hi, The following produces almost exactly what I needed. The problems are that the 'panel.dotplot' call (commented) generates the error 'Error in NextMethod("[") : argument "subscripts" is missing, with no default'. The other problem is that the colors alternate between the levels of the 'site' variable, rather than 'year'. barley$yield2

Hi, I'm having a bit of trouble with 'scales="free"' in the segplot() function of latticeExtra. Say we need panels for each year, showing only those counties that are represented in each one: ---<--------------------cut here---------------start------------------->--- library(latticeExtra) data(USCancerRates) uscr.w <- subset(USCancerRates, state ==

Hi, The following code produces a plot on X11: splom(~iris[1:4], groups = Species, data = iris, panel = panel.superpose, key = list(title = "Three Varieties of Iris", columns = 3, points = list(pch = super.sym$pch[1:3], col = super.sym$col[1:3]), text = list(c("Setosa", "Versicolor", "Virginica")))) However,

Hi, I'm usually confused about when to use 'slot' or '@'. I've frequently read that it's always preferable to use accessor functions, so I would think the '@' operator should be avoided. However, ?slot contains the following advise: "Generally, the only reason to use the functional form rather than the simpler operator is _because_ the slot name has to

Hello, I've encountered the need to cut some chron objects of the form: R> mychron <- chron(sort(runif(10, 0, 10))) R> mychron [1] (01/01/70 16:36:20) (01/02/70 00:08:46) (01/03/70 16:54:49) [4] (01/04/70 06:45:00) (01/07/70 06:21:24) (01/07/70 18:28:44) [7] (01/08/70 00:47:05) (01/08/70 05:11:44) (01/10/70 01:07:53) [10] (01/10/70 17:46:53) into arbitrary (e.g. a given number

Hello, I'm trying to build a frequency table for a vector, broken down by the combination of factors. For further analyses, I need to have the result arranged in a new data frame, with the upper limit of the histogram's breaks, the per bin count, and the factors identifying each record. My problem is including the latter: ---<---------------cut