R tcrossprod具有独特的组合

Question

可能有一个版本的tcrossprod实现了这一点，但我无法找到它。从下面的示例中，如果订单无关紧要，如何只获得第一次出现组合的行？即（1,2）对我来说与（2,1）相同。

a <- c(1,2,3,4)
b <- c(10,5,4,10)
df<- data.frame(a,b)

melt(tcrossprod(df$b,1/df$b))

> melt(tcrossprod(df$b,1/df$b))
   Var1 Var2 value
1     1    1  1.00
2     2    1  0.50
3     3    1  0.40
4     4    1  1.00
5     1    2  2.00
6     2    2  1.00
7     3    2  0.80
8     4    2  2.00
9     1    3  2.50
10    2    3  1.25
11    3    3  1.00
12    4    3  2.50
13    1    4  1.00
14    2    4  0.50
15    3    4  0.40
16    4    4  1.00

Answer 1

拥有m <- melt(tcrossprod(df$b,1/df$b))，您只需执行以下操作：

subset(m,X1>X2)

#   X1 X2 value
#2   2  1   0.5
#3   3  1   0.4
#4   4  1   1.0
#7   3  2   0.8
#8   4  2   2.0
#12  4  3   2.5