similar to: Regular expressions: offsets of groups

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

What causes the error report: logical(0) to arise in the rms function lrm? Here's my data: But both the dependent and the independent variable seem fine... > str(AABB) 'data.frame': 1176425 obs. of 9 variables: $ sex : int 1 1 0 1 1 0 0 0 0 0 ... $ faint : int 0 0 0 0 0 0 0 0 0 0 ... Here's the simplified model and error AABB$model1 < lrm (faint ~ sex)

Good evening, I understand that NTFS & VMDK do not relate to Solaris or ZFS, but I was wondering if anyone has any experience of checking the alignment of data blocks through that stack? I have a VMware ESX 4.0 host using storage presented over NFS from ZFS filesystems (recordsize 4KB). Within virtual machine VMDK files, I have formatted NTFS filesystems, block size 4KB. Dedup is turned on.

We are tried to count how many files belong to certain group. Our system administrator told us "non-owner" can easy change file group name to another. I have been tried several combination and never successful (only ROOT can change file group to other name). Does anyone know how "no-owner" can change file group name? Thanks.

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

Is there a cleaner way of combining two expressions. This example works and gives what I want plot(1:10) aa <- expression(alpha==.05) bb <- expression(beta ==.80) aabb <- expression(alpha==.05 ~ ", " ~ beta ==.80) text(5, 10, aa) text(5, 9, bb) text(5, 8, aabb) text(5,1, parse(text=paste(deparse(aa[[1]]), deparse(bb[[1]]), sep="~"))) text(5,2,

On 23/11/2014 09:39, peter dalgaard wrote: > >> On 23 Nov 2014, at 01:05 , Henrik Bengtsson <hb at biostat.ucsf.edu> wrote: >> >> On Sat, Nov 22, 2014 at 12:42 PM, Duncan Murdoch >> <murdoch.duncan at gmail.com> wrote: >>> On 22/11/2014, 2:59 PM, Stuart Ambler wrote: >>>> A colleague?s R program behaved differently when I ran it, and we