r - 将向量中的唯一值对相乘并对结果求和

Question

我想乘以一个向量的唯一对，然后对其求和，不包括由相同元素组成的对，这样对于 c(1:4):

(1*2) + (1*3) + (1*4) + (2*3) + (2*4) + (3*4) == 35

下面的代码适用于上面的例子:

x <- c(1:4)
bar <- NULL
for( i in 1:length(x)) { bar <- c( bar, i * c((i+1) : length(x)))}
sum(bar[ 1 : (length(bar) - 2)])

但是，我的实际数据是有理数向量，而不是整数，因此循环的 (i+1) 部分将不起作用。有没有办法在 i 之后查看集合的下一个元素，例如j，这样我就可以写 i * c((j : length(x))?

我知道 for 循环通常不是最有效的方法，但我想不出如何通过 apply 等实现这一点。欢迎提供示例， 也。感谢您的帮助。

Answer 1

循环的替代方法是使用 combn 并使用 FUN 参数乘以组合。然后求和结果:

sum(combn(x = 1:4, m = 2, FUN = function(x) x[1] * x[2]))
# [1] 35

根据@bgoldst 的建议，在FUN 中使用prod 会更好:

sum(combn(x = 1:4, m = 2, FUN = prod))