我想在数据框中计算组平均值，并在包含这些组平均值的原始数据框中创建一个新列。 （我正在进行重复性研究，我希望在新列中插入，单位和通道内的测量值的平均值，以便我可以减去它并计算残差。）
我的数据：
> head(mytestdata,15)
Insertion Measurement Unit Channel Value
1 1 1 A5 10 9.41
2 1 1 A5 11 9.51
3 1 1 A5 12 10.59
4 1 1 A5 13 9.45
5 1 2 A5 10 9.42
6 1 2 A5 11 9.03
7 1 2 A5 12 10.62
8 1 2 A5 13 9.39
9 1 3 A5 10 9.38
10 1 3 A5 11 9.87
11 1 3 A5 12 11.34
12 1 3 A5 13 9.59
13 2 1 A5 10 12.10
14 2 1 A5 11 11.28
15 2 1 A5 12 12.95
具体来说，我想计算每次插入的平均值，单位和通道，并将其作为meanValue添加到数据框。然后从Value中减去meanValue以获得Residual。
应该是这样的：
Insertion Measurement Unit Channel Value meanValue
1 1 1 40 10 11.79 11.56
2 1 1 40 11 11.01 11.38
3 1 1 40 12 10.86 11.19
4 1 1 40 13 10.29 10.91
5 1 2 40 10 11.47 11.56
6 1 2 40 11 11.84 11.38
7 1 2 40 12 11.39 11.19
8 1 2 40 13 11.25 10.91
9 1 3 40 10 11.42 11.56
10 1 3 40 11 11.28 11.38
11 1 3 40 12 11.31 11.19
12 1 3 40 13 11.18 10.91
13 2 1 40 10 10.97 11.55
14 2 1 40 11 11.78 11.87
15 2 1 40 12 11.48 11.25
我知道如何让这个组意味着使用by，aggregate等，这会得到第二个列表或表格，其中包含值。我也相信我可以使用一些复杂的循环程序得到我想要的东西，但是我想把它们放回到优雅的单行或双行解决方案中的相同数据框中，我认为必须有一个这样做的方法，但经过几天的搜索，我找不到它。我不想要一个繁琐的解决方案，因为当我扩展到更多数据时，我希望它能很好地工作。
答案 0 :(得分：3)
您可以使用ave来计算groupmeans：
df$MeanValue <- with(df, ave(Value, Insertion, Unit, Channel, FUN = mean))
然后计算残差：
df$Residual <- df$Value - df$MeanValue
df
# Insertion Measurement Unit Channel Value MeanValue Residual
#1 1 1 A5 10 9.41 9.403333 0.006666667
#2 1 1 A5 11 9.51 9.470000 0.040000000
#3 1 1 A5 12 10.59 10.850000 -0.260000000
#4 1 1 A5 13 9.45 9.476667 -0.026666667
#5 1 2 A5 10 9.42 9.403333 0.016666667
#6 1 2 A5 11 9.03 9.470000 -0.440000000
#7 1 2 A5 12 10.62 10.850000 -0.230000000
#8 1 2 A5 13 9.39 9.476667 -0.086666667
#9 1 3 A5 10 9.38 9.403333 -0.023333333
#10 1 3 A5 11 9.87 9.470000 0.400000000
#11 1 3 A5 12 11.34 10.850000 0.490000000
#12 1 3 A5 13 9.59 9.476667 0.113333333
#13 2 1 A5 10 12.10 12.100000 0.000000000
#14 2 1 A5 11 11.28 11.280000 0.000000000
#15 2 1 A5 12 12.95 12.950000 0.000000000
或者您可以使用dplyr
library(dplyr)
df %>% group_by(Insertion, Unit, Channel) %>% mutate(MeanValue = mean(Value), Residual = Value - MeanValue)
答案 1 :(得分：2)
使用data.table
library(data.table)
setDT(mytestdata)[, c("MeanValue", "Residual") := {m= mean(Value);list(m, Value-m)}, by=list(Insertion, Unit, Channel)]
mytestdata
# Insertion Measurement Unit Channel Value MeanValue Residual
# 1: 1 1 A5 10 9.41 9.403333 0.006666667
# 2: 1 1 A5 11 9.51 9.470000 0.040000000
# 3: 1 1 A5 12 10.59 10.850000 -0.260000000
# 4: 1 1 A5 13 9.45 9.476667 -0.026666667
# 5: 1 2 A5 10 9.42 9.403333 0.016666667
# 6: 1 2 A5 11 9.03 9.470000 -0.440000000
# 7: 1 2 A5 12 10.62 10.850000 -0.230000000
# 8: 1 2 A5 13 9.39 9.476667 -0.086666667
# 9: 1 3 A5 10 9.38 9.403333 -0.023333333
# 10: 1 3 A5 11 9.87 9.470000 0.400000000
# 11: 1 3 A5 12 11.34 10.850000 0.490000000
# 12: 1 3 A5 13 9.59 9.476667 0.113333333
# 13: 2 1 A5 10 12.10 12.100000 0.000000000
# 14: 2 1 A5 11 11.28 11.280000 0.000000000
# 15: 2 1 A5 12 12.95 12.950000 0.000000000