按行顺序替换上三角矩阵元素

Question

我有一个方阵如下

> ex_mat
          [,1]      [,2]      [,3]      [,4]
[1,] 0.4270634 2.1920890 0.5647472 1.7149861
[2,] 2.0556220 1.1157322 2.6723637 0.3155507
[3,] 1.2252602 0.1063053 0.6396099 0.7903348
[4,] 0.3614062 1.1118661 0.5000143 0.2491543

我用这个按行顺序（从最大到最小）对矩阵的上部非对角线部分进行了排序。

> rank(-(t(ex_mat)[lower.tri(ex_mat)]))
[1] 2 5 3 1 6 4

我想用上面获得的秩替换“ex_mat”的上部非对角线元素。我使用 https://statisticsglobe.com/modify-diagonal-lower-upper-triangular-part-matrix-r 作为示例代码。虽然我的排名是正确的，但代码似乎是按列顺序插入排名。

> ex_mat_new <- ex_mat
> ex_mat_new[upper.tri(ex_mat_new)] <- rank(-(t(ex_mat)[lower.tri(ex_mat)]))
> ex_mat_new
          [,1]      [,2]      [,3]      [,4]
[1,] 0.4270634 2.0000000 5.0000000 1.0000000
[2,] 2.0556220 1.1157322 3.0000000 6.0000000
[3,] 1.2252602 0.1063053 0.6396099 4.0000000
[4,] 0.3614062 1.1118661 0.5000143 0.2491543

我该如何解决这个问题？ [1,4] 和 [2,3] 关闭。谢谢。

Answer 1

值插入按列顺序（在矩阵中，data.frame）。我们可以在 lower.tri 上赋值，然后得到 transpose

ex_mat_new[lower.tri(ex_mat_new)] <- rank(-(t(ex_mat)[lower.tri(ex_mat)]))
ex_mat_new <- t(ex_mat_new)
ex_mat_new[lower.tri(ex_mat_new)] <- ex_mat[lower.tri(ex_mat)]

-输出

> ex_mat_new
          [,1]      [,2]      [,3]      [,4]
[1,] 0.4270634 2.0000000 5.0000000 3.0000000
[2,] 2.0556220 1.1157322 1.0000000 6.0000000
[3,] 1.2252602 0.1063053 0.6396099 4.0000000
[4,] 0.3614062 1.1118661 0.5000143 0.2491543

或者这可以用 replace

在一行中完成

t(replace(t(ex_mat), lower.tri(ex_mat), rank(-(t(ex_mat)[lower.tri(ex_mat)]))))

-输出

       [,1]      [,2]      [,3]      [,4]
[1,] 0.4270634 2.0000000 5.0000000 3.0000000
[2,] 2.0556220 1.1157322 1.0000000 6.0000000
[3,] 1.2252602 0.1063053 0.6396099 4.0000000
[4,] 0.3614062 1.1118661 0.5000143 0.2491543

数据

ex_mat <- structure(c(0.4270634, 2.055622, 1.2252602, 0.3614062, 2.192089, 
1.1157322, 0.1063053, 1.1118661, 0.5647472, 2.6723637, 0.6396099, 
0.5000143, 1.7149861, 0.3155507, 0.7903348, 0.2491543), .Dim = c(4L, 
4L), .Dimnames = list(NULL, NULL))

按行顺序替换上三角矩阵元素

Replacing upper triangular matrix elements in row order

r

matrix

数据