Dear Ben,Thanks for replying but its still not working.your code was>lapply(res,colMeans)but i want to use "res1" instead of "res". when i did use it, i got same error.eliza > To: r-h...@stat.math.ethz.ch > From: bbol...@gmail.com > Date: Mon, 24 Dec 2012 00:31:41 +0000 > Subject: Re: [R] colmeans not working > > eliza botto <eliza_botto <at> hotmail.com> writes: > > > Dear useRs,You must all the planning for the christmas, but i am > > stucked in my office on the following issue i had a file containg > > information about station name, year, month, day, and discharge > > information. i opened it by using following command > > > > dat1<-read.table("EL.csv",header=TRUE, sep=",",na.strings="NA") > > You can probably use > > dat1 <- read.csv("EL.csv") > > (although you may have to double-check some of the other > default differences between read.csv and read.table, e.g. > quote and comment.char arguments) > > > then by using following codes suggested by arun and rui i managed to obtain > > an > output > > library(reshape2) > res <- lapply(split(dat1,dat1$st), > function(x) dcast(x,month~year,mean,value.var="discharge")) > > [snip] > > res1 <- lapply(res, function(x)x[,-1]) > > (c() is redundant here) > > > $EE > > 2005 2006 2008 2009 > > 1 1.7360776 0.8095275 1.6369044 0.8195241 > > 2 0.6962079 3.8510720 0.4319758 2.3304495 > > 3 1.0423625 2.7687266 0.2904245 0.7015527 > > 4 2.4158326 1.2315324 1.4287387 1.5701019 > > > > $WW > > 2008 2009 2010 > > 1 1.4737028 2.314878 2.672661 > > 2 1.6700918 2.609722 2.112421 > > 3 3.2387775 7.305766 6.939536 > > 4 6.7063592 18.745256 13.278218 > > > > > > Now you just need > > lapply(res,colMeans) > > ______________________________________________ > 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.