Суммирование двойной переменной

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