我试图使用R来对具有三个索引的函数进行三次求和。
我可以使用以下代码轻松地在Mathematica中执行此操作：
out = Sum[B[G[[k]]] * A[G[[k]], G[[j]], G[[i]]] * prod[G[[j]],G[[i]]],
{k,1,Length[G]},{j,1,Length[G]},{i,1,Length[G]}]
其中G是一个矩阵，B[.]，A[.]和prod[.]都是预定义的函数。
在Mathematica中，Sum[f,{k,k_min, k_max}, {j,j_min, j_max}, {i, i_min, i_max}]将评估三重和，Sum（k = k_min，k_max）[Sum（j = j_min，j_max）[Sum（i = i_min，i_max）[f]]]，其中f是一种功能。
我现在正试图在R中这样做并且遇到了很多困难。我尝试按以下方式将sapply与sum一起使用，但它似乎也无法正常工作。
guts.i.j.k <- function(i,j,k) B.k(k) * A.i.j.k(i, j, k) * prod.i.j(i, j)
innermostSum <- function(j,k) sum(sapply(1:length(G), FUN=guts.i.j.k, j=j, k=k ))
middleSum <- function(k) sum(sapply(1:length(G), FUN=innermostSum, k=k ))
outsideSum <- sum(sapply(1:length(G), FUN=middleSum))
return(outsideSum)
非常感谢任何帮助。