我试图找到一种更有效的方式来编写这段代码, 我认为应用,mapply和扫描,但我认为没有办法重写它...
points.proj是m x k矩阵,data.proj是n x k矩阵。
所以基本上我想用尖端对points.mat的每个元素应用fun。其他矩阵data.mat ...的相同列数也应该是m×k矩阵。
for(i in 1:m){
for(j in 1:k){
Bounds[i,j] <- fun(points.proj[i,j],data.proj[,j])
}}
fun <- function(a,b) sum(a<b)
> points.proj
[,1] [,2]
[1,] 6 5
[2,] 7 6
[3,] 8 5
> data.proj
[,1] [,2]
[1,] 8 3
[2,] 2 0
[3,] 9 4
[4,] 6 7
[5,] 2 9
> Bounds
[,1] [,2]
[1,] 2 2
[2,] 2 2
[3,] 1 2
感谢您的帮助