双变量求和

Question

所以一段时间以来，我一直在尝试在RStudio上进行两次求和，但进展却很快。...有人知道如何求解这些方程式吗？

到目前为止，我为第一次求和所做的代码如下：

IndexStart = 1

i = seq(IndexStart, 17, 1)

j = seq(IndexStart, 13, 1)

resultb = sum(i*j)

print(resultb)

第二，由于pi和pj具有不同的值，因此我想将p1和p2列为单独的变量，也许使用上面的代码样式？

任何输入都会有所帮助

Answer 1

可以在outer的帮助下计算第一次求和。

i <- 1:17
j <- 1:13
sum(outer(i, j))
#[1] 13923

对第二个求和进行类似的操作。
由于您没有发布测试数据，因此我将创建测试数据。

set.seed(1)    # Make the rsults reproducible
p.i <- runif(2)
p.j <- runif(2)

p.ij <- outer(p.i, p.j)
logp.ij <- log(p.ij)
sum(p.ij*logp.ij)
#[1] -1.325546