У меня есть много подматриц, которые я хотел бы собрать вместе, чтобы создать только одну матрицу.
Цель этого кода - разделить основную матрицу на подматрицы и вычислить пропорциональность между каждой из них.эти матрицы.
Я получил много подматриц, но я не знаю, как собрать их обратно.
Приветствия
Это мой код:
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
}