Question

我正在尝试以分组方式将采样函数应用于数据帧，在该数据帧中，应从每个组（如果组大小小于n 或所有组成员）中采样n个样本。

我首先尝试使用dplyr

library(dplyr)
mtcars %>% group_by(cyl) %>% sample_n(2)

这在n小于所有组的大小时有效，但是当我选择n大于组的大小时并不能容纳整个组（请注意，一个圆柱组中有7辆车）：

mtcars %>% group_by(cyl) %>% sample_n(8)
Error: `size` must be less or equal than 7 (size of data), 
set `replace` = TRUE to use sampling with replacement

我试图通过创建一个适应的group_n函数来解决此问题，例如：

sample_n_or_all <- function(tbl, n) {
  if (nrow(tbl) < n)return(tbl)
  sample_n(tbl, n)
}

但是使用我的自定义函数（mtcars %>% group_by(cyl) %>% sample_n_or_all(8)）会产生相同的错误。

有人建议我如何调整自己的功能，以便将其应用于每个组吗？还是该问题的另一种解决方案？

Answer 1

我们可以检查组中的行数，并将其相应地传递给sample_n。

library(dplyr)
n <- 8

temp <- mtcars %>% group_by(cyl) %>% sample_n(if(n() < n) n() else n) 
temp

#    mpg   cyl  disp    hp  drat    wt  qsec    vs    am  gear  carb
#   <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
# 1  21.4     4 121     109  4.11  2.78  18.6     1     1     4     2
# 2  27.3     4  79      66  4.08  1.94  18.9     1     1     4     1
# 3  24.4     4 147.     62  3.69  3.19  20       1     0     4     2
# 4  22.8     4 108      93  3.85  2.32  18.6     1     1     4     1
# 5  26       4 120.     91  4.43  2.14  16.7     0     1     5     2
# 6  33.9     4  71.1    65  4.22  1.84  19.9     1     1     4     1
# 7  30.4     4  75.7    52  4.93  1.62  18.5     1     1     4     2
# 8  30.4     4  95.1   113  3.77  1.51  16.9     1     1     5     2
# 9  21       6 160     110  3.9   2.62  16.5     0     1     4     4
#10  17.8     6 168.    123  3.92  3.44  18.9     1     0     4     4
# … with 13 more rows

之后，我们可以检查每个组中的行数。

table(temp$cyl)

#4 6 8 
#8 7 8 

table(mtcars$cyl)

# 4  6  8 
#11  7 14

Answer 2

我们可以在不使用pmin的逻辑条件的情况下做到这一点

library(dplyr)
tmp <- mtcars %>%
         group_by(cyl) %>%
         sample_n(pmin(n(), n))
# A tibble: 23 x 11
# Groups:   cyl [3]
#     mpg   cyl  disp    hp  drat    wt  qsec    vs    am  gear  carb
#   <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
# 1  33.9     4  71.1    65  4.22  1.84  19.9     1     1     4     1
# 2  27.3     4  79      66  4.08  1.94  18.9     1     1     4     1
# 3  21.4     4 121     109  4.11  2.78  18.6     1     1     4     2
# 4  30.4     4  75.7    52  4.93  1.62  18.5     1     1     4     2
# 5  21.5     4 120.     97  3.7   2.46  20.0     1     0     3     1
# 6  32.4     4  78.7    66  4.08  2.2   19.5     1     1     4     1
# 7  30.4     4  95.1   113  3.77  1.51  16.9     1     1     5     2
# 8  26       4 120.     91  4.43  2.14  16.7     0     1     5     2
# 9  17.8     6 168.    123  3.92  3.44  18.9     1     0     4     4
#10  21       6 160     110  3.9   2.62  16.5     0     1     4     4
# … with 13 more rows

-检查

table(tmp$cyl)
# 4 6 8 
# 8 7 8

自定义分组的dplyr函数（sample_n）

2 个答案: