比较r中的样本频率

Question

我想比较来自两个不同观察结果的样本频率。问题是第一个不包含第二个数字的整个范围。我怎么能在不编写for循环的情况下根据count返回的x值对它们进行组合？ 这是一个MWE澄清：

library(plyr)
a <- c(5, 4, 5, 7, 3, 5, 6, 5, 5, 4, 5, 5, 4, 5, 4, 7, 2, 4, 4, 5, 3, 6, 5, 6, 4, 4, 5, 4, 5, 5, 6, 7, 4)
b <- c(1, 3, 4, 6, 2, 7, 7, 4, 3, 6, 6, 3, 6, 6, 5, 6, 6, 5)

a.count <- count(a)
b.count <- count(b)

我想要的结果看起来应该是这样的：

   freq.a freq.b                   
1  1         
2  1       1                                       
3  3       2   
4  2      10                                              
5  2      13
6  7       4                                            
7  2       3

Answer 1

如果你把数据放在长格式中（每次观察一行，有一个变量来自哪个样本），那么你就可以制作一个列联表：

    data.frame(v=df.a, s='a') %>% rbind(data.frame(v=df.b, s='b')) %>%
      xtabs(f=~v+s)

产地：

   s
v    a  b
  1  0  1
  2  1  1
  3  2  3
  4 10  2
  5 13  2
  6  4  7
  7  3  2

Answer 2

df <- merge(a.count, b.count, by ='x', all=TRUE)[2:3]
names(df) <- c('freq.a', 'freq.b')
df

  freq.a freq.b
1     NA      1
2      1      1
3      2      3
4     10      2
5     13      2
6      4      7
7      3      2