用 NA 填充两个矩阵的缺失数据

您可以使用它来查找匹配的列和行，否则它将添加带有 NA 的列或行。用于对每个矩阵执行操作，以将结果生成为列表。您可以使用以下代码：matchlapply

l <- list(t1, t2)

lapply(l, \(x) {
  # Your required columns and rows names
  colrown=c("a","b","c","d")
  # columns part
  coln=colnames(x)
  x=x[,match(colrown,coln)]
  colnames(x)=colrown
  # Rows part
  rown=rownames(x)
  x=x[match(colrown,rown),]
  rownames(x)=colrown
  return(x)
})
#> [[1]]
#>    a  b  c  d
#> a  1  0  1 NA
#> b  0  0  1 NA
#> c  1  0  1 NA
#> d NA NA NA NA
#> 
#> [[2]]
#>    a  b  c  d
#> a  1 NA  1  0
#> b NA NA NA NA
#> c  0 NA  0  1
#> d  0 NA  0  1

^{创建于 2022-09-24 使用 reprex v2.0.2}

数据

#first matrix
t1 = matrix(
  c(1, 0, 1, 0, 0, 1, 1, 0, 1),
  nrow = 3,  
  ncol = 3,        
  byrow = TRUE         
)
rownames(t1) <- c("a","b", "c")
colnames(t1) <- c("a","b", "c")

#second matrix
t2 = matrix(
  c(1, 1, 0, 0, 0, 1, 0, 0, 1),
  nrow = 3,  
  ncol = 3,        
  byrow = TRUE         
)
rownames(t2) <- c("a","c", "d") 
colnames(t2) <- c("a","c", "d")

t1
#>   a b c
#> a 1 0 1
#> b 0 0 1
#> c 1 0 1
t2
#>   a c d
#> a 1 1 0
#> c 0 0 1
#> d 0 0 1

如果每个矩阵中的行名和列名始终相同（此处没有硬编码的行/列名）：

nam <- union(rownames(t1), rownames(t2))
m <- array(dim=rep(length(nam), 2),
           dimnames=list(nam, nam))
lapply(list(t1, t2), function(x) {
  m[rownames(x), colnames(x)] <- x
  m
})

# [[1]]
#    a  b  c  d
# a  1  0  1 NA
# b  0  0  1 NA
# c  1  0  1 NA
# d NA NA NA NA
# 
# [[2]]
#    a  b  c  d
# a  1 NA  1  0
# b NA NA NA NA
# c  0 NA  0  1
# d  0 NA  0  1