1

我有一个填充了下三角形的矩阵，并想用下三角形的匹配值填充矩阵的上三角形。

我一直在尝试的代码是：

    r <-  read.table("rcorrected.csv", header = TRUE, sep = ",", check.names = FALSE)
    m <- as.matrix(r)


m[upper.tri(m, diag = FALSE)] <- m[lower.tri(m, diag= FALSE)]

输出：

Warning message:
In m[upper.tri(m, diag = FALSE)] <- m[lower.tri(m, diag = FALSE)] :
  number of items to replace is not a multiple of replacement length


    structure(m)
            1    2    3    4    5    6    7    8    9   10   11   12
 [1,]  1 1.00   NA   NA   NA   NA   NA   NA   NA   NA   NA   NA   NA
 [2,]  2 0.50 1.00   NA   NA   NA   NA   NA   NA   NA   NA   NA   NA
 [3,]  3 0.57 0.50 1.00   NA   NA   NA   NA   NA   NA   NA   NA   NA
 [4,]  4 0.00 0.04 0.16 1.00   NA   NA   NA   NA   NA   NA   NA   NA

dput(head(m[, c(2,1)]))
    structure(c(1, 0.5, 0.57, 0, 0, 0.23, 1, 2, 3, 4, 5, 6), .Dim = c(6L, 
    2L), .Dimnames = list(NULL, c("1", "")))

矩阵 m 是否被读取为列表？

4

2 回答 2

2

在这里，我们在分配diag = TRUE的 rhs 上分配不相等数量的观察值，并且diag = TRUE应该lower.tri在函数内部而不是函数外部

m[upper.tri(m, diag = FALSE)] <- m[lower.tri(m, diag= FALSE)]
于 2020-01-09T19:11:59.167 回答
1

如果m是一个矩阵，如最后注释中所示的矩阵，下三角形填充，那么这也填充了上三角形部分。 

m2 <- m + t(m)
m2
##      [,1] [,2] [,3] [,4]
## [1,]    0    5    9   13
## [2,]    5    0   10   14
## [3,]    9   10    0   15
## [4,]   13   14   15    0

如果 的对角线m不是 0，那么这样做：

m + t(m) - diag(diag(m))

笔记

m <- matrix(c(0, 5, 9, 13, 0, 0, 10, 14, 0, 0, 0, 15, 0, 0, 0, 0), 4)
m
##      [,1] [,2] [,3] [,4]
## [1,]    0    0    0    0
## [2,]    5    0    0    0
## [3,]    9   10    0    0
## [4,]   13   14   15    0

如果我们从一个数据框开始，DF那么首先使用as.matrix它来转换它。

DF <- structure(list(V1 = c(0, 5, 9, 13), V2 = c(0, 0, 10, 14), V3 = c(0, 
0, 0, 15), V4 = c(0, 0, 0, 0)), class = "data.frame", row.names = c(NA, 
-4L))
m <- as.matrix(DF)
于 2020-01-09T20:57:28.777 回答