図を作る

  • アメリカ人類遺伝学会のメモはこちらにするとして、お絵かきネタまではそちらに書かなくてもよいので、こちらのブログに
  • ベン図というのがあります。集合の図です。
  • 3群までのベン図は簡単
  • 学会では4群、5群が簡単そうに描かれていました
  • Rに関数があるはず…
library(gplots) # packageを読み込んで
# 適当な名前を作る関数を作って
oneName <- function() paste(sample(LETTERS,5,replace=TRUE),collapse="")
# 名前を1000個作ります
geneNames <- replicate(1000, oneName())
                       
## 
# 5群(A,B,C,D,E)に上で作った名前から要素数を指定してサンプリングします
# 複数の群に帰属する名前も生じます
GroupA <- sample(geneNames, 400, replace=FALSE)
GroupB <- sample(geneNames, 750, replace=FALSE)
GroupC <- sample(geneNames, 250, replace=FALSE)
GroupD <- sample(geneNames, 300, replace=FALSE)
GroupE <- sample(geneNames, 500, replace=FALSE)

input  <-list(GroupA,GroupB,GroupC,GroupD,GroupE)
input
# 帰属関係をベン図にしましょう
venn(input)
  • この方法で描けるのは、5群まで
  • それより多い群については、こんな感じで描図が話題になります