有效的方法来做这个双重求和？

Question

我对以下代码的效率有疑问：

g=0

for (i in 1:n){

  for (j in 1:m){

    g=g+(as.numeric(x[i]>y[j])*as.numeric(z[i]<=a)*as.numeric(h[j]>a))

  }
}

感谢您的帮助。

Answer 1

t(h > a) %*% sapply(x, ">", y) %*% (z <= a)