ID の 2 つの文字ベクトルがあります。
2 つの文字ベクトルを比較したいと思います。特に、次の図に興味があります。
- A と B の両方にある ID の数
- A にあるが B にはない ID の数
- B にあるが A にはない ID の数
ベン図も描いてみたいです。
ID の 2 つの文字ベクトルがあります。
2 つの文字ベクトルを比較したいと思います。特に、次の図に興味があります。
ベン図も描いてみたいです。
以下に、試してみる基本事項をいくつか示します。
> A = c("Dog", "Cat", "Mouse")
> B = c("Tiger","Lion","Cat")
> A %in% B
[1] FALSE TRUE FALSE
> intersect(A,B)
[1] "Cat"
> setdiff(A,B)
[1] "Dog" "Mouse"
> setdiff(B,A)
[1] "Tiger" "Lion"
同様に、次のように簡単にカウントを取得できます。
> length(intersect(A,B))
[1] 1
> length(setdiff(A,B))
[1] 2
> length(setdiff(B,A))
[1] 2
私は通常、大規模なセットを扱っているので、ベン図の代わりに表を使用します。
xtab_set <- function(A,B){
both <- union(A,B)
inA <- both %in% A
inB <- both %in% B
return(table(inA,inB))
}
set.seed(1)
A <- sample(letters[1:20],10,replace=TRUE)
B <- sample(letters[1:20],10,replace=TRUE)
xtab_set(A,B)
# inB
# inA FALSE TRUE
# FALSE 0 5
# TRUE 6 3
intersectとsetdiffの代わりに、共通要素の%in%と boolean ベクトルを使用する別の方法です。2 つのリストではなく 2つのベクトルを実際に比較したいと思います。リストは任意の型の要素を含む可能性のある R クラスですが、ベクトルには常に 1 つの型の要素のみが含まれるため、真に等しいものを簡単に比較できます。ここでは、要素が文字列に変換されます。これは、存在していた最も柔軟性のない要素タイプであったためです。
first <- c(1:3, letters[1:6], "foo", "bar")
second <- c(2:4, letters[5:8], "bar", "asd")
both <- first[first %in% second] # in both, same as call: intersect(first, second)
onlyfirst <- first[!first %in% second] # only in 'first', same as: setdiff(first, second)
onlysecond <- second[!second %in% first] # only in 'second', same as: setdiff(second, first)
length(both)
length(onlyfirst)
length(onlysecond)
#> both
#[1] "2" "3" "e" "f" "bar"
#> onlyfirst
#[1] "1" "a" "b" "c" "d" "foo"
#> onlysecond
#[1] "4" "g" "h" "asd"
#> length(both)
#[1] 5
#> length(onlyfirst)
#[1] 6
#> length(onlysecond)
#[1] 4
# If you don't have the 'gplots' package, type: install.packages("gplots")
require("gplots")
venn(list(first.vector = first, second.vector = second))
前述のように、R でベン図をプロットするには複数の選択肢があります。gplots を使用した出力を次に示します。
With sqldf: Slower but very suitable for data frames with mixed types:
t1 <- as.data.frame(1:10)
t2 <- as.data.frame(5:15)
sqldf1 <- sqldf('SELECT * FROM t1 EXCEPT SELECT * FROM t2') # subset from t1 not in t2
sqldf2 <- sqldf('SELECT * FROM t2 EXCEPT SELECT * FROM t1') # subset from t2 not in t1
sqldf3 <- sqldf('SELECT * FROM t1 UNION SELECT * FROM t2') # UNION t1 and t2
sqldf1 X1_10
1
2
3
4
sqldf2 X5_15
11
12
13
14
15
sqldf3 X1_10
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
上記の回答の 1 つと同じサンプル データを使用します。
A = c("Dog", "Cat", "Mouse")
B = c("Tiger","Lion","Cat")
match(A,B)
[1] NA 3 NA
この関数は、 内のすべての値match
の位置を含むベクトルを返します。したがって、の 2 番目の要素は の 3 番目の要素です。他の試合はありません。B
A
cat
A
B
A
とで一致する値を取得するには、次のB
ようにします。
m <- match(A,B)
A[!is.na(m)]
"Cat"
B[m[!is.na(m)]]
"Cat"
A
との一致しない値を取得するにはB
:
A[is.na(m)]
"Dog" "Mouse"
B[which(is.na(m))]
"Tiger" "Cat"
さらに、length()
一致する値と一致しない値の合計数を取得するために使用できます。