Question

对于矩阵（as.matrix），如何生成一个表格，其中行等于矩阵的行？

>table(matrix)

和

>hist(matrix)

显示矩阵中每个唯一数据值的累积和，但我想要一个表，其中行与每个矩阵行的值相同，表列是矩阵中每个唯一数据值的总和。

示例矩阵：

   1  2  3  4 
a  5  5  4  6    
b  5  5  5  5     
c  8  7  6  6   
d  2  6  6  6     
e  7  7  5  4

所需的输出表：

   2  4  5  6  7  8
a  0  1  2  1  0  0
b  0  0  4  0  0  0
c  0  0  0  2  1  1
d  1  0  0  3  0  0
e  0  1  1  0  2  0

Answer 1

另一种方法是将您的matrix转换为长data.frame（使用stack），此时您可以轻松使用table：

这是您的数据：

mymat <- structure(c(5L, 5L, 8L, 2L, 7L, 5L, 5L, 7L, 6L, 7L, 4L, 5L, 6L, 
            6L, 5L, 6L, 5L, 6L, 6L, 4L), .Dim = c(5L, 4L), .Dimnames = list(
              c("a", "b", "c", "d", "e"), c("1", "2", "3", "4")))

这就像长data.frame：

一样

head(stack(data.frame(t(mymat))))
#   values ind
# 1      5   a
# 2      5   a
# 3      4   a
# 4      6   a
# 5      5   b
# 6      5   b

以下是我们如何使用它来创建所需的表格：

with(stack(data.frame(t(mymat))), table(ind, values))
#    values
# ind 2 4 5 6 7 8
#   a 0 1 2 1 0 0
#   b 0 0 4 0 0 0
#   c 0 0 0 2 1 1
#   d 1 0 0 3 0 0
#   e 0 1 1 0 2 0

Answer 2

## source data
x=as.matrix(read.table(text="
   1  2  3  4 
a  5  5  4  6    
b  5  5  5  5     
c  8  7  6  6   
d  2  6  6  6     
e  7  7  5  4
"))

# result

table(rep(rownames(x),ncol(x)),c(x))

#   2 4 5 6 7 8
# a 0 1 2 1 0 0
# b 0 0 4 0 0 0
# c 0 0 0 2 1 1
# d 1 0 0 3 0 0
# e 0 1 1 0 2 0

Answer 3

我也使用了apply：

t(apply(mat, 1, function(x) table(factor(x, levels = unique(sort(c(mat)))))))

R > mat  = matrix(sample(1:8, 20, replace = T), 5, 4)
R > mat
     [,1] [,2] [,3] [,4]
[1,]    5    6    1    4
[2,]    4    3    4    8
[3,]    4    8    4    3
[4,]    3    3    5    1
[5,]    1    1    3    1
R > t(apply(mat, 1, function(x) table(factor(x, levels = unique(sort(c(mat)))))))
     1 3 4 5 6 8
[1,] 1 0 1 1 1 0
[2,] 0 1 2 0 0 1
[3,] 0 1 2 0 0 1
[4,] 1 2 0 1 0 0
[5,] 3 1 0 0 0 0

Answer 4

您可以在行上使用apply，然后使用mapply和ifelse语句来取回矩阵。

假设X是你的矩阵：

# this will get you the values, just not in a nice matrix
tables.list <- apply(X, 1, table)

# unique values
vals <- sort(unique(c(X)))

# this will get you the matrix
results <- t(mapply(function(v, t)
   ifelse(v %in% names(t), t[as.character(v)], 0), list(vals), tables.list ))

# give it names
dimnames(results) <- list(rownames(X), vals)

results

#   2 4 5 6 7 8
# a 0 1 2 1 0 0
# b 0 0 4 0 0 0
# c 0 0 0 2 1 1
# d 1 0 0 3 0 0
# e 0 1 1 0 2 0

R表由矩阵行组成

4 个答案: