请教高手如何在R中快速的计算下面的求和

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问有没有人知道如何避免循环来计算$$\sum_{1\leq i<j<k\leq n}a_{ij}a_{ik}$$
其中$\mathbf{A}=(a_{ij})_{1\leq i,j\leq n}$是对称阵，即$a_{ij}=a_{ji}$.

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：请教高手 请问有没有 有没有人 Sum 有没有 如何

你试试这个 可以不可以

1. dat=function(data,i,j){
   
2. sum(t(sapply(1:i,function(i)combn(data[i,1:j],2,FUN=function(x)prod(x)))))
   
3. }
   
4. data=matrix(1:16,4,4)                       
   
5. dat(data,3,3)                        

复制代码

楼上厉害！学到东西了。我run了下，似乎有点小问题。修改了下：

dat=function(data,i,j){
    sum(unlist(sapply(1:i,function(i)combn(data[i,(i+1):j],2,FUN=function(x)prod(x)))))
  }                       

n<- 5
data=t(matrix(1:n^2,n,n))
dat(data,nrow(data)-2,ncol(data))

xinxin329 发表于 2014-6-17 18:49
楼上厉害！学到东西了。我run了下，似乎有点小问题。修改了下：

dat=function(data,i,j){

哈哈 其实我也想提醒下楼主的  1<=i < j < k<n 让楼主自己修改下标问题。
数学上面    我不专业

非常感谢楼上两位