按组求和大于当前值的值

我测量了不同地块中树木的断面积。这是一个小例子，有两个地块，每个地块有 4 棵树：

Plot    Tree    BasalArea
1         1         4
1         2         5
1         3         7
1         4         3
2         1         4
2         2         6
2         3         9
2         4         5

在每个图中，我想计算基础面积大于焦点树的树木的基础面积总和。

例如，地块 1 中的树 1 的面积为 4。该地块内有 2 棵树的面积大于树 1：树 2 和树 3 的面积分别为 5 和 7。因此，树 1 的“BA_Larger”为 5 + 7 = 12。

同一地块中的树 2 的底面积 = 5。在地块 1 中，只有一棵树的面积比树 2 的面积大：树 3 的面积为 7。因此，树 2 的“BA_Larger”为 7。

最后，数据框应该是这样的：

Plot    Tree    BasalArea   BA_Larger
1         1        4          12
1         2        5           7
1         3        7           0
1         4        3          16
2         1        4          20
2         2        6           9
2         3        9           0
2         4        5          15

数据集非常大。我尝试计算“BA_Larger”，但没有成功。非常感谢任何帮助。

The baseR 解决方案ave():

within(df, BA_Larger <- ave(BasalArea, Plot, FUN = function(x) sapply(x, function(y) sum(x[x > y]))))

With a tidyverse风格，你也可以使用map_int() or map_dbl() from purrr.

library(dplyr)
library(purrr)

df %>%
  group_by(Plot) %>%
  mutate(BA_Larger = map_int(BasalArea, ~ sum(BasalArea[BasalArea > .]))) %>%
  ungroup()

Output

# # A tibble: 8 x 4
#    Plot  Tree BasalArea BA_Larger
#   <int> <int>     <int>     <int>
# 1     1     1         4        12
# 2     1     2         5         7
# 3     1     3         7         0
# 4     1     4         3        16
# 5     2     1         4        20
# 6     2     2         6         9
# 7     2     3         9         0
# 8     2     4         5        15

Data

df <- structure(list(Plot = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), Tree = c(1L, 
2L, 3L, 4L, 1L, 2L, 3L, 4L), BasalArea = c(4L, 5L, 7L, 3L, 4L, 
6L, 9L, 5L)), class = "data.frame", row.names = c(NA, -8L))