如何将子矩阵汇总成一个大矩阵?

时间:2019-04-03 09:40:15

标签: r matrix

为了创建一个矩阵,我想将许多子矩阵聚集在一起。

此代码的目标是将主矩阵划分为子矩阵,并计算每个矩阵之间的比例。

我有很多子矩阵,但我不知道如何将它们重新聚集在一起。

欢呼

这是我的代码:

data <- matrix(sample(1002*1002), 1002, 1002)

ncol<-ncol(data)
rest<-ncol%%100
large<-ncol-rest
blocks<-ncol%/%100 


ngroup <- rep(1:blocks, each = 100)

if (rest > 0) ngroup <- c(ngroup, rep(blocks+ 1, rest))

split <- split(1:ncol, ngroup)


#that gives me all the pairwise comparison between the blocks .



combs <- expand.grid(1:length(split), 1:length(split)) 
combs <- t(apply(combs, 1, sort))
combs <- unique(combs) 



for(i in 1:nrow(combs)){

  batch1 <- split[[combs[i,1]]]
  batch2 <- split[[combs[i,2]]]
  dat.i <- cbind(data[,batch1], data[,batch2])
  rho.i <- cor_rho(dat.i) ## cor_rho is my function
}

0 个答案:

没有答案