使用data.table按组进行子集化

had*_*ley 51 r greatest-n-per-group data.table

假设我有一个包含一些棒球运动员的数据表:

library(plyr)
library(data.table)

bdt <- as.data.table(baseball)
Run Code Online (Sandbox Code Playgroud)

对于每个玩家(由id给出),我想找到与他们玩最多游戏的年份相对应的行.这在plyr中很简单:

ddply(baseball, "id", subset, g == max(g))
Run Code Online (Sandbox Code Playgroud)

data.table的等效代码是什么?

我试过了:

setkey(bdt, "id") 
bdt[g == max(g)]  # only one row
bdt[g == max(g), by = id]  # Error: 'by' or 'keyby' is supplied but not j
bdt[, .SD[g == max(g)]] # only one row
Run Code Online (Sandbox Code Playgroud)

这有效:

bdt[, .SD[g == max(g)], by = id] 
Run Code Online (Sandbox Code Playgroud)

但它比plyr快30%,这表明它可能不是惯用语.

edd*_*ddi 69

这是快速的data.table方式:

bdt[bdt[, .I[g == max(g)], by = id]$V1]
Run Code Online (Sandbox Code Playgroud)

这避免了构造.SD,这是表达式中的瓶颈.

编辑:实际上,OP速度慢的主要原因不仅仅在于它.SD,而是它以特定的方式使用它 - 通过调用[.data.table,此时有很大的开销,所以在循环中运行它(当一个人做a by)累积一个非常大的惩罚.

  • 最终将优化原始方法.参见[FR 2330](https://r-forge.r-project.org/tracker/index.php?func=detail&aid=2330&group_id=240&atid=978)优化`.SD [i]`查询以保持优雅但让它更快不变. (6认同)
  • +1我打赌Hadley想要以编程方式做这件事,在这种情况下他想要使用这种语法,`bdt [bdt [,.I [g == max(g)],by = id] [ ,V1]]`对吗? (5认同)
  • @joran我正在手动构建调用,所以这并不重要 (2认同)
  • 该问题链接自此处从R-Forge移至GitHub [#613](https://github.com/Rdatatable/data.table/issues/613) (2认同)