Hi,
I would like to extract monthly data automatically using a code. Given are my
sample data.
I know how to do one by one:
jan_data1 <- Pooraka_data[Pooraka_data$Month=="1",4]
feb_data1 <- Pooraka_data[Pooraka_data$Month=="2",4]
mar_data1 <- Pooraka_data[Pooraka_data$Month=="3",4]
apr_data1 <- Pooraka_data[Pooraka_data$Month=="4",4]
...
then i try:
p0_mn <- function(dt)
{ p0 <- sum(dt==0)/length(dt)
mn <- mean(dt)
c(p0=p0, mn=mn)
}
rbind(p0_mn(jan_data2),p0_mn(feb_data1),
p0_mn(mar_data1),p0_mn(apr_data1),p0_mn(may_data1),p0_mn(jun_data1),p0_mn(jul_data1),
p0_mn(aug_data1),p0_mn(sep_data1),p0_mn(oct_data1),p0_mn(nov_data1),p0_mn(dec_data1))
)
> head(Pooraka_data,20); tail(Pooraka_data,20)
Year Month Day Amount
1 1901 1 1 0.0
2 1901 1 2 3.0
3 1901 1 3 0.0
4 1901 1 4 0.5
5 1901 1 5 0.0
6 1901 1 6 0.0
7 1901 1 7 0.3
8 1901 1 8 0.0
9 1901 1 9 0.0
10 1901 1 10 0.0
11 1901 1 11 0.5
12 1901 1 12 1.8
13 1901 1 13 0.0
14 1901 1 14 0.0
15 1901 1 15 2.5
16 1901 1 16 0.0
17 1901 1 17 0.0
18 1901 1 18 0.0
19 1901 1 19 0.0
20 1901 1 20 0.0
Year Month Day Amount
32858 1990 12 17 0.0
32859 1990 12 18 0.0
32860 1990 12 19 0.8
32861 1990 12 20 0.0
32862 1990 12 21 0.0
32863 1990 12 22 8.0
32864 1990 12 23 0.0
32865 1990 12 24 0.0
32866 1990 12 25 0.0
32867 1990 12 26 0.0
32868 1990 12 27 0.4
32869 1990 12 28 0.0
32870 1990 12 29 0.0
32871 1990 12 30 0.0
32872 1990 12 31 0.0
32873 1991 1 1 0.0
32874 1991 1 2 0.0
32875 1991 1 3 0.0
32876 1991 1 4 0.0
32877 1991 1 5 5.4
thank you.
[[alternative HTML version deleted]]
