Найдите строки матрицы, содержащие строки другой матрицы.

Question

Предположим, что у меня есть матрица m x n M1 и k x l матрица M2 с l < = n. Я хочу найти все эти строки в M1, которые содержат в себе некоторую строку M2.

Для примера рассмотрим следующую ситуацию:

> M1 <- matrix(c(1,2,3,4,5,6,7,8,9), nrow = 3, ncol = 3, byrow = TRUE) 
> M2 <- matrix(c(1,3,8,9), nrow = 2, ncol = 2, byrow = TRUE) 
> M1 
     [,1] [,2] [,3] 
[1,] 1 2 3 
[2,] 4 5 6 
[3,] 7 8 9 
> M2 
     [,1] [,2] 
[1,] 1 3 
[2,] 8 9 

Затем ряды один и три М1 выполнить условие, в строке один содержит 1 и 3 и последней строки 8 и 9.

Так как для достижения этого эффективным способом? Я написал код с использованием циклов, но поскольку я работаю с очень большими матрицами, это решение занимает много времени.

Answer 1

Более общий пример:

M1 <- matrix(c(1,2,3,1,2,3,4,5,6,7,8,9,10,11,12,13,1,2), nrow = 6, ncol = 3, byrow = TRUE) 
M2 <- matrix(c(1,2,6,9, 10,11,16,17, 19, 2), nrow = 5, ncol =2, byrow = TRUE) 

Во-первых, использовать match, чтобы найти индексы значений соответствия в M1.

ind <- match(M1, M2) 

Теперь, используя моды оператора %% с индексами и количеством строк, вы найдете строки. Это работает, потому что индексы всегда будут содержать строку M2, а также общее количество строк, поэтому числа в одной строке возвратят тот же результат.

rows <- ind %% nrow(M2) 

Затем m представляет собой матрицу, содержащую номер строки совпадающих значений между M1 и M2. Линии будут выбраны, только если один и тот же индекс появится в той же строке 2 раза (или, более общо, число раз, равное числу столбцов в M2). Это гарантирует, что строка M1 будет рассмотрена только в том случае, если она содержит все элементы строки в M2.

m <- matrix(rows, nrow = nrow(M1)) 
matchRows <- apply(m, 1, duplicated, incomparables = NA) 
M1rows <- which(colSums(matchRows)==ncol(M2)-1) 

Answer 2

Этот метод будет проверять каждую строку из М2 и возвращает индекс строки из M1, если он содержится или НС в случае, если это не

M1 <- matrix(c(1,2,3,4,5,6,7,8,9), nrow = 3, ncol = 3, byrow = TRUE) 
> M1 
    [,1] [,2] [,3] 
[1,] 1 2 3 
[2,] 4 5 6 
[3,] 7 8 9 

M2 <- matrix(c(1,3,8,5,4,5,1,2), nrow = 4, ncol = 2, byrow = TRUE) 
> M2 
    [,1] [,2] 
[1,] 1 3 
[2,] 8 5 
[3,] 4 5 
[4,] 1 2 

y = apply(M2,1,function(x){ 
    z = unique(which(M1 %in% x)%%nrow(M1)) 
    ifelse(length(z)==1,ifelse(z==0,nrow(M1),z),NA) 
}) 

> y 
[1] 1 NA 2 1 

Это означает, что строка 2 от М2 не в M1, а строки 1 и 4 из M2 находятся в строке 1 в M1. Также строка 3 в M2 находится в строке 2 в M1.