Jak znaleźć wspólne elementy z wielu wektorów?

159

Czy ktoś może mi powiedzieć, jak znaleźć wspólne elementy z wielu wektorów?

a <- c(1,3,5,7,9)
b <- c(3,6,8,9,10)
c <- c(2,3,4,5,7,9)

Chcę uzyskać wspólne elementy z powyższych wektorów (np: 3 i 9)

r vector r-faq

— Chares
źródło

45

Nie jest dobrym pomysłem używanie cjako nazwy zmiennej ...

— Marek

4

dlaczego to list jak inne?

— Mostafa

12

@DimitriPetrenko, ponieważ możesz deklarować listy za pomocą c(1,2...).

— Mathias711,

333

Może istnieć sprytniejszy sposób, aby to zrobić, ale

intersect(intersect(a,b),c)

wykona robotę.

EDYCJA: Sprytniej i wygodniej, jeśli masz dużo argumentów:

Reduce(intersect, list(a,b,c))

— bnaul
źródło

16

+1 za przypomnienie nam o Reducei poprawne użycie wielkich liter w R.

— mariotomo

8

Warto zauważyć, że dotyczy intersectto operacji na zbiorach. Jeśli masz elementy powtarzające się w wektorach, utracisz tę informację, ponieważ wektory zostaną zamienione w zbiory przed przecięciem. Np intersect(c(1,1,2,3), c(1,1,3,4))spowodowałoby c(1,3), a może nie chciał wynik c(1,1,3).

— Giora Simchoni

1

@GioraSimchoni, jak możesz uzyskać c (1,1,3), jeśli naprawdę tego chcesz?

— StatsSorceress

@StatsSorceress Załóżmy, że chcesz "przecięcie zachowujące duplikaty" wektorów składających się z dodatnich liczb całkowitych, wszystkie na liście L. Następujący kod działa:

N <- max(unlist(L)); LT <- lapply(L, tabulate, nbins = N); v <- do.call(pmin, LT); unlist(sapply(1:N, function(x) rep(x, v[x])))

Innym sposobem na zrobienie tego byłoby użycie matchfunkcji wraz z ujemnym indeksowaniem do iteracyjnego usuwania z każdego z wektorów każdy element dodany do "jądra".

— Montgomery Clift

24

To już dobra odpowiedź, ale można to zrobić na kilka innych sposobów:

unique(c[c%in%a[a%in%b]])

lub,

tst <- c(unique(a),unique(b),unique(c))
tst <- tst[duplicated(tst)]
tst[duplicated(tst)]

Można oczywiście pominąć uniquepołączeń, jeśli wiesz, że nie ma powtarzających się wartości ciągu a, blub c.

— James
źródło

7

intersect_all <- function(a,b,...){
  all_data <- c(a,b,...)
  require(plyr)
  count_data<- length(list(a,b,...))
  freq_dist <- count(all_data)
  intersect_data <- freq_dist[which(freq_dist$freq==count_data),"x"]
  intersect_data
}


intersect_all(a,b,c)

UPDATE EDIT Prostszy kod

intersect_all <- function(a,b,...){
  Reduce(intersect, list(a,b,...))
}

intersect_all(a,b,c)

— Abhishek K Baikady
źródło