小白求教：计算距离的一个for循环老是显示脚标有误

1. for(i in 1 : nrow(a)){
   
2.      for(j in 1 : nrow(b)){
   
3.           #让b的每一条数据轮流与a计算距离
   
4.           a_b[i,j] <- as.matrix(dist(rbind(a[i,], b[j,]), method = "euclidean"))
   
5.    }
   
6. }

复制代码

菜鸟求助a是50条数据（鸢尾花的数据），b也是50条数据，想让a的每条数据与b的每条数据算一下距离，然后生成一个矩阵a_b，例如a的第一条数据和b的第一条数据距离就是a_b[1,1]，a的第二条数据与b的第一条数据的距离就是a_b[2,1]。但是现在报误：incorrect number of subscripts on matrix，不知道如何修改，请大神赐教

下面给了我的方法和你的方法，结果的矩阵都是第m行第n列的数值=a的第n行和b的第m行的距离

1. mydist <- function(a, b) {
   
2.   apply(a, 1, function(vec, mat) {
   
3.     apply(mat, 1, function(x){
   
4.       dist(rbind(x, vec))
   
5.     })
   
6.   }, mat = b)
   
7. }
   
8. 
   
9. yourdist <- function(a, b) {
   
10.   result <- matrix(NA, nrow = nrow(b), ncol = nrow(a))
    
11.   for (i in 1:nrow(a))
    
12.     for (j in 1:nrow(b))
    
13.       result[j, i] <- dist(rbind(a[i, ], b[j, ]))
    
14.   return(result)
    
15. }

复制代码

推荐mydist，速度应该能快10倍左右

1. > n <- 100
   
2. > a <- data.frame(x = runif(n), y = runif(n))
   
3. > b <- data.frame(x = runif(n), y = runif(n))
   
4. > system.time(mydist(a, b))
   
5. 用户 系统 流逝
   
6. 0.25 0.00 0.25
   
7. > system.time(yourdist(a, b))
   
8. 用户 系统 流逝
   
9. 3.34 0.00 3.50

复制代码

a_b[i,j]里面应该是个数值，不用as.matrix吧