search for: spluque

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 came across a case where there's a discrepancy between minimum and maximum values reported by 'summary' and the 'min' and 'max' functions: ---<---------------cut here---------------start-------------->--- R> str(tt) num [1:1397] 1952 1970 1976 1967 1946 ... R> summary(tt) Min. 1st Qu. Median Mean 3rd Qu. Max. 1920 1960 1970

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

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

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, Does anybody have experience building windows packages in Debian? I think the mingw32 package is the ticket for it. In the past, I've done this using the procedure on the R News article (Building Microsoft Windows versions of R and R packages under Intel Linux), but it seems it should be possible with mingw32 directly. Any tips on what to read? Thanks. Cheers, -- Seb

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

Hi, Recently during an upgrade I was asked to use labels or UUID in /etc/fstab. I've done that and everything went fine, except for the only CD/DVD ROM device in my system. According to KInfoCenter (the KDE devide information application) the this device has the label: CD/DVD Drive. So I placed this line in /etc/fstab: LABEL="CD/DVD Drive" /media/cdrom0 udf,iso9660 user,noauto 0

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,

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'm trying to fit the following model to data using 'nls': y = alpha_1 * beta_1 * exp(-beta_1 * x) + alpha_2 * beta_2 * exp(-beta_2 * x) and the call I've been using is: nls(y ~ alpha_1 * beta_1 * exp(-beta_1 * x) + alpha_2 * beta_2 * exp(-beta_2 * x), start=list(alpha_1=4, alpha_2=2, beta_1=3.5, beta_2=2.5), trace=TRUE, control=nls.control(maxiter =

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

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

Hi, I have no problems whatsoever accessing and using my GNU Linux user directory from Windows Vista, but the reverse is not true. My /etc/samba/smb.conf has this (via testparm -s): ---<---------------cut here---------------start-------------->--- ~$ testparm -s /etc/samba/smb.conf Load smb config files from /etc/samba/smb.conf Processing section "[homes]" Processing section

Hi, How can we set the environment for the minuslog function in mle()? The call in this code fails because the "ll" function cannot find the object 'y'. Modifying from the example in ?mle: library(stats4) ll <- function(ymax=15, xhalf=6) { -sum(stats::dpois(y, lambda=ymax/(1+x/xhalf), log=TRUE)) } fit.mle <- function(FUN, x, y) { loglik.fun <- match.fun(FUN)

Dear Ubunutu R users, Vincent Goulet and myself would like to announce the availability of R 2.11.0 for Ubuntu's latest release, Lucid Lynx. The packages should be on or arriving soon to your favorite mirror. R 2.11.0 packages also exist on CRAN for all Ubuntu releases going back to hardy. As Johannes Ranke pointed out for the Debian release, R packages installed under R < 2.10.0