ログインでの二重合計R

Question

このRの問題を解決するのに誰も助けてくれますか？私はログ演算子で二重合計を実行する必要があり、私はそれを行う方法を見つけることができません。この合計は、確率分布の対数尤度式の一部である。私はこれが最初は簡単な質問であることを知っていますが、私は解決なしでそれに時間を費やしています。前もって感謝します。

総和式は：と思う

これは、合計するデータを準備するためにRコードのチャンクです（）。これについて

set.seed(123) 
y <- rpois(100, 5) 
a <- 5; b <- 1; p = 0.1 
# ni <- frequency of repetitions from 100 Poisson distributed numbers 
ni <- table(y);ni 
i <- 2:length(ni) 
j <- 1:(length(i)-1) 
#Sum <- sum(????) 

Answer 1

set.seed(123) 
y <- rpois(100, 5) 
a <- 5; b <- 1; p = 0.1 
n <- as.vector(table(y)) 


result_list <- sapply(2:length(ni), function(i) { 
    sapply(1:(i-1), function(j) { 
     n[i]*log(a + b*i+p*j) 
    }) 
}) 

result <- sum(unlist(result_list)) 

どのように？