如何将 R 中的两个总和更改为两个 for 循环？

Question

编辑：

我有一个名为y的向量，它有n数字元素和一个名为w的矩阵n*n ，其中它的元素是数字。

y<-c(312,218,259,76,263) 

> w
     [,1] [,2] [,3] [,4] [,5]
[1,]    0    1    0    1    1
[2,]    1    0    1    1    1
[3,]    0    1    0    0    1
[4,]    1    1    0    0    0
[5,]    1    1    1    0    0

我想对上述数据使用以下公式： 我使用函数 external 和sum在outer中编写了以下代码：

dy<-y-mean(y)
M<-outer(dy,dy,"*")
a<-sum(w*M)
s0<-sum(w)
n<-length(y)
v<-n/sum(dy^2)
L<-v*a/s0

L的output为-0.1873

现在，我想将上面的代码更改为使用两个for循环编写它（不使用outer和sum ）。 实际上，我想使用两个for循环而不是使用-0.1873和sum来制作公式，然后再次找到L as outer 。

预先感谢您的帮助。

Answer 1

首先计算该公式的所有分子和分母，然后从中计算 I。

这应该可以解决问题：

##Fraction: n1/d1 * n2/d2

d1 = 0
n1 = n
d2 = 0
n2 = 0
my = mean(y)
for(i in 1:n){
  d1 = d1 + (my-y[i])**2
  for(j in 1:n){
    n2 = n2+w[i,j]*(y[i]-my)*(y[j]-my)
    d2 = d2+w[i,j]
  }
}

I = n1/d1 * n2/d2