Why it is so slow? A little research located a mail group posting from a Aug. 2011 where @hadley, the package author, states
This is a drawback of the way that ddply always works with data
frames. It will be a bit faster if you use summarise instead of
data.frame (because data.frame is very slow), but I’m still thinking
about how to overcome this fundamental limitation of the ddply
approach.
As for being efficient plyr code I didn’t know either. After a bunch of param testing and bench-marking it looks like we can do better.
The summarize()
in your command is a just helper function, pure and simple. We can replace it with our own sum function since it isn’t helping with anything that isn’t already simple and the .data
and .(price)
arguments can be made more explicit. The result is
ddply( dd[, 2:3], ~price, function(x) sum( x$volume ) )
The summarize
may seem nice, but it just isn’t quicker than a simple function call. It makes sense; just look at our little function versus the code for summarize
. Running your benchmarks with the revised formula yields a noticeable gain. Don’t take that to mean you’ve used plyr incorrectly, you haven’t, it just isn’t efficient; nothing you can do with it will make it as fast as other options.
In my opinion the optimized function still stinks as it isn’t clear and must be mentally parsed along with still being ridiculously slow compared with data.table ( even with a 60% gain ).
In the same thread mentioned above, regarding the slowness of plyr, a plyr2 project is mentioned. Since the time of the original answer to the question the plyr author has released dplyr
as the successor of plyr. While both plyr and dplyr are billed as data manipulation tools and your primary stated interest is aggregation you may still be interested in your benchmark results of the new package for comparison as it has a reworked backend to improve performance.
plyr_Original <- function(dd) ddply( dd, .(price), summarise, ss=sum(volume))
plyr_Optimized <- function(dd) ddply( dd[, 2:3], ~price, function(x) sum( x$volume ) )
dplyr <- function(dd) dd %.% group_by(price) %.% summarize( sum(volume) )
data_table <- function(dd) dd[, sum(volume), keyby=price]
The dataframe
package has been removed from CRAN and subsequently from the tests, along with the matrix function versions.
Here’s the i=5, j=8
benchmark results:
$`obs= 500,000 unique prices= 158,286 reps= 5`
test elapsed relative
9 data_table(d.dt) 0.074 1.000
4 dplyr(d.dt) 0.133 1.797
3 dplyr(d.df) 1.832 24.757
6 l.apply(d.df) 5.049 68.230
5 t.apply(d.df) 8.078 109.162
8 agg(d.df) 11.822 159.757
7 b.y(d.df) 48.569 656.338
2 plyr_Optimized(d.df) 148.030 2000.405
1 plyr_Original(d.df) 401.890 5430.946
No doubt the optimizing helped a bit. Take a look at the d.df
functions; they just can’t compete.
For a little perspective on the slowness of the data.frame structure here are micro-benchmarks of the aggregation times of data_table and dplyr using a larger test dataset (i=8,j=8
).
$`obs= 50,000,000 unique prices= 15,836,476 reps= 5`
Unit: seconds
expr min lq median uq max neval
data_table(d.dt) 1.190 1.193 1.198 1.460 1.574 10
dplyr(d.dt) 2.346 2.434 2.542 2.942 9.856 10
dplyr(d.df) 66.238 66.688 67.436 69.226 86.641 10
The data.frame is still left in the dust. Not only that, but here’s the elapsed system.time to populate the data structures with the test data:
`d.df` (data.frame) 3.181 seconds.
`d.dt` (data.table) 0.418 seconds.
Both creation and aggregation of the data.frame is slower than that of the data.table.
Working with the data.frame in R is slower than some alternatives but as the benchmarks show the built in R functions blow plyr out of the water. Even managing the data.frame as dplyr does, which improves upon the built-ins, doesn’t give optimal speed; where as data.table is faster both in creation and aggregation and data.table does what it does while working with/upon data.frames.
In the end…
Plyr is slow because of the way it works with and manages the data.frame manipulation.
[punt:: see the comments to the original question].
## R version 3.0.2 (2013-09-25)
## Platform: x86_64-pc-linux-gnu (64-bit)
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] microbenchmark_1.3-0 rbenchmark_1.0.0 xts_0.9-7
## [4] zoo_1.7-11 data.table_1.9.2 dplyr_0.1.2
## [7] plyr_1.8.1 knitr_1.5.22
##
## loaded via a namespace (and not attached):
## [1] assertthat_0.1 evaluate_0.5.2 formatR_0.10.4 grid_3.0.2
## [5] lattice_0.20-27 Rcpp_0.11.0 reshape2_1.2.2 stringr_0.6.2
## [9] tools_3.0.2