similar to: Applying functions to partitions

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

Hi list, I have a function that detects saccadic eye movements in a time series of eye positions sampled at a rate of 250Hz. This function needs three vectors: x-coordinate, y-coordinate, trial-id. This information is usually contained in a data frame that also has some other fields. The names of the fields are not standardized. > head(eyemovements) time x y trial 51

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

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

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

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

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

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

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

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

UNOFFICIAL Firstly thanks for the help with my previous problem building SAMBA. The UNOFFICIAL in the subject heading is added automatically by our email system. I'm getting the following error when trying to join a 2003 server domain. ... Adding CN=TITUS,OU=Domain Controllers,DC=SSUNIT050,DC=local Adding

dear helplist, my student has fifty trees, numbered one to fifty, and a vector recording which tree a certain possum slept in on 12 nights. R> c [1] 3 14 17 22 26 26 17 40 43 25 46 46 R> Thus it slept in tree #3 on Monday, then tree #14 on Tues, and so on. I wish to test the null hypothesis that the animal chooses trees randomly; try R> table(c) c 3 14 17 22 25 26 40 43 46 1 1

Hi My config: Samba 2.2.3 & Linux Debian Testing kernel 2.4.20 There are two network cards on my pc : -one links to the network R1 (138.231.121.12) -one links to the network R2 (192.168.0.1) I would like to open two shares: -the first one on R1 in the workgroup wkgrp1 -the second one on R2 in the workgroup wkgrp2 I read that, thanks to the directives "netbios

Hi all, I'm unable to compile R under an HPC system (which like many HPC systems doesn't always have the latest and greatest software). The version of GCC is 4.1.2 under CentOS 5.7. Running ./configure gives the following output: R is now configured for x86_64-unknown-linux-gnu Source directory: . Installation directory: /home/tbarik C compiler: gcc

had to fix the subject... was getting under my skin! yeah, could you put up the FLAC version of the worst track that is less than 20 megs compressed? (I'll have to grab it with a 56k modem). by worst I mean the one where shorten beats flac by the most. also: 1. what version of shorten are you using? 2. what command-line options for flac and shorten did you use on this track? thanks, Josh

Estimado, te agradeceria me envies el tutorial en espaniol para R. Estoy dando mis primeros pasos con esta aplicacion. Desde ya muchas gracias Danilo Ceschin Ph.D IGBMC 1 rue Laurent Fries 67404 ILLKIRCH CEDEX - FRANCE tel 33 3 88 65 3457 email danilo at titus.u-strasbg.fr