search for: rois

Hi everybody, I'm having problems with the barchart command in the lattice package. I'm creating barcharts from matrices with with anything from 20 to 71 columns. When I leave the column names alone, that is they are set in the read.table command (and inherited by subsequent commands) the correct number barcharts is created by the "print(barchart(...))" command.

Hello- I have a script which steps through a series of subjects, and for the subjects I remove outlying values. After removing these outliers, I specify a cutoff, keeping only values over a certain value (e.g., 1.96). I want to populate a matrix with a statistic of the values that make the cutoff (for example, the mean). However, in some subjects, after outliers and the cutoff are specified,

Hi, I'm currently working on porting some SAS scripts to R, and hence need to do the same calculation (and get the same results) as SAS in order to make the transition easier for users of the script. In the script, I'm dealing with a two-factorial repeated-measures anova. I'll try to give you a short overview of the setup: - two between-cell factors: facBetweenROI (numbering

Hallo! I have run a EEG experiment and got the following data group: 1 2 1 2 1 2 1 2 ... as factor, 2 levels between subjects fixed effect (patient vs control) subj: 1 2 ... 14 1 2 ... 14 as factor 7 patients 7 control random effect condition: 1 1 ... 2 2 ... 1 1 ... 2 2 as factor, 2 levels within subjects, ie every subject worked on every cond fixed effect roi: 1 ... 2 ... 3 ... 4 ... as factor,

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 all, I wonder if you could help me with (what is surely) a simple R question? I've written a simple R script (attached) to perform multiple ANOVAs on the columns in a table loaded in from a a file (also attached). I loop through the list of columns for which I want to perform the ANOVAs, storing the column name in a variable. When I try to perform the ANOVA using the variable name as a

How do we specify a contrast interaction matrix for an ANOVA model? We have a two-factor, repeated measures design, with Cue Direction (2) x Brain Hemisphere(2) Each of these has 2 levels, 'left' and 'right', so it's a simple 2x2 design matrix. We have 8 subjects in each cell (a balanced design) and we want to specify the interaction contrast so that: CueLeft>CueRght

Hi, I'm still trying to figure out that GLM procedure in SAS. Let's start with the simple example: PROC GLM; MODEL col1 col3 col5 col7 col9 col11 col13 col15 col17 col19 col21 col23 =/nouni; repeated roi 6, ord 2/nom mean; TITLE 'ABDERUS lat ACC 300-500'; That's the same setup that I had in my last email. I have three factors: facSubj,facCond and facRoi. I had this pretty

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 (fi...

Hi all- I am stumped. The code in A. returns errors at lines 14 and 15 and fails to load series1 and series2. However, in B., if temp1 and temp2 are called again (which returns a "Table exists" error; see lines 14-17 in B.) series1 and series2 load correctly. Any ideas? Also-I am open to any suggestions to improve the code as I am a horrific programmer. Thanks A. 1 >

Dear Community, I have a problem I am trying to code in R. Suppose there are 7 products and each have a ROI value attached to it. There are 400 representatives who calls 150,000 customers with these product detailing and achieve sales through the calls. There is a cost per call and revenue earned from the call. Constraint is revenue should be greater than the cost. Now there can be at most 7

I'v set the default_mail_env to maildir:/var/maildir/%d/%u where /var/maildir/ is the root dorectory for my maildir tree and I'm getting the folowing error in the maillog : Jan 3 15:16:58 websrv pop3(roi): unlink(/var/maildir//roi/cur/new:2,T) failed: Is a directory Notice that the %d is missing. If I change to 'default_mail_env = maildir:/var/maildir/example.tld/%u' then

Hi everybody, I'm using the gmodels package to convert human readable contrasts into the format required by R and would be grateful if someone could confirm for me whether I've got the contrasts right in the sample code below. I'm working on the assumption that the contrasts are index according to the way that levels reports them for a factor. In my case

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

Dear All, I am doing experiments in live plant tissue using a laser confocal microscope. The method is called "fluorescence recovery after photo-bleaching" (FRAP) and here follows a short summary: 1. Record/ measure fluorescence intensity in a defined, round region of interest (ROI, in this case a small spot) to determine the initial intensity value before the bleaching. This

Hi all, I've done some anovas, found some significance and proceeded to do the post hoc tests to determine the source of the significance. I used the multcmp package as suggested by Andy Liaw (thank you). I am, however, unsure how to interpret the results of the call to simint. Could someone help me with this? I've attached the code and the results it produced. There are 3 groups

...epsiHF <- (n * (k-1) * epsiGG - 2) / ((k-1) * ((n-1) - (k-1)*epsiGG)) c(epsiGG,epsiHF) } # three factors, facROI,facCond,facSubj # facROI,facCond are main effects, facSubj is "repeatedness" # G-G and H-F corrections for a main effect # we do the gghf stuff for the ROI, which means ROIs in columns, # subjects in rows mtx <- NULL for (iROI in 1:length(unique( facROI ))) { for (iSubj in 1:length(unique( facSubj ))) { mtx <- c(mtx, mean(vecData[facROI==unique(facROI)[iROI] & facSubj==unique(facSubj)[iSubj]]) ) } } mtx <- matrix(mtx...

...epsiHF <- (n * (k-1) * epsiGG - 2) / ((k-1) * ((n-1) - (k-1)*epsiGG)) c(epsiGG,epsiHF) } # three factors, facROI,facCond,facSubj # facROI,facCond are main effects, facSubj is "repeatedness" # G-G and H-F corrections for a main effect # we do the gghf stuff for the ROI, which means ROIs in columns, # subjects in rows mtx <- NULL for (iROI in 1:length(unique( facROI ))) { for (iSubj in 1:length(unique( facSubj ))) { mtx <- c(mtx, mean(vecData[facROI==unique(facROI)[iROI] & facSubj==unique(facSubj)[iSubj]]) ) } } mtx <- matrix(mtx...

Trying to install wanpipe 3.5.18. No errors during compile. But when I reach the point where wanpipe and dahdi_cfg is started, I encountered an error. Starting WAN Router... Loading WAN drivers: wanpipe done. Starting up device: wanpipe1 wanconfig: WAN device wanpipe1 driver load failed !! : ioctl(wanpipe1,ROUTER_SETUP) failed: : 22 - Invalid

So I have setup my samba box as a PDC and everything is working fine. There is one thing that keeps bugging me: When I run locallly: smbclient -L localhost -U% I get the following: ----------------------------- Domain=[MYDOMAIN] OS=[Unix] Server=[Samba 3.0.14a-2] Sharename Type Comment --------- ---- ------- IPC$ IPC IPC Service