Amit Patel wrote:
--- On Fri, 15/5/09, Amit Patel amitrh...@yahoo.co.uk wrote:
From: Amit Patel amitrh...@yahoo.co.uk
Subject: Help with loops
To: r-help@r-project.org
Date: Friday, 15 May, 2009, 12:17 PM
Hi
I am trying to create a loop which averages replicates in
my data.
The original data has many rows. and consists of 40 column
zz[,2:41] plus row headings in zz[,1]
I am trying to average each set of values (i.e. zz[1,2:3]
averaged and placed in average_value[1,2] and so on.
below is my script but it seems to be stuck in an endless
loop
Any suggestions??
If you want row-wise means, rowMeans does it without any loop, but
instead you could loop over the column indices or go without any
explicit loop as in:
averaged_value - zz[,1:21]
temp - as.matrix(zz[,2:41])
dim(temp) - c(nrow(temp),2,20)
# and then either:
temp - aperm(temp, c(1,3,2))
averaged_value[,2:21] - rowMeans(temp, dim=2)
#
averaged_value[,2:21] - apply(temp, c(1,3), mean)
Uwe Ligges
for (i in 1:length(zz[,1])) {
#calculates Meanss
#Sample A
average_value[i,2] - rowMeans(zz[i,2:3])
average_value[i,3] - rowMeans(zz[i,4:5])
average_value[i,4] - rowMeans(zz[i,6:7])
average_value[i,5] - rowMeans(zz[i,8:9])
average_value[i,6] - rowMeans(zz[i,10:11])
#Sample B
average_value[i,7] - rowMeans(zz[i,12:13])
average_value[i,8] - rowMeans(zz[i,14:15])
average_value[i,9] - rowMeans(zz[i,16:17])
average_value[i,10] - rowMeans(zz[i,18:19])
average_value[i,11] - rowMeans(zz[i,20:21])
#Sample C
average_value[i,12] - rowMeans(zz[i,22:23])
average_value[i,13] - rowMeans(zz[i,24:25])
average_value[i,14] - rowMeans(zz[i,26:27])
average_value[i,15] - rowMeans(zz[i,28:29])
average_value[i,16] - rowMeans(zz[i,30:31])
#Sample D
average_value[i,17] - rowMeans(zz[i,32:33])
average_value[i,18] - rowMeans(zz[i,34:35])
average_value[i,19] - rowMeans(zz[i,36:37])
average_value[i,20] - rowMeans(zz[i,38:39])
average_value[i,21] - rowMeans(zz[i,40:41])
}
thanks
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.