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

2531

收藏 2014-06-17

请问有没有人知道如何避免循环来计算$$\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}$.

扫码加我拉你入群

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

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

全部回复

jmpamao

2014-6-17 16:17:47

你试试这个可以不可以

复制代码

扫码加我拉你入群

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

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

xinxin329

2014-6-17 18:49:18

楼上厉害！学到东西了。我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))

扫码加我拉你入群

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

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

jmpamao

2014-6-17 18:54:04

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

dat=function(data,i,j){

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

扫码加我拉你入群

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

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

楚韵荆风

2014-6-17 21:08:08

非常感谢楼上两位

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群