R中的双重求和

Question

我必须在R：

中实现这个总和

我试图通过以下方式实现第一个指数：

f <- rep(0,48)
for(i in 1:48) {
f[i] <- exp(-sum(lambda[which(lambda<lambda[i])])-(i-1))}

其中48是lambda的长度。我认为这个表达式看起来不错，但我现在的问题是......如何添加第二个求和？由于双重差异，我不知道任何想法。 你有什么建议吗？ 谢谢。

Answer 1

我不确定你的公式中有几个括号，但是类似的东西应该可以解决这个问题：

exp(-sum(lambda[20:j]) + (j-20)) * sum(unlist(lapply(20:j, function(i) lambda[i] * exp(sum(lambda[20:i]) + i-20) )))

在R中，使用矢量化比使用循环更容易。