矩阵 - 镜像下三角形到上三角形

Posted 2021-03-29

如何用diag翻转lower.tri，没有最后一行翻到upper.tri？

m <- structure(c(1, 2, 3, 4, 5, NA, 66, 77, 88, 99, NA, NA, 333, 444, 
555, NA, NA, NA, 101, 202, NA, NA, NA, NA, 99999), .Dim = c(5L, 
5L), .Dimnames = list(c("v_2", "v_3", "v_3a", "v_3b", "v_4"), 
    c("v_0", "v_2", "v_3", "v_3a", "v_3b")))

好像：

     v_0 v_2 v_3 v_3a  v_3b
v_2    1  NA  NA   NA    NA
v_3    2  66  NA   NA    NA
v_3a   3  77 333   NA    NA
v_3b   4  88 444  101    NA
v_4    5  99 555  202 99999

期望的结果：

     v_0 v_2 v_3 v_3a  v_3b
v_2    1   1   2    3     4
v_3    2  66  66   77    88
v_3a   3  77 333  333   444
v_3b   4  88 444  101   101
v_4    5  99 555  202 99999

相当肯定的功能lower.tri()和upper.tri()可能会有所帮助。

答案

你可以做：

tm <- t(m)[,-nrow(m)] # or t(m[-nrow(m),])
m[upper.tri(m)] <- tm[upper.tri(tm, diag = TRUE)]
# > m
#      v_0 v_2 v_3 v_3a  v_3b
# v_2    1   1   2    3     4
# v_3    2  66  66   77    88
# v_3a   3  77 333  333   444
# v_3b   4  88 444  101   101
# v_4    5  99 555  202 99999