我想通过 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
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
我想一定有一些更好的方法可以做到这一点。你有什么建议?