获取矩阵条目的排名？

Question

假设一个矩阵：

> a <- matrix(c(100, 90, 80, 20), 2, 2)
> a
     [,1] [,2]
[1,]  100   80
[2,]   90   20

假设我想将矩阵的元素转换为行列：

>rank.a <- rank(a)
> rank.a
[1] 4 3 2 1

这将返回一个向量，即矩阵结构丢失。是否可以对矩阵进行排序，使得输出的形式为：

     [,1] [,2]
[1,]   4    2 
[2,]   3    1

Answer 1

@EDi 的答案的另一种选择是复制 a ，然后将 rank(a) 的输出直接分配到 a 副本的元素中：

> a <- matrix(c(100, 90, 80, 20), 2, 2)
> rank.a <- a
> rank.a[] <- rank(a)
> rank.a
     [,1] [,2]
[1,]    4    2
[2,]    3    1

这样您就不必通过询问输入矩阵的维度来重建矩阵。

请注意（正如 @Andrie 在评论中提到的）只有在想要保留原始 a 时才需要复制 a。需要注意的要点是，由于 a 已经具有适当的维度，因此我们可以将其视为向量，并将 a 的内容替换为 a 的秩向量代码>a。

Answer 2

为什么不将向量转换回矩阵，其尺寸与原始矩阵相同？

> a <- matrix(c(100, 90, 80, 20, 10, 5), 2, 3)
> a
     [,1] [,2] [,3]
[1,]  100   80   10
[2,]   90   20    5
> rank(a)
[1] 6 5 4 3 2 1
> rmat <- matrix(rank(a), nrow = dim(a)[1], ncol = dim(a)[2])
> rmat
     [,1] [,2] [,3]
[1,]    6    4    2
[2,]    5    3    1

Answer 3

@Gavin Simpson 有一个非常好的和优雅的解决方案！但有一个警告：

矩阵的类型将保持不变或扩大。大多数情况下您不会注意到，但请考虑以下事项：

a <- matrix( sample(letters, 4), 2, 2)
rank.a <- a
rank.a[] <- rank(a)
typeof(rank.a) # character

由于矩阵一开始就是字符，因此 rank 值（双精度）被强制转换为字符串！

这是一种更安全的方法，只需复制所有属性：

a <- matrix( sample(letters, 4), 2, 2)
rank.a <- rank(a)
attributes(rank.a) <- attributes(a)
typeof(rank.a) # double

或者，作为单行使用 struct 仅复制相关属性（但需要更多输入）：

a <- matrix( sample(letters, 4), 2, 2)
rank.a <- structure(rank(a), dim=dim(a), dimnames=dimnames(a))

当然，dimnames在这种特殊情况下可以被忽略。