2

如果我有以下数据:

mat1           <- matrix( c(0,2,3,1,0,1,1,1,1), nrow=3 )
rownames(mat1) <- LETTERS[1:3]
colnames(mat1) <- LETTERS[1:3]

mat1
#  A B C
#A 0 1 1
#B 2 0 1
#C 3 1 1

如何到达data.frame我知道哪些行矩阵条目大于列条目的地方?例如,我想要的结果mat1是某种data.frame指示:

A < B 
A < C 
B > A
B = C
C > A
C = B

我一直在 upper.tri使用减法来搞乱等,但这无济于事。

4

1 回答 1

5

要将每个元素mat1[i,j] 与对角线另一侧的对应元素进行比较mat1[j,i],可以使用转置t

r <- sign( mat1 - t(mat1) )
r
#   A  B  C
# A 0 -1 -1
# B 1  0  0
# C 1  0  0

paste( 
  colnames(r)[col(r)], 
  ifelse(r>0, ">", ifelse(r<0, "<", "=")), 
  rownames(r)[row(r)] 
)
# [1] "A = A" "A > B" "A > C" "B < A" "B = B" "B = C" "C < A" "C = B" "C = C"
于 2013-05-19T15:13:00.250 回答