按组查找向量中最常见的组合

Question

我有一个包含两列的表格，即id 和item：

df <- data.frame(id=c(1,1,2,2,2,2,3,3,3,4,4,4,4,4),item=c(1,2,3,1,2,3,4,1,2,3,1,2,1,2))

我想为每个id 找到 3 个项目的最频繁组合（顺序无关紧要）。所以基本上，n 选择r 其中n = number of items within id 和r = 3。每个 id 的项目数量各不相同 - 有些超过 3 个，有些则更少。

我是 R 新手，阅读了有关 combn 和 expand.grid 的信息，但我不知道如何在我的情况下使用它们（在每个 id 中工作）。

“Find most frequent combination of values in a data.frame”是我找到的最接近的问题。

编辑：基于示例的预期答案是组合“1、2、3”，它出现在 id 2 和 4 中。

Answer 1

这是一个使用dplyr的想法

df1 <- df %>% 
        group_by(id) %>% 
        arrange(item) %>% 
        summarise(new = paste(unique(combn(item, length(unique(item)), toString)), collapse = '/'))
df1
# A tibble: 4 × 2
#     id                                             new
#  <dbl>                                           <chr>
#1     1                                            1, 2
#2     2                     1, 2, 3 / 1, 3, 3 / 2, 3, 3
#3     3                                         1, 2, 4
#4     4 1, 1, 2 / 1, 1, 3 / 1, 2, 2 / 1, 2, 3 / 2, 2, 3

names(sort(table(unlist(strsplit(df1$new, '/'))), decreasing = TRUE)[1])
#[1] "1, 2, 3"

Answer 2

library(dplyr)
grouped <- df %>% group_by(id,item) %>% summarize(count = n()) %>% arrange(desc(count))

瞧。最高计数从高到低排序。

编辑：刚刚意识到我没有完全回答你的问题。我希望我给了你一个好的开始。

Answer 3

我认为这是你想要使用 base R 的（不需要包）：

a <- aggregate(item~id, df, unique)
a <- lapply(a$item, 'length<-', max(lengths(a$item)))
m <- matrix(unlist(a), ncol=3, byrow = T)
m <- t(apply(m,1,function(x) sort(x,na.last = T)))

#     [,1] [,2] [,3]
#[1,]    1    2   NA
#[2,]    1    2    3
#[3,]    1    2    4
#[4,]    1    2    3

一旦我们得到矩阵m，矩阵中出现频率最高的行就是你想要的：

t <- table(apply(m, 1, paste, collapse = "/")) 
as.numeric(strsplit(names(which.max(t)), "/")[[1]]) 

#[1] 1 2 3