这样的事情应该可以完成(在基础 R 中):
transform(mtcars, dif_mpg=mpg-ave(mpg, cyl, FUN=mean))
ave 在cyl 定义的mpg 的子组上计算FUN。 transform 允许您向数据框添加/修改列,还可以在数据框的上下文中评估表达式(因此您不必输入mtcars$mpg 等)。以下是结果的前 6 行:
mpg cyl disp hp drat wt qsec vs am gear carb dif_mpg
Mazda RX4 21.0 6 160.0 110 3.90 2.620 16.46 0 1 4 4 1.25714286
Mazda RX4 Wag 21.0 6 160.0 110 3.90 2.875 17.02 0 1 4 4 1.25714286
Datsun 710 22.8 4 108.0 93 3.85 2.320 18.61 1 1 4 1 -3.86363636
Hornet 4 Drive 21.4 6 258.0 110 3.08 3.215 19.44 1 0 3 1 1.65714286
Hornet Sportabout 18.7 8 360.0 175 3.15 3.440 17.02 0 0 3 2 3.60000000
Valiant 18.1 6 225.0 105 2.76 3.460 20.22 1 0 3 1 -1.64285714
其他替代方案包括 dplyr 包(如 David Robinson 所示)、data.table:
library(data.table)
(data.table(mtcars, keep.rownames=T)[, dif_mpg:=mpg - mean(mpg), by=cyl])
还有plyr(虽然你应该使用dplyr而不是plyr,因为它更快):
library(plyr)
ddply(mtcars, "cyl", transform, dif_mpg=mpg-mean(mpg))