[R] extracting row means from a list
Andy Bunn
abunn at whrc.org
Tue Jul 19 17:13:55 CEST 2005
Hello: I'm reading in a series of text files (100 files that are each 2000
rows by 6 columns). I wish to combine the columns (6) of each file (100) and
get the row mean. I'd like to end up with a data.frame of 2000 rows by 6
columns.
foo <- list()
for(i in 1:10){
# The real data are read in from a series of numbered text files
foo[[i]] <- data.frame(x1 = rnorm(100), x2 = rnorm(100), x3 =
rnorm(100),
x4 = rnorm(100), x5 = rnorm(100), x6 =
rnorm(100))
}
str(foo)
# by hand
mean.x1 <-
rowMeans(cbind(foo[[1]][,1],foo[[2]][,1],foo[[3]][,1],foo[[4]][,1],foo[[5]][
,1]),
foo[[6]][,1],foo[[7]][,1],foo[[8]][,1],foo[[9]][,1
],foo[[10]][,1]))
mean.x2 <-
rowMeans(cbind(foo[[1]][,2],foo[[2]][,2],foo[[3]][,2],foo[[4]][,2],foo[[5]][
,2]),
foo[[6]][,2],foo[[7]][,2],foo[[8]][,2],foo[[9]][,2
],foo[[10]][,2]))
# and so on to column 6
mean.x6 <-
rowMeans(cbind(foo[[1]][,6],foo[[2]][,6],foo[[3]][,6],foo[[4]][,6],foo[[5]][
,6]),
foo[[6]][,6],foo[[7]][,6],foo[[8]][,6],foo[[9]][,6
],foo[[10]][,6]))
I've implemented this with nested loops that create temporary variables and
calc the mean, but the approach is clunky. E.g.,
# nested loops
for(i in 1:ncol(foo[[1]])){
for(j in 1:length(foo)){
# etc ...
}
}
Is there a way to build a better mouse trap?
TIA, Andy
Thanks, Andy
More information about the R-help
mailing list