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" <[email protected]>
To: arun <[email protected]>
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 [[email protected]]
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" <[email protected]>
To: "[email protected]" <[email protected]>
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.
______________________________________________
[email protected] 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.
