使用行名和列名比较 R 中的矩阵

Question

我试图根据行和列名比较两个5x5矩阵。请注意，在矩阵B中，第2行和第4行的名称已交换。

a <- matrix(1:25, nrow = 5)
colnames(a) = c("col1", "col2", "col3", "col4", "col5")
rownames(a) = c("row1", "row2", "row3", "row4", "row5")
a

b <- matrix(1:25, nrow = 5)
colnames(b) = c("col1", "col2", "col3", "col4", "col5")
rownames(b) = c("row1", "row4", "row3", "row2", "row5")
b

c <- b-a
c

     col1 col2 col3 col4 col5
row1    0    0    0    0    0
row4    0    0    0    0    0
row3    0    0    0    0    0
row2    0    0    0    0    0
row5    0    0    0    0    0

计算差异C = B- a时，代码返回零的矩阵，就好像它在寻找矩阵中元素的位置，而不是行和列索引。我希望第4行的结果是：

d <- a["row4",] - b["row4",]
col1 col2 col3 col4 col5 
   2    2    2    2    2

有什么想法如何做到这一点？

Answer 1

以下示例将使用公共行和列从a中减去b。

ro <- intersect(rownames(a), rownames(b))
co <- intersect(colnames(a), colnames(b))

a[ro,co] - b[ro,co]

##>      col1 col2 col3 col4 col5
##> row1    0    0    0    0    0
##> row2   -2   -2   -2   -2   -2
##> row3    0    0    0    0    0
##> row4    2    2    2    2    2
##> row5    0    0    0    0    0

Answer 2

a - b[sort(rownames(b)),]

     col1 col2 col3 col4 col5
row1    0    0    0    0    0
row2   -2   -2   -2   -2   -2
row3    0    0    0    0    0
row4    2    2    2    2    2
row5    0    0    0    0    0

如果列名的顺序也不同，则可以在逗号后添加sort（colnames（b））。

顺便说一句，我建议不要使用c作为r中对象的名称，因为c是创建向量的内置功能。