我想在R中按照列名和行名合并矩阵:
我希望得到以下样式的m3:
我猜想应该有更好的方法做到这一点。你有什么建议吗?
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
我猜想应该有更好的方法做到这一点。你有什么建议吗?