Группа в петле в тидиверсе

Question

Могу ли я сгруппировать цикл в тидиверсе?

Большая задача - заменить переменную группировки на NA, если в группе мало наблюдений.Я хочу объединить небольшие группы в группу NA.

Однако приведенный ниже код не позволит мне group_by(x), где x - это циклическая переменная.

library(tidyverse)

for (x in c("cyl", "gear")) {
    mtcars %>%
        add_count(x) %>% 
        mutate(x = ifelse(n() < 10, NA, x))
}

Я получаю следующую ошибку.

Error in grouped_df_impl(data, unname(vars), drop) : 
  Column `x` is unknown

Answer 1

Вы также можете использовать mutate_at с table

library(tidyverse)

mtcars %>% 
  mutate_at(vars(cyl, gear), ~ {
            t <- table(.)
            ifelse(. %in% names(t[t < 10]), NA, .)})

Эту функцию можно упростить до одной строки с помощью purrr::keep

mtcars %>% 
  mutate_at(vars(cyl, gear),
            ~ ifelse(. %in% names(keep(table(.), `<`, 10)), NA, .))

Или, если вы оказалисьработая с таблицей data.table, вы можете использовать «обновление соединения» для подмножества групп с низким числом, а затем назначить NA этому подмножеству

library(data.table)
dt <- as.data.table(mtcars)

for(x in c('cyl', 'gear'))
  dt[dt[, .N, x][N < 10], on = x, (x) := NA]

Это приведет к тому же результату

all.equal(
  dt, 
  mtcars %>% 
    mutate_at(vars(cyl, gear),
    ~ ifelse(. %in% names(keep(table(.), `<`, 10)), NA, .)) %>% 
    setDT
)
# [1] TRUE

Answer 2

Вы имеете в виду что-то подобное?

library(dplyr)
for (x in c("cyl", "gear")) {
  col <- sym(x)
  mtcars <- mtcars %>%
    add_count(!!col) %>% 
    mutate(!!col := ifelse(n < 10, NA, !!col)) %>% 
    select(-n)
}
mtcars
#> # A tibble: 32 x 11
#>      mpg   cyl  disp    hp  drat    wt  qsec    vs    am  gear  carb
#>    <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
#>  1  21      NA  160    110  3.9   2.62  16.5     0     1     4     4
#>  2  21      NA  160    110  3.9   2.88  17.0     0     1     4     4
#>  3  22.8     4  108     93  3.85  2.32  18.6     1     1     4     1
#>  4  21.4    NA  258    110  3.08  3.22  19.4     1     0     3     1
#>  5  18.7     8  360    175  3.15  3.44  17.0     0     0     3     2
#>  6  18.1    NA  225    105  2.76  3.46  20.2     1     0     3     1
#>  7  14.3     8  360    245  3.21  3.57  15.8     0     0     3     4
#>  8  24.4     4  147.    62  3.69  3.19  20       1     0     4     2
#>  9  22.8     4  141.    95  3.92  3.15  22.9     1     0     4     2
#> 10  19.2    NA  168.   123  3.92  3.44  18.3     1     0     4     4
#> # ... with 22 more rows

^{Создано в 2018-12-08 пакетом Представить (v0.2.1)}

(Не самый простой синтаксис, язнаю ....)