在R中，通过列名和行名合并矩阵

Question

我想在R中按照列名和行名合并矩阵：

m1 = matrix(c(1,2,3, 11,12,13), nrow = 2, ncol = 3, byrow = TRUE,
           dimnames = list(c("r1", "r2"),
                           c("a", "b", "c")))
m2 = matrix(c(4, 5, 0, 2,3,4), nrow = 2, ncol = 3, byrow = TRUE,
        dimnames = list(c("r2", "r3"),
                        c("d", "b", "c")))

检查m1:

> m1
    a  b  c
r1  1  2  3
r2 11 12 13

检查 m2：

> m2
   d b c
r2 4 5 0
r3 2 3 4

我希望得到以下样式的m3：

> m3
    a  b  c d
r1  1  2  3 0
r2 11 17 13 4
r3  0  3  4 2

我没有找到一种优雅的方法来实现此目的。使用plyr包中的rbind.fill.matrix函数，我可以间接地获得m3。

require(plyr)
m3 = rbind.fill.matrix(m1, m2)
rownames(m3) = c(rownames(m1), rownames(m2))
m3[is.na(m3)]=0 # replace na with zero
m3 = t(sapply(by(m3,rownames(m3),colSums),identity)) # aggregate matrix by rownames

我猜想应该有更好的方法做到这一点。你有什么建议吗？

- Frank Wang

---

你从哪里得到m3中的17？ - Rich Scriven

@RichardScriven 17来自于第r2行和b列的位置，具体而言，17 = m1 [2, 2] + m2 [1, 2]。 - Frank Wang

Answer 1

以下内容看起来是有效的：

tmp = rbind(as.data.frame(as.table(m1)), as.data.frame(as.table(m2)))
#tmp = aggregate(Freq ~ Var1 + Var2, tmp, sum)  #unnecessary   
xtabs(Freq ~ Var1 + Var2, tmp)
#    Var2
#Var1  a  b  c  d
#  r1  1  2  3  0
#  r2 11 17 13  4
#  r3  0  3  4  2

编辑：如@AnandaMahto所指出的那样，xtabs是一个“列联表”而不是“重塑数据”的函数，因此默认情况下会进行求和。

- alexis_laz

---

在这种情况下，您可以跳过聚合步骤，因为xtabs将负责对值进行求和。 - A5C1D2H2I1M1N2O1R2T1

Answer 2

我使用了这段代码：

m1 = m1[sort(rownames(m1)),sort(colnames(m1))]
m2 = m2[sort(rownames(m2)),sort(colnames(m2))]

nr = unique(c(rownames(m1),rownames(m2)))
nc = unique(c(colnames(m1),colnames(m2)))

m3 = matrix(0,nr=length(nr),nc=length(nc),dimnames=list(nr,nc))
m3[rownames(m3)%in%rownames(m1),colnames(m3)%in%colnames(m1)]=m1

m3[rownames(m3)%in%rownames(m2),colnames(m3)%in%colnames(m2)]=m3[rownames(m3)%in%rownames(m2),colnames(m3)%in%colnames(m2)]+m2

- mcereda