为了创建一个矩阵,我想将许多子矩阵聚集在一起。
此代码的目标是将主矩阵划分为子矩阵,并计算每个矩阵之间的比例。
我有很多子矩阵,但我不知道如何将它们重新聚集在一起。
欢呼
这是我的代码:
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
}