R 根据确切的行列名称插入矩阵值

发布于05月07日

我有一个具有相对较大的行名和列名的空白矩阵(6000乘以200).然后，我有多个矩阵，它们表示这个矩阵的子集(通常大约是1500乘以120)，还包含行名和列名.我想根据行/列名的组合将给定子集中单元格的信息插入到空白矩阵中.

以下是数据和所需输出的示例:

adj1 <- matrix(0, 3, 4)
colnames(adj1) <- c("g1", "g2", "g3", "g4")
rownames(adj1) <- c("Tim", "John", "Sarah")

adj2 <- matrix(0, 2, 2)
colnames(adj2) <- c("g1", "g4")
rownames(adj2) <- c("Tim", "John")
adj2[1,2] <- 1
adj2

combined_adj <- matrix(0,3,4)
colnames(combined_adj) <- c("g1","g2","g3","g4")
rownames(combined_adj) <- c("Tim","John", "Sarah")
combined_adj[1,4] <- 1
combined_adj

我在某种程度上设法通过列名进行了交集，但不是通过行/列名称的确切组合.

推荐答案

Here is a function that allows you to replace values by row and column names.
Its arguments are

x矩阵越大；
y具有新值的矩阵；
z目标矩阵；如果缺少，将创建一个新的零矩阵，其维度和行/列名称为x.

matreplace <- function(x, y, z){
  if(missing(z)) {
    z <- x
    z[] <- 0
  }
  ij <- Map(intersect, dimnames(z), dimnames(y))
  i <- ij[[1L]]
  j <- ij[[2L]]
  z[i, j] <- y[i, j]
  z
}

adj1 <- matrix(0, 3, 4)
colnames(adj1) <- c("g1", "g2", "g3", "g4")
rownames(adj1) <- c("Tim", "John", "Sarah")

adj2 <- matrix(0, 2, 2)
colnames(adj2) <- c("g1", "g4")
rownames(adj2) <- c("Tim", "John")
adj2[1,2] <- 1

combined_adj <- matreplace(adj1, adj2)

adj2[2, 1] <- pi
matreplace(adj1, adj2, combined_adj)
#>             g1 g2 g3 g4
#> Tim   0.000000  0  0  1
#> John  3.141593  0  0  0
#> Sarah 0.000000  0  0  0

^{创建于2023-05-07年第reprex v2.0.2页}