对于示例数据框:
df1 <- structure(list(id = 1:21, region = structure(c(1L, 1L, 1L, 1L,
2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L,
4L), .Label = c("a", "b", "c", "d"), class = "factor"), weight = c(0.35,
0.65, 0.99, 1.5, 3.2, 2.1, 1.3, 3.2, 1.3, 2, 0.6, 0.6, 0.6, 0.45,
1, 1.2, 1.4, 2, 1.3, 1, 2), condition = c(0L, 1L, 0L, 1L, 0L,
0L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L
)), .Names = c("id", "region", "weight", "condition"), class = "data.frame", row.names = c(NA,
-21L))
我希望按区域排除结果变量中没有最高或最低数量1的区域。例如,我通常会这样做:
summary <- setDT(df)[,.(.result = weighted.mean((condition==1),
w = weight)*100), by = region]
哪会给我: 总结
region .result
1: a 61.60458
2: b 39.69466
3: c 50.56180
4: d 61.03896
因此,我将从数据帧df
中对区域c和d进行子集化是否可以一步完成此操作而无需手动查看摘要数据框?
答案 0 :(得分:3)
我的理解是您希望排除所有不是最高值和最低值的值。它不能作为一个衬里完成,但是如果你添加以下内容,你应该得到你想要的东西:
incl <- summary[c(which.min(.result), which.max(.result)),region]
newdf <- df1[region %in% incl,]
newdf
id region weight condition
1: 5 b 3.20 0
2: 6 b 2.10 0
3: 7 b 1.30 0
4: 8 b 3.20 1
5: 9 b 1.30 0
6: 10 b 2.00 1
7: 1 a 0.35 0
8: 2 a 0.65 1
9: 3 a 0.99 0
10: 4 a 1.50 1