search for: mmuurr

hi all - i'm just wondering what sort of code people write to essentially performa an aggregate call, but with different functions being applied to the various columns. for example, if i have a data frame x and would like to marginalize by a factor f for the rows, but apply mean() to col1 and median() to col2. if i wanted to apply mean() to both columns, i would call: aggregate(x, list(f),

hi folks, i use ESS mode in emacs often to interact with R, and while i know how to save a session transcript, i'm wondering how to save just the history of the commands (i.e. identically to how R gives the history save option from its native CLI). if you use ESS with an R process, commands entered from a separate buffer and sent to the process do not get registered with the R history, and

hi all - i'm having trouble using lme to specify a mixed effects model. i'm pretty sure this is quite easy for the experienced anova-er, which i unfortunately am not. i have a data frame with the following columns: col 1 : "Score1" (this is a continuous numeric measure between 0 and 1) col 2 : "Score2" (another continuous numeric measure, this time bounded between 0

i've written a function to coerce a matrix (e.g. from numeric to logical), but i'd like to know if someone has a more elegant method for this: > m <- matrix(c(0, 1, 1, 0), ncol = 2) > m <- as.logical(m) > m [1] FALSE TRUE TRUE FALSE i'd like 'm' to still be a matrix with the original dimensions. in my function to do this, i coerce 'm' to a logical,

hi all - i've seen versions of this question before, but none seem to get directly at my solving my (probably very simple) issue: i simply want to annotate the tick marks on an axis with (superscripted) 10^x notation, and tried this: axis(1, at = axTicks(1), as.expression(substitute(10^foo, list(foo = axTicks(1)))) thinking the as.expression/substitute would create the appropriate

Hi all, I have a question about finalizers... I have a package that manages state for a few connections, and I'd like to ensure that these connections are 'cleanly' closed upon either (i) R quitting or (ii) an unloading of the package. So, in a pared-down example package with a single R file, it looks something like: ##### BEGIN PACKAGE CODE ##### .CONNS <- new.env(parent =

i'm somehow embarrassed to even ask this, but is there any built-in method for doing this: my_list <- list() my_list[[1]] <- matrix(1:20, ncol = 5) my_list[[2]] <- matrix(20:1, ncol = 5) now, knowing that these matrices are identical in dimension, i'd like to unfold the list to a 2x4x5 (or some other permutation of the dim sizes) array. i know i can initialize the array, then

this may be slightly off-topic, as it doesn't pertain directly to the R application, but some of the documentation. when reading R's fullrefman.pdf (available from http://cran.r-project.org/doc/manuals/fullrefman.pdf) in Mac OS X's preview.app (version 4.1, on Mac OS 10.5.x), if i try to do a keyword search within the document, the indexing step freezes about 2/3 the way through the

hi all - i'm having some difficulty figuring out how to convert between "user units" (which i can't find a definition for in the plotrix package) and either (a) device units (e.g. inches with PDFs) or (b) user coordinates along any particular axis. as an example, suppose i set up a PDF device with inches, the device has both outer and inner magins, and the plot region has

Hi - since Rgraphviz was officially moved over to Bioconductor, this might be a misguided post, but it's worth a shot: I'm plotting a graph using Rgraphviz and in an attempt to force the edges to be behind the nodes (on a fairly complex graph), I set the "outputorder" graph attribute to "edgesfirst" (as is described in the graphviz docs here ->

hi all - i'm trying to 'R CMD build' a package, but i have what appears to be a bootstrapping problem: i've included a vignette in my package, with R code interwoven (and built using Sweave), but in this documentation i have a code line: > library(MyPackage) now, when trying to build a .tar.gz install-able version of my package in a clean state, i remove the original

i've noticed a strange problem when plotting a stepfun. according to the documentation, the xlim parameter should bound the range of the function being plotted, and is returned as the extreme two values (i.e. first and last) in the vector t from the plot.stepfun call. instead, it plots beyond the desired range (although the limits are preserved for the viewing space). to reproduce: foo

hi all - quick question: i have a matrix m, say nrow=5, ncol=4. in a function i'd like to retrieve certain rows or columns from m, but which rows/cols are not known ahead of time. the function should return a sub-matrix (i.e. still of class 'matrix'). when selecting a single column (or row), the indexing operation converts the object of class 'matrix' into a vector of

hi all - quick question: i have a matrix 'y' of response values, with two explanatory variables 'x1' and 'x2'. tested values of 'x1' and 'x2' are sitting in two vectors 'x1' and 'x2'. i want to learn model parameters without "unrolling" the matrix of response values. example below: # some fake data for the example x1 <- 1:5 x2

Hi all, this is probably some tricky configuration file problem, but I figure someone here might have come across this in the past: In short, I've been using Rscript to run my scripts, usually successfully. But recently I've run into a strange problem, and the only function that causes an error is "print". I can create a file that contains only a print command, and it causes

hi, i'm trying to prepend some plain (i.e. unevalutated) text to a typographically evaluated R expression in a legend(...) call. here's a working legend(...) call that is close to what i'd like (where x and y are returned from an lm(...) call): legend("topleft", legend = c(bquote(R^2 == .(summary(x)$r.squared)), bquote(R^2 == .(summary(y)$r.squared))), fill =

hi all -- this might not be the correct list for this question/discussion, though R-help didn't seem like the correct venue, either, so... i'm looking for just some extra clarification of how local variables are defined/bound, beyond the simple cases given in the Language document. the particular instance is when there is variable assignment inside a function. normally, this creates a

Hi all --- this is less a specific question and more general regarding S3 classes. I've noticed that quite a few very common default implementations of generic functions (e.g. `unique`, `[`, `as.data.frame`) strip away class information. In some cases, it appears conditionals have been created to re-assign the class, but only for a few special types. For example, in `unique.default`, if the

hi, i'm having some trouble getting a package to load a shared library object in .onLoad(...) i have a shared object file, say "mylib.so". if i start an R session, and via the CLI specify the actual library via: > dyn.load("mylib.so") everything works quite well (i.e. i can then follow with some .Call (...) methods) now, i'd like to include this shared library in

hi all - i'm working on an R package that makes use of my own shared library written in C. but i also am making use of another C-written library. (my package is for facilitating biological namespace translations via online (i.e. up-to-date) biological databases.) problem is, the library i'm using is not a standard library (i.e. i doubt it will be installed on most users' machines). i