[R] Oh apply functions, how you confuse me
ROLL Josh F
JRoll at lcog.org
Thu Sep 1 21:55:52 CEST 2011
Dang Jim this looks to do the trick though I never heard of a data.table, interesting, I will explore more. Thanks you very much.
-----Original Message-----
From: jim holtman [mailto:jholtman at gmail.com]
Sent: Thursday, September 01, 2011 11:20 AM
To: ROLL Josh F
Cc: r-help at r-project.org
Subject: Re: [R] Oh apply functions, how you confuse me
Is this close to what you are asking for:
> require(data.table)
> Dt.. <- data.table(Df..)
> R <- Dt..[,
+ list(
+ sum = sum(Volume)
+ , weight = sum(Volume * Mph) / sum(Volume)
+ )
+ , by = list(Min5Break, Day, Hour, Dir)
+ ]
> R
Min5Break Day Hour Dir sum weight
1 1 0 NB 730.8880 32.60224
2 1 0 NB 766.4083 35.88443
3 1 0 SB 776.7592 32.66822
4 1 0 SB 768.0923 33.55988
5 1 0 NB 767.5472 36.00546
6 1 0 NB 767.6600 30.38747
7 1 0 SB 814.9662 31.88483
8 1 0 SB 795.4855 30.91495
9 1 0 NB 828.4439 31.57477
10 1 0 NB 797.7522 29.49832
11 1 0 SB 826.5165 32.74487
12 1 0 SB 824.0942 36.28309
1 2 1 NB 830.0683 29.59320
2 2 1 NB 838.8179 34.59878
3 2 1 SB 877.3518 30.77636
4 2 1 SB 838.9765 30.90577
5 2 1 NB 736.6560 30.54381
6 2 1 NB 772.3622 31.40094
7 2 1 SB 819.2347 29.22674
8 2 1 SB 840.9048 32.59222
9 2 1 NB 818.8383 37.55142
10 2 1 NB 783.8896 32.54565
11 2 1 SB 699.0401 30.76466
12 2 1 SB 773.5594 35.87076
cn Min5Break Day Hour Dir sum weight
>
On Thu, Sep 1, 2011 at 1:11 PM, LCOG1 <jroll at lcog.org> wrote:
> Hi guys,
> I have a crap load of data to parse and have enjoyed creating a script
> that takes this data and creates a number of useful graphics for our
> area. I am unable to figure out one summary though and its all cause
> I dont fully understand the apply family of functions. Consider the following:
>
>
>
> #Create data
> Df..<-rbind(data.frame(Id=1:1008,Dir=rep(c("NB","NB","SB","SB"),252),M
> ph=runif(1008,0,65),
> Volume=runif(1008,0,19),Hour=rep(00,1008),Min5Break=rep(1:12,84),Day=r
> ep(1,1008)),
> data.frame(Id=2009:2016,Dir=rep(c("NB","NB","SB","SB"),252),Mph=runif(
> 1008,0,65),
> Volume=runif(1008,0,19),Hour=rep(01,1008),Min5Break=rep(1:12,84),Day=r
> ep(2,1008)))
>
> #Example calc
> Results_<-list()
>
> #Sum Volume by 5 minute break by Day by Direction
> Results_$FiveMin.Direction<-tapply(Df..$Volume,list(Df..$Min5Break,Df.
> .$Day,Df..$Hour,Df..$Dir),sum)
>
> The data is a snap shot of what im working with and I am trying to get
> to something similar to the last line where the volumes are summed.
> What i want to do is to do a weighted average for the speed by 5
> minute break. So for all the speeds and volumes in a given hour of 5
> minute break(12 per hour), i would want to
>
> sum(Volumes[1:12]*Speed[1:12]) / sum(Volumes[1:12]
>
> The output resembling the one from the above but having these weighted
> values. I am assuming the sum function in the above would be replaced
> by a function doing the calculation but I am still not sure how to do
> this using apply functions, so perhaps this isnt the best option.
>
> Hope this is clear and hope you guys(and of course ladies) can offer
> some guidance.
>
> Cheers,
> Josh
>
>
>
>
> --
> View this message in context:
> http://r.789695.n4.nabble.com/Oh-apply-functions-how-you-confuse-me-tp
> 3784212p3784212.html Sent from the R help mailing list archive at
> Nabble.com.
>
> ______________________________________________
> R-help at 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.
>
--
Jim Holtman
Data Munger Guru
What is the problem that you are trying to solve?
More information about the R-help
mailing list