r - 获取矩阵每一列的列联表

Question

假设我有一个数字矩阵。矩阵有dim(X)=(200,5)，每个元素都在 1 到 5 之间。

我想知道每列中每个数字的计数。看起来像的东西

  X1 X2 X3 X4 X5
1 #  #  #  #  #
2 #  #  #  #  #
3 #  #  #  #  #
4 #  #  #  #  #
5 #  #  #  #  #

因为有 200 行，所以每列的总和应该是 200。

table看起来很有希望，但它只返回整个矩阵的计数，而不是列。我怎样才能做到这一点？

Answer 1

一般来说，我会这样做。例如

• 当你有一个字母矩阵时；
• 当您仍然有一个整数矩阵但它们不连续时，例如 5、7、10、11、20。

--

cX <- c(X)
k <- sort(unique(cX))
## as if we have a matrix of factors
XX <- matrix(match(cX, k), dim(X)[1], dimnames = list(k, 1:dim(X)[2]))
## aligned column-wise contingency table
tab <- apply(XX, 2, tabulate)
## aligned column-wise proportion table
prop <- tab / colSums(tab)[col(tab)]

---

我放弃了我最初的答案

lapply(data.frame(X), table)
apply(X, 2, table)

或第二个版本（一个更强大的版本，但与第一个解决方案一样低效）：

k <- sort(unique(c(X)))
apply(X, 2, function (u) table(factor(u, levels = k)) )

上面的新答案对于您的示例来说有点“矫枉过正”，但在实践中更有用（我认为）。

Answer 2

tabulate在基础 R 中怎么样：

apply(m,2,tabulate)

#    [,1] [,2] [,3] [,4] [,5]
#[1,]   39   47   38   42   34
#[2,]   41   43   41   36   39
#[3,]   46   33   38   44   39
#[4,]   35   31   40   41   53
#[5,]   39   46   43   37   35

或table：

apply(m,2,table)

数据

set.seed(1)
m <- t(replicate(200,sample(5),))