set.seed(123)
B = matrix( c(5, 3, 3, 1, 5, 1,3,1,NA,NA), nrow=5, ncol=2)
m1<-matrix(nrow=5,ncol=2,data=runif(10))
m2<-matrix(nrow=5,ncol=2,data=runif(10))
m2[1,2]=NA; ml <- list(m1, m2)
ind <- sapply(unique(c(B)), function(x) which(B == x, arr.ind = TRUE))
re <- lapply(ind, function(x) lapply(ml, function(y) y[x]))
res=lapply(re, function(x) c(t(do.call(cbind, x))))
Run Code Online (Sandbox Code Playgroud)
但我不知道哪个对应哪个.例如:res[[1]]代表5, 3, or 1 in B什么?有没有办法t从B的通信类(数字)中命名输出colms ?
我们可以用 complete.cases
lapply(re, function(x) {
v1 <- c(t(do.call(cbind, x)))
v1[complete.cases(v1)]})
Run Code Online (Sandbox Code Playgroud)