r - 有效地求和通过不同矩阵关联的整个行

Question

我的代码的目标是创建一个新的矩阵BySite，其中的行是LOLZ行的总和（我正在使用colSums()），这些行与XY矩阵相关。因此，如果两个XY行匹配，那么我希望新矩阵具有一行，该行将相关的（根据XY）LOLZ行相加。

BySite<-matrix(0,dim(unique(XY))[1],dim(LOLZ)[2]);
BySite[1,]<-LOLZ[1,];
for (i in 1:dim(unique(XY))[1]){
  tree<-which(apply(XY, 1, function(x) all(x == unique(XY)[i,]))==TRUE)
  if (length(tree)==1){
    BySite[i,]<-LOLZ[tree,]
  }
  else {
    BySite[i,]<-colSums(LOLZ[tree,])
  }
cat(i," *sigh*")
}

XY是纬度和经度矩阵，因此可以使用：

XY<-matrix(c(rep(c(55:74),100),rep(c(101:120),100)),2000,2)

对于LOLZ

 LOLZ<-matrix(abs(rnorm(2000*50)),2000,50)

我的代码有效，但速度很慢，正在寻找一些建议以使其更快。

Answer 1

这是尝试使用较小的数据集：

XY = cbind(rep(50:55, 3), rep(100:105, 3))
set.seed(007); LOLZ = matrix(sample(1:5, 18 * 5, T), 18, 5)

paste_XY = paste(XY[, 1], XY[, 2], sep = "; ") #or apply(XY, 1, paste, collapse = "; ")
ans = rowsum(LOLZ, paste_XY)

#after running your code to build "BySite"
sum(ans != BySite)
#[1] 0