将嵌套求和转换为R型矩阵乘法
我尝试在R中实现嵌套求和。for循环实现为:将嵌套求和转换为R型矩阵乘法,r,nested,matrix-multiplication,R,Nested,Matrix Multiplication,我尝试在R中实现嵌套求和。for循环实现为: sum = 0 for(i in 1:n){ for(j in 1:n){ for(k in 1:n){ sum = sum + w[i,j]*w[j,k] } } } 其中,w为对称方阵,n为行数(或列数)。 我正在努力实施。(因此,我不允许在此处书写Latex或添加图像。) 上面的嵌套for循环会花费很长时间。如何以R方式有效地实现此功能?尝试以下方法: Sum2 <- s
sum = 0
for(i in 1:n){
for(j in 1:n){
for(k in 1:n){
sum = sum + w[i,j]*w[j,k]
}
}
}
其中,w
为对称方阵,n
为行数(或列数)。
我正在努力实施。(因此,我不允许在此处书写Latex或添加图像。)
上面的嵌套for循环会花费很长时间。如何以R方式有效地实现此功能?尝试以下方法:
Sum2 <- sum(w %*% w)
all.equal(Sum, Sum2)
## [1] TRUE
Sum2试试这个:
Sum2 <- sum(w %*% w)
all.equal(Sum, Sum2)
## [1] TRUE
Sum2您能提供一个示例和期望的输出吗?您能提供一个示例和期望的输出吗?