similar to: RE-Row-wise two sample T-test on subsets of a matrix

Hello all, I am trying to run a two sample t-test on a matrix which is a 196002*22 matrix. I want to run the t-test, row-wise, with the first 11 columns being a part of the first group and columns 12-22 being a part of the second group. I tried running something like (temp.matrix being my 196002*22 matrix) t.test(temp.matrix[,1:11],temp.matrix[,12:22],paired=TRUE) or somthing like

Hello all, I am trying to create a matrix of 1s and -1s without any repetitions for a specified number of columns. e.g. 1s and -1s for 3 columns can be done uniquely in 2^3 ways. -1 -1 -1 -1 -1 1 -1 1 -1 -1 1 1 1 -1 -1 1 -1 1 1 1 -1 1 1 1 and for 4 columns in 2^4 ways and so on. I finally used the function combn([0 1],3) that I found at the following link

Hello all again, I want to do bitwise addition in R. I am trying to generate a matrix 0000 0001 0010 .... .... 1111 I know the other ways of generating this matrix but I need to look at bitwise addition. Any suggestions??? thanks a lot Nameeta ------------------------------------------------- This email is intended only for the use of the individual or...{{dropped}}

Hello R group, Given a correlation matrix, I would like to obtain the best subset of pairs in the matrix of some size > n such that the mean of r for that subset is a maximum compared to any other possible subset of size > n. I've been looking at the deal and subselect packages but they don't seem to do what I need. Does anyone have any suggestions? Thanks in advance, Ryan

Hi, Having several daily wind speed time series I want to extract those consecutive days over and below certain values (i.e. 5 < x <8) Don't know which funtion to use (aggregate, lapply?) and how to do it. Thanks in advance Antonio

Greetings, I'm learning R and I'm stuck on a basic concept: how to specify a logical condition once and then perform multiple transformations under that condition. The program below is simplified to demonstrate the goal. Its results are exactly what I want, but I would like to check the logical state of gender only once and create both (or any number of) scores at once.

Hi! I've made great progress in my R programming, but I am again stuck on a beginner's problem. I would like to start R with a command line that loads a file, and if possible, executes a function. Can anyone give me an example of how to do this? For example, in lisp, I would say: $ lisp -load toto.lisp -eval '(do-something $PORT)' to load the file "toto.lisp", then

i have the following code below and at the end there are some plotting statements. it actualy looks quite nice when you run it but there is just one strange thing happening that don't know how to fix. the three things being plotted are aggfxdata[,logbidask] which has its own set of times ( it's a zoo object ). rollmeanlogbidask which has a subset of the times that aggfxdata has.

I think I asked a similar question 3 years ago to the Splus list and I think the answer was no or noone answered so noone should spend more than 5 minutes on this because it could definitely be a waste of time. My question is whether the function below can be rewritten without a for loop. apply is fine if it can be done that way but i doubt it. I call it a lot and would prefer to not loop.

A development note: In the function "install.packages", it would be helpful to those of us who have atypical installations and install manually from ZIP files to have pkgnames <- sub("\\.zip$", "", pkgnames) replaced with pkgnames <- sub("\\.zip$", "", pkgnames, ignore.case = TRUE) because the contributed zipfiles are ZIPfiles. The

A development note: In the function "install.packages", it would be helpful to those of us who have atypical installations and install manually from ZIP files to have pkgnames <- sub("\\.zip$", "", pkgnames) replaced with pkgnames <- sub("\\.zip$", "", pkgnames, ignore.case = TRUE) because the contributed zipfiles are ZIPfiles. The

I have comma delimited asci data with each row being in the format : 2006-01-24 02:41:24.00011,1.22930000,5,1.22950000,7 . . . . and i'm trying to use read.zoo ( which is similar to read.table ) to read in the data. the data goes all the way out to milliseconds and i can't figure out what to put for the format field. if i put nothing, then read.zoo gets rid of the the whole time and

I did a ?do.call but i don't think i understand it. if a, b,c,d are numeric vectors then could someone explain the difference between do.call("cbind",list(a,b,c,d)) and cbind(a,b,c,d). or point to an archive on it. the return value of cbind is a matrix or dataframe depending on what is sent in but i don't understand wheen it would be useful to use do.call. i realize it

It seems like this should be pretty straight forward, but for some reason the answer escapes me. I have a normal distribution S made up of two normal distributions C and C-bar. I need to find the percentage of the area of S that both C and C-bar occupy. Any suggestions? Thanks in advance, Alan Gibson

I have two numeric vectors each of length 17 and each is named the exact same way. so obsnum ppppp ppppm pppmp . dot dot dot...... temp1 is 1417 52 63 85 obsnum ppppp ppppm pppmp . dot dot dot...... temp2 is 1213 41 50 97 what i want to have is a resultant matrix with 2 rows and 16 columns where the 16

i am looping and creating plots which are coming to the screen. i am in linux and remember ( in a previopus life ) i used to use a command so that the next graph n + 1 didn't appear on the screen until i hit the return key after graph n appeared. i thought i remeber using the command unix but when i type unix at the r prompt, it gices me system. it's probably a deprecated command then but

as jim pointed out ( i think we were figuring this out simultaneously. thanks a lot jim ), it looks it does have something to do with the fact that it's a zoo object because below i consider two cases. in the first case, fxdatab is a zoo object and i get the length of temp to be 1. in the second case, fxdatac is a matrix and the results in temp are correct. gabor : i really hate to bother

I familarized myelf with kalmanlike and structts which are approaches for building and estimating ( and forecasting ) state space models ( or the equivalent arima models ). back in 2003, gavin simpson wrote an email describing the west and harrison apprach to estimate state space models and asked if anything was out there for using that approach. the goals of this approach are the same as kalman

hi everyone : i'm trying to rotate the x axis labels in my barlot ( yes, since yesterday ) and my prob12matrix is 2 rows and 16 columns and i followed the instructions in the archive "rotate barplot labels". i was really looking forward to seeing my labels turned but the plot just came back with no labels ? can anyone see what might be going on ? it worked when i followed the

Given: > mitest <- matrix(1:16,ncol=4) > dimnames(mitest)[[1]] <- c("a","b","c","d") > dimnames(mitest)[[2]] <- c("a","b","c","d") > mitest a b c d a 1 5 9 13 b 2 6 10 14 c 3 7 11 15 d 4 8 12 16 I can do: > mitest[cbind(c(1,2,3),c(2,3,2))] [1] 5 10 7 but using the names does not