将一个矩阵的行名和列名与另一个矩阵匹配以进行逐个元素的计算

Question

我有两个矩阵，我需要进行行名和列名匹配来进行逐个元素的计算。第一个计算是phij/pij，第二个是((phij-pij)^2)/pij 长矩阵 phij 具有由破折号分隔的行名，例如啊-咩。列名没有破折号。我需要将破折号之后的行名部分（即 Baa 和 phij 矩阵中的列名与较小矩阵 pij 的行名和列名匹配。

我尝试使用 for 循环，但它与实际的行名和列名不匹配，而是查找序列中的位置。

LR<-phij
ChiSq<-phij

ROWS <- data.frame(ROW0=rownames(phij),
                     ROW1=substr(rownames(phij),regexpr("-", rownames(phij))+1,nchar(rownames(phij))))
  
  COLNAMES  <- c(colnames(phij))
  for(rowN in 1:length(ROWS$ROW0)){
    for(colN in COLNAMES){
      LR[ROWS$ROW0[rowN],colN]<-LR[ROWS$ROW0[rowN],colN]/pij[ROWS$ROW1[rowN],colN]
      ChiSq[ROWS$ROW0[rowN],colN]<-((ChiSq[ROWS$ROW0[rowN],colN]-pij[ROWS$ROW1[rowN],colN])^2)/pij[ROWS$ROW1[rowN],colN]
    }
  }

数据：

phij：

structure(c(0.111111111111111, 0.2, 0, 0, 0, 0.25, 0, 0.666666666666667, 
0.166666666666667, 0, 0, 0.666666666666667, 0.5, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.4, 0.333333333333333, 0, 1, 
0, 0.166666666666667, 0, 0.571428571428571, 0, 0, 0, 0.4, 0.272727272727273, 
0, 0, 0, 0, 0, 0, 0, 0, 0.222222222222222, 0.6, 0, 0, 0.25, 0, 
0, 0.333333333333333, 0.333333333333333, 0, 0, 0.333333333333333, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0.5, 1, 0, 0, 0, 0, 0, 0.2, 0.666666666666667, 
0, 0, 0, 0.166666666666667, 0, 0.142857142857143, 1, 0, 0, 0.2, 
0.272727272727273, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0.25, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0.222222222222222, 0.2, 0, 0, 0, 0, 
0, 0, 0.166666666666667, 0, 0, 0, 0.5, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0.2, 0, 0, 0, 0, 0.333333333333333, 0, 0.142857142857143, 
0, 0, 0, 0.2, 0.181818181818182, 0, 0, 0, 0, 0, 0, 0, 0, 0.444444444444444, 
0, 0, 0, 0.75, 0.5, 0, 0, 0.333333333333333, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0.5, 0, 0, 0, 1, 1, 0, 0.2, 0, 0, 0, 1, 0.333333333333333, 
0, 0.142857142857143, 0, 0, 1, 0, 0.272727272727273, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), class = "table", .Dim = c(49L, 
7L), .Dimnames = list(i = c("A-A", "A-Aa", "A-Aaa", "A-B", "A-Ba", 
"A-Baa", "A-Caa", "Aa-A", "Aa-Aa", "Aa-Aaa", "Aa-B", "Aa-Ba", 
"Aa-Baa", "Aa-Caa", "Aaa-A", "Aaa-Aa", "Aaa-Aaa", "Aaa-B", "Aaa-Ba", 
"Aaa-Baa", "Aaa-Caa", "B-A", "B-Aa", "B-Aaa", "B-B", "B-Ba", 
"B-Baa", "B-Caa", "Ba-A", "Ba-Aa", "Ba-Aaa", "Ba-B", "Ba-Ba", 
"Ba-Baa", "Ba-Caa", "Baa-A", "Baa-Aa", "Baa-Aaa", "Baa-B", "Baa-Ba", 
"Baa-Baa", "Baa-Caa", "Caa-A", "Caa-Aa", "Caa-Aaa", "Caa-B", 
"Caa-Ba", "Caa-Baa", "Caa-Caa"), j = c("A", "Aa", "Aaa", "B", 
"Ba", "Baa", "Caa")))

pij：

structure(c(0, 0, 0, 0, 0, 0, 0, 0, 0.608695652173913, 0.323529411764706, 
0.129032258064516, 0.176470588235294, 0.125, 0, 0, 0.173913043478261, 
0.323529411764706, 0.258064516129032, 0.294117647058824, 0.25, 
0, 0, 0.0869565217391304, 0.235294117647059, 0.419354838709677, 
0.352941176470588, 0.25, 0, 0, 0.130434782608696, 0.117647058823529, 
0.161290322580645, 0.117647058823529, 0.25, 0, 0, 0, 0, 0.032258064516129, 
0.0588235294117647, 0.125, 0, 0, 0, 0, 0, 0, 0, 0), class = "table", .Dim = c(7L, 
7L), .Dimnames = list(i = c("Aaa", "Aa", "A", "Baa", "Ba", "B", 
"Caa"), j = c("Aaa", "Aa", "A", "Baa", "Ba", "B", "Caa")))

Answer 1

您可以创建一个与phij 具有相同维度的pij 的新矩阵，然后执行您想要的计算。

new_pij <- pij[sub('-.*', '', rownames(phij)), colnames(phij)]

你可以这样做：

phij/new_pij

和

((phij-new_pij)^2)/new_pij