similar to: What's the best way to tell a function about relevant fields in data frames

I have a data frame called s3. This data frame has a column called saccade which has two levels 1 and -1. > head(s3$saccade, 100) [1] NA NA NA NA -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 [26] -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 1 1 1 [51] 1 1 1 1 1 1 1 1 1 1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 [76] -1 -1 -1 -1 -1

Dear list! > gregexpr("a+(b+)", "abcdaabbc") [[1]] [1] 1 5 attr(,"match.length") [1] 2 4 What I want is the offsets of the matches for the group (b+), i.e. 2 and 7, not the offsets of the complete matches. Is there a way in R to get that? I know about gsubgn and strapply, but they only give me the strings matched by groups not their offsets. I could write

I'd like do a simple xyplot with customized order of panels and try to understand how to use index.cond for that. Several attempts didn't deliver the correct results. Now, I noticed the following: > p <- xyplot(dur~roi|trial, data) > p$index.cond [[1]] [1] 1 2 3 4 5 6 7 8 9 10 These numbers are "valid indexing vector for the integer vector

Hi, I have a data frame and want to merge adjacent rows if some condition is met. There's an obvious solution using a loop but it is prohibitively slow because my data frame is large. Is there an efficient canonical solution for that? > head(d) rt dur tid mood roi x 55 5523 200 4 subj 9 5 56 5523 52 4 subj 7 31 57 5523 209 4 subj 4 9 58 5523 188 4 subj 4 7

Hi list! I have a large matrix which I'd like to partition into blocks and for each block I'd like to compute the mean. Following a example where each letter marks a block of the partition: a a a d g g a a a d g g a a a d g g b b b e h h b b b e h h c c c f i i I'm only interested in the resulting matrix of means. How can this be done efficiently?

Hi list! An operation that I often need is splicing two vectors: > splice(1:3, 4:6) [1] 1 4 2 5 3 6 For numeric vectors I use this hack: splice <- function(x, y) { xy <- cbind(x, y) xy <- t(xy) dim(xy) <- length(x) * 2 return(xy) } So far, so good (?). But I also need splicing for factors and I tried this: splice <- function(x, y) { xy <-

I have a data frame that is the result of a cast (reshape) operation. I deleted the variable column and tried to melt the resulting data frame. Depending on which method I use to delete the column I get different error messages when melting: > head(tinfos) vpn group trial_no item relation trial_type rt variable # 1 102 2 1 4351 diag1 distractor 8471 fix_d 27

Hi list, I'm using the package XML to create a simple XML document. Unfortunately constructing the XML tree is extremely slow. My code (see below) adds only about 100 nodes per second on an Intel i5 machine. There's clearly something wrong but I don't see what. Here's a sample of the XML document: <?xml version="1.0" encoding="utf-8"?>

Hi list, I run R on Linux and OSX. On both systems I use R version 2.9.2 (2009-08-24) and reshape version: 0.8.2 (2008-11-04). When I do a melt with na.rm=T on a data frame I get different results on these systems: library(reshape) x <- read.table(textConnection("char trial wn p E10I13D0 4 r E10I13D0 4 a E10I13D0 4 c E10I13D0 4 t E10I13D0 4 i E10I13D0 4 c E10I13D0 4 e E10I13D0

When loading mclust, it shows a license agreement message. This message shows up in my document when I use Sweave. I did the following: <<echo=FALSE,include=FALSE>>= library(mclust) @ Is this a problem with mclust, Sweave or with me? How can it be fixed? Thanks for any suggestions! Titus

Dear list! I have a script that processes a large number of data files. When one file fails to process correctly, I want the script to write a message and to continue with the next file. I achieved this with tryCatch: for (f in files) tryCatch({heavy.lifting(f)}, error=function(e) log.error.to.file(e)) I also want to log warning messages and tried something like this: for (f in

Hi list! I have a data frame called fix and a list of index vectors called rois: > head(rois, 3) [[1]] [1] 2 1 [[2]] [1] 3 [[3]] [1] 6 7 28 26 27 24 25 The part that's causing the issue is the following line: lapply(rois, function(roi) fix$x[roi] <- 100) So for every index vector I'd like to set the respective entries in the data frame (fix) to 100. I

Hi list, I want to place a plot (eg "plot(1:10)") and a xyplot side by side in one figure. The two columns should have different widths. If I'm not missing something, this rules out the standard way: par(mfcol=n). I tried layout(), but it doesn't seem to like it when grid scribbles into its views (the plot works but the view for the xyplot is left empty). I would

Hello, I have run an eye-tracking experiment for which I now like to analyse the saccades. Participants looked from a fixation cross (ia = 5) to the target area (ia = 4) in following example of a data frame. ia = 9 stands for everything else. A saccade is indicated by saccade = 1. Sometimes the saccade just ends before the target area (see below). This is due to the parameters that determine a

Thomas, unfortunately it does slow down the machine. I am working in a chroot enviroment and there are other users on the same box. Therefore it's of special importance that I don't consume the cpu by myself. Titus Am Sonntag den, 15. September 2002, um 23:37, schrieb Thomas Vander Stichele: > Apart from what probably is a bad assumption in the code, does it > actually

linux-2.6 9c501935a3cdcf6b1d35aaee3aa11c7a7051a305 cleaned <linux/socket.h> from allmost any userspace export. thus define the stuff that used to be there in klibc socket.h Signed-off-by: maximilian attems <max at stro.at> --- usr/include/sys/socket.h | 37 +++++++++++++++++++++++++++++++++++++ 1 files changed, 37 insertions(+), 0 deletions(-) diff --git

Well, I figured out how to make the linker to the right thing with regards to producing a dynamic section for the core... it involves using the options "-pie -Bsymbolic" ... which means we're building a relocatable binary even though we don't want to relocate it, and -Bsymbolic is there to force the linker to resolve all symbols at link time. -hpa

This is OT but relevant to the 'tech depth' of this list. The Always-On Top 100 of 2006 were announced this evening a streaming webcast of some of the CEO's will be online tomorrow at http://www.alwayson-network.com/comments.php?id=14349_0_1_0_C If you want to see a full list of the companies you can see it here

Our asterisk system gains access to the PSTN through a voip provider. We have no digium or other telephony hardware in our system. Do the zttest results still matter to us? Our results were as follows: --- Results after 1007 passes --- Best: 100.000000 -- Worst: 99.780273 -- Average: 99.975763

Patch "ltc/gm107: use nvkm_mask to set cbc_ctrl1" sets the 3rd bit of the CTRL1 register instead of writing it entirely in gm107_ltc_cbc_clear(). As a counterpart, gm107_ltc_cbc_wait() must also be modified to wait on that single bit only, otherwise a timeout may occur if some other bit of that register is set. This happened at least on GM206 when running glmark2-drm. While we are at