подмножество на матрица в for цикъл с логически оператор или (|)

Опитвам се да направя подгрупи на матрица. Всяко подмножество зависи от набор от логически условия. Например първото подмножество на първия списък (R[[1]][[1]]) съответства на последната колона на матрицата, равна на 1 или 3 и т.н. Бих искал да използвам това в цикъл.

Това е пример за данните, които имам:

l = list(
  list(c(1, 3),c(2, 3)),
  list(c(1,5),c(2,5,6,7),c(3,6,8),c(4,7,8))
)

a1 <- array(1, dim=c(12,10))
a2 <- array(1, dim=c(18,10))
v1 <- c(1,1,3,1,2,3,1,1,3,1,2,2)
v2 <- c(1,1,3,1,5,2,4,5,5,5,7,8,1,4,5,8,2,4)<br/>
m1 <- as.matrix(cbind(a1,v1))
m2 <- as.matrix(cbind(a2,v2))
m <- list(m1,m2)

Това е кодът, който използвам:

R <- rep(list(list()),length(m))
  for(i in 1:length(m)) {
    for(j in 1:length(l[[i]])) {
      R[[i]][[j]] <- m[[i]][m[[i]][,dim(m[[i]])[2]] == l[[i]][[j]],]
}}

Получавам следните резултати, които не са това, което искам:

[[1]]
[[1]][[1]]
[1,] 1 1 1 1 1 1 1 1 1 1  1
[2,] 1 1 1 1 1 1 1 1 1 1  3
[3,] 1 1 1 1 1 1 1 1 1 1  1
[[1]][[2]]
[1,] 1 1 1 1 1 1 1 1 1 1  2
[2,] 1 1 1 1 1 1 1 1 1 1  3
[3,] 1 1 1 1 1 1 1 1 1 1  2
[[2]]
[[2]][[1]]
[1,] 1 1 1 1 1 1 1 1 1 1  1
[2,] 1 1 1 1 1 1 1 1 1 1  5
[3,] 1 1 1 1 1 1 1 1 1 1  5
[4,] 1 1 1 1 1 1 1 1 1 1  1
[[2]][[2]]
[1,] 1 1 1 1 1 1 1 1 1 1  5
[2,] 1 1 1 1 1 1 1 1 1 1  2
[[2]][[3]]
 1  1  1  1  1  1  1  1  1  1  8 
[[2]][[4]]
[1,] 1 1 1 1 1 1 1 1 1 1  4
[2,] 1 1 1 1 1 1 1 1 1 1  7
[3,] 1 1 1 1 1 1 1 1 1 1  8

Това е резултатът, който бих искал да получа вместо това:

[[1]]
[[1]][[1]]
 [1,] 1 1 1 1 1 1 1 1 1 1  1
 [2,] 1 1 1 1 1 1 1 1 1 1  1
 [3,] 1 1 1 1 1 1 1 1 1 1  3
 [4,] 1 1 1 1 1 1 1 1 1 1  1
 [5,] 1 1 1 1 1 1 1 1 1 1  3
 [6,] 1 1 1 1 1 1 1 1 1 1  1
 [7,] 1 1 1 1 1 1 1 1 1 1  1
 [8,] 1 1 1 1 1 1 1 1 1 1  3
 [9,] 1 1 1 1 1 1 1 1 1 1  1
[[1]][[2]]
[1,] 1 1 1 1 1 1 1 1 1 1  3
[2,] 1 1 1 1 1 1 1 1 1 1  2
[3,] 1 1 1 1 1 1 1 1 1 1  3
[4,] 1 1 1 1 1 1 1 1 1 1  3
[5,] 1 1 1 1 1 1 1 1 1 1  2
[6,] 1 1 1 1 1 1 1 1 1 1  2
[[2]]
[[2]][[1]]
 [1,] 1 1 1 1 1 1 1 1 1 1  1
 [2,] 1 1 1 1 1 1 1 1 1 1  1
 [3,] 1 1 1 1 1 1 1 1 1 1  1
 [4,] 1 1 1 1 1 1 1 1 1 1  5
 [5,] 1 1 1 1 1 1 1 1 1 1  5
 [6,] 1 1 1 1 1 1 1 1 1 1  5
 [7,] 1 1 1 1 1 1 1 1 1 1  5
 [8,] 1 1 1 1 1 1 1 1 1 1  1
 [9,] 1 1 1 1 1 1 1 1 1 1  5
[[2]][[2]]
[1,] 1 1 1 1 1 1 1 1 1 1  5
[2,] 1 1 1 1 1 1 1 1 1 1  2
[3,] 1 1 1 1 1 1 1 1 1 1  5
[4,] 1 1 1 1 1 1 1 1 1 1  5
[5,] 1 1 1 1 1 1 1 1 1 1  5
[6,] 1 1 1 1 1 1 1 1 1 1  7
[7,] 1 1 1 1 1 1 1 1 1 1  5
[8,] 1 1 1 1 1 1 1 1 1 1  2
[[2]][[3]]
[1,] 1 1 1 1 1 1 1 1 1 1  3
[2,] 1 1 1 1 1 1 1 1 1 1  8
[3,] 1 1 1 1 1 1 1 1 1 1  8
[[2]][[4]]
[1,] 1 1 1 1 1 1 1 1 1 1  4
[2,] 1 1 1 1 1 1 1 1 1 1  7
[3,] 1 1 1 1 1 1 1 1 1 1  8
[4,] 1 1 1 1 1 1 1 1 1 1  4
[5,] 1 1 1 1 1 1 1 1 1 1  8
[6,] 1 1 1 1 1 1 1 1 1 1  4

r list for-loop logical-operators
person Charlie    schedule 09.08.2013    source източник

Отговори (2)


arrow_upward
1
arrow_downward

Търсите ли това?

Ако искате същото да приложи всички елементи от списък m:

k<-lapply(m,function(x) x[x[,11] %in% c(1,3),])

> k
[[1]]
                          v1
 [1,] 1 1 1 1 1 1 1 1 1 1  1
 [2,] 1 1 1 1 1 1 1 1 1 1  1
 [3,] 1 1 1 1 1 1 1 1 1 1  3
 [4,] 1 1 1 1 1 1 1 1 1 1  1
 [5,] 1 1 1 1 1 1 1 1 1 1  3
 [6,] 1 1 1 1 1 1 1 1 1 1  1
 [7,] 1 1 1 1 1 1 1 1 1 1  1
 [8,] 1 1 1 1 1 1 1 1 1 1  3
 [9,] 1 1 1 1 1 1 1 1 1 1  1

[[2]]
                         v2
[1,] 1 1 1 1 1 1 1 1 1 1  1
[2,] 1 1 1 1 1 1 1 1 1 1  1
[3,] 1 1 1 1 1 1 1 1 1 1  3
[4,] 1 1 1 1 1 1 1 1 1 1  1
[5,] 1 1 1 1 1 1 1 1 1 1  1
person Metrics    schedule 09.08.2013
comment
Благодаря за вашият отговор. Това не беше точно това, което търсех, но ми помогна да намеря решението, което публикувах по-долу. - person Charlie; 10.08.2013

arrow_upward
1
arrow_downward

Ето едно решение, което работи за мен с %in%:

R <- rep(list(list()),length(m)) <br/>
for(i in 1:length(m)) { <br/>
for(j in 1:length(l[[i]])) { <br/>
R[[i]][[j]] <- m[[i]][m[[i]][,dim(m[[i]])[2]] %in% l[[i]][[j]],] <br/>
}}
person Charlie    schedule 10.08.2013