r-如何计算有条件的平均值

Question

下面是生成可复制数据帧的脚本：

id <- c(1:20)
a <- as.numeric(round(runif(20,-40,40),2))
b <- as.numeric(round(a*1.4+60,2))
df <- as.data.frame(cbind(id, a, b))

我想计算"a"在不同条件下"b"的平均值。例如，当-40=<a<0；当0＝<a<40.

非常感谢！

Answer 1

这里有一个快速的data.table解决方案（假设coef是a）

library(data.table)
setDT(df)[, .(MeanASmall = mean(b[-40 <= a & a < 0]),
              MeanABig = mean(b[0 <= a & a <= 40]))]
#    MeanASmall MeanABig
# 1:   33.96727    89.46

如果a范围有限，你可以用基本R也快速完成

sapply(split(df, df$a >= 0), function(x) mean(x$b))
#     FALSE     TRUE 
#  33.96727 89.46000 

Answer 2

以下解决方案可以做到这一点：

子集

ndf1<-subset(df, a>=-40 & a<=0)
ndf2<-subset(df, a>=0 & a<=40)
mean(ndf1[,3])
mean(ndf2[,3])

或者更简单的

mean(df[a>=-40 & a<=0, 3]) 
mean(df[a>=0 & a<=40, 3]) 

使用ddply

library(plyr)
ddply(df, .(a>=-40 & a<=0), summarize, mean=mean(b))
ddply(df, .(a>=0 & a<=40), summarize, mean=mean(b))