正如标题所述,我需要根据另一个矩阵 (B) 中给出的列和行名称,将对称矩阵 (A) 中的值替换为 0。例如,
A <-matrix(c(0,1,2,3,4,1,0,1,2,3,2,1,0,1,2,3,2,1,0,3,4,3,2,3,0), 5)
colnames(A) <- rownames(A) <-c(LETTERS[1:5])
A
A B C D E
A 0 1 2 3 4
B 1 0 1 2 3
C 2 1 0 1 2
D 3 2 1 0 3
E 4 3 2 3 0
B <- matrix(c("A","B","C","B","D","E","E","C"),4)
B #reference matrix
[,1] [,2]
[1,] "A" "D"
[2,] "B" "E"
[3,] "C" "E"
[4,] "B" "C"
Run Code Online (Sandbox Code Playgroud)
我试过
A[B[,1],B[,2]]<-0 #for the above-diagonal half but it didn't work as I expected:
A
A B C D E
A 0 1 0 0 0
B 1 0 0 0 0
C 2 1 0 0 0
D 3 2 1 0 3
E 4 3 2 3 0
Run Code Online (Sandbox Code Playgroud)
具体来说,我需要将A[A,D], A[B,E], A[C,E]and替换A[B,C]为 0(对角线上方),并用A[D,A]etc.(对角线下方)替换。
像这样的东西?:
A[rbind(B, B[,2:1])] <- 0 # B[,2:1] swaps rows/columns reference
A
# A B C D E
#A 0 1 2 0 4
#B 1 0 0 2 0
#C 2 0 0 1 0
#D 0 2 1 0 3
#E 4 0 0 3 0
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
119 次 |
| 最近记录: |