Jorge you are awesome, This code is very efficient. it's exactly what I wanted to do and it's very short and not complicated.
Thanks a lot Jorge Ivan Velez wrote: > > Dear Alex, > Is this what you want? > > my=read.table(textConnection(" > X8 X9 X10 X102 X110 X177 X283 X284 X286 X292 X297 X306 X308 X314 > 0 1 0 0 0 1 0 0 0 0 0 0 0 0 > 0 0 0 0 1 0 0 0 0 0 0 0 1 0 > 0 1 0 0 0 0 0 0 0 0 0 0 1 0 > 1 0 0 0 1 0 0 0 0 1 0 0 0 > 0"),header=TRUE) > closeAllConnections() > rownames(my)=paste('s',1:4,sep="") > > # Five samples > res=replicate(5,my[,sample(colnames(my),5)],simplify=FALSE) > > # mean for each sample > Means=do.call(rbind,lapply(res,function(x) rowMeans(x))) > rownames(Means)=paste('sample',1:5,sep="") > Means > > # Global mean for sample > colMeans(Means) > > > HTH, > > Jorge > > > On Thu, Oct 16, 2008 at 9:10 AM, Alex99 <[EMAIL PROTECTED]> wrote: > >> >> >> Hi everyone, >> I have dataset which I take random samples of it 5 times. each time I get >> the mean for rows for each sample. >> at the end I need to calculate the Average of all means for each sample >> and >> each row. to clear it up I give an example: >> say this is my dataset. >> X8 X9X10X12 X13 X14 X15 X16X17X18X19 X20 X21 X22 >> s1 0 1 0 0 0 1 0 0 0 0 0 0 0 0 >> s2 0 0 0 0 1 0 0 0 0 0 0 0 1 0 >> s3 0 1 0 0 0 0 0 0 0 0 0 0 1 0 >> s4 1 0 0 0 1 0 0 0 0 1 0 0 0 0 >> >> I get a random sample and mean of row for that sample(5 times , but I >> just >> put 2 of them here to give you an idea) >> >> X12 X9 X10 >> s1 0 1 0 >> s2 0 0 0 >> s3 0 1 0 >> s4 0 0 0 >> s1 s2 s3 s4 >> 0.3333333 0.0000000 0.3333333 0.0000000 >> >> X10 X18 X8 >> s1 0 0 0 >> s2 0 0 0 >> s3 0 0 0 >> s4 0 1 1 >> s1 s2 s3 s4 >> 0.0000000 0.0000000 0.0000000 0.6666667 >> >> This is the code I used: >> for(i in 1:5) >> { >> temp<-sample(A3,3, replace=F) >> Avg=rowMeans(temp) >> show(temp) >> show(Avg) >> } >> >> Now, the problem is how can I save the result for each row(s1,s2,s3,s4) >> so >> that I can get the grand average from 5 runs?I thought about using a >> vector >> in the "for" loop but it's no good, it over right so basically I only get >> the means for last sample. any idea how to do it? >> Thanks a lot >> >> -- >> View this message in context: >> http://www.nabble.com/saving-result-of-a-%22for%22-loop-tp20013519p20013519.html >> Sent from the R help mailing list archive at Nabble.com. >> >> ______________________________________________ >> 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. >> > > [[alternative HTML version deleted]] > > ______________________________________________ > 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. > > -- View this message in context: http://www.nabble.com/saving-result-of-a-%22for%22-loop-tp20013519p20015092.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.