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的长度。我认为这个表达式看起来不错，但是我的问题是，现在我可以加上第二个is..How吗？我不知道，因为有两个指数。你有什么建议吗？谢谢。

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中，使用向量化比循环更容易。