HI, Try this: dat1<- read.csv("data2.csv",header=FALSE,sep="\t") dim(dat1) #[1] 345156 5
res<- aggregate(V5~V2+V3+V4,data=dat1,FUN=mean) dim(res) #[1] 16436 4 head(res) # V2 V3 V4 V5 #1 1961 1 1 0.51761905 #2 1962 1 1 0.04666667 #3 1963 1 1 1.14761905 #4 1964 1 1 0.05190476 #5 1965 1 1 0.51428571 #6 1966 1 1 0.50952381 tail(res) # V2 V3 V4 V5 #16431 2000 12 31 0.6428571 #16432 2001 12 31 0.1600000 #16433 2002 12 31 0.3314286 #16434 2003 12 31 0.2271429 #16435 2004 12 31 3.7423810 #16436 2005 12 31 1.1371429 A.K. ----- Original Message ----- From: "Zilefac, Elvis" <aez...@mail.usask.ca> To: arun <smartpink...@yahoo.com> Cc: Sent: Thursday, May 16, 2013 2:25 AM Subject: RE: Loop through a simulation Dear AK, Here is data from 21 simulations or so. Thanks so much. ________________________________________ From: arun [smartpink...@yahoo.com] Sent: Thursday, May 16, 2013 12:04 AM To: Zilefac, Elvis Subject: Re: Loop through a simulation Dear Zilefac, I am not able to download your dataset. It says "An error occured while downloading". Could you send a smaller dataset? For example the data between 1961-1965 or so. ________________________________ From: "Zilefac, Elvis" <aez...@mail.usask.ca> To: "smartpink...@yahoo.com" <smartpink...@yahoo.com> Sent: Thursday, May 16, 2013 1:05 AM Subject: Loop through a simulation Dear AK, Please help again. I have this data file to loop through. The columns are as follows: [simulation, year, month, day, rainfall_values] I have daily data from 1961 to 2005 . I have made 63 stochastic simulations each from 1961 to 2005. Now, I would like to obtain the values for each rainfall_value, which is the mean from 63 such values. For example, rainfall_value of row one should be the mean(simulation 1 row 1, simulation 2 row 1,...simulation 63 row 1) etc. In the end, you will have daily data from 1961 to 200, same length as simulation 1 ,for example. Thanks so much for this great help. Atem. ______________________________________________ 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.