Guest User

Untitled

a guest
Dec 10th, 2018
93
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.38 KB | None | 0 0
  1. library(VennDiagram)
  2.  
  3. BM1 <- c('Cistaceae', 'Oleaceae','Fagaceae','Rhamnaceae','Asteraceae','Asteraceae')
  4. BM2 <- c('Boraginaceae','Poaceae','Poaceae')
  5. BM3 <- c('Poaceae','Asteraceae','Chenopodiaceae','Asteraceae','Poaceae','Solanaceae','Solanaceae','Asteraceae','Rubiaceae')
  6. BM4<- c('Amaranthaceae','Chenopodiaceae','Convolvulaceae','Convolvulaceae','Cucurbitaceae','Cucurbitaceae','Cucurbitaceae','Cucurbitaceae','Boraginaceae','Boraginaceae','Solanaceae')
  7.  
  8. v2 <- venn.diagram(list(Oak=BM1, Wasteland=BM2,Edge=BM3, Crop=BM4),
  9. fill = c("red", "blue","green","purple"),
  10. alpha = c(0.5, 0.5, 0.5,0.5), cat.cex = 1.5, cex=1.5,
  11. filename=NULL)
  12. grid.newpage()
  13. grid.draw(v2)
  14.  
  15.  
  16. v2[[7]]$label <- paste(setdiff(BM1,union(BM2,BM3)),collapse = "")
  17. inters <- intersect(BM1,intersect(BM2,BM3))
  18. v2[[8]]$label <- paste(setdiff(intersect(BM1,BM2), inters),collapse = "")
  19. v2[[9]]$label <- paste(setdiff(BM2,union(BM1,BM3)),collapse = "")
  20. v2[[10]]$label <- paste(setdiff(intersect(BM1,BM3), inters),collapse = "")
  21. v2[[11]]$label <- paste(ww[[1,5]], collapse = "")
  22. v2[[12]]$label <- paste(setdiff(intersect(BM3,BM2), inters),collapse = "")
  23. v2[[13]]$label <- paste(setdiff(BM3,union(BM1,BM2)),collapse = "")
  24.  
  25.  
  26. grid.newpage()
  27. grid.draw(v2)
  28.  
  29. library(VennDiagram)
  30.  
  31. BM1 <- c('Cistaceae', 'Oleaceae','Fagaceae','Rhamnaceae','Asteraceae','Asteraceae')
  32. BM2 <- c('Boraginaceae','Poaceae','Poaceae')
  33. BM3 <- c('Poaceae','Asteraceae','Chenopodiaceae','Asteraceae','Poaceae','Solanaceae','Solanaceae','Asteraceae','Rubiaceae')
  34. BM4<- c('Amaranthaceae','Chenopodiaceae','Convolvulaceae','Convolvulaceae','Cucurbitaceae','Cucurbitaceae','Cucurbitaceae','Cucurbitaceae','Boraginaceae','Boraginaceae','Solanaceae')
  35.  
  36. lista_conjuntos <- list(Oak=BM1, Wasteland=BM2,Edge=BM3, Crop=BM4)
  37.  
  38. v2 <- venn.diagram(lista_conjuntos,
  39. fill = c("red", "blue","green","purple"),
  40. alpha = c(0.5, 0.5, 0.5,0.5), cat.cex = 1.5, cex=.8,
  41. filename=NULL)
  42.  
  43. labels <- lapply(calculate.overlap(lista_conjuntos), function(x) {paste(unique(x), collapse='n')})
  44.  
  45. tail(labels, 3)
  46. $`a14`
  47. [1] ""
  48.  
  49. $a1
  50. [1] "Rubiaceae"
  51.  
  52. $a3
  53. [1] "AmaranthaceaenConvolvulaceaenCucurbitaceae"
  54.  
  55. min(which(unlist(lapply(v2, function(x) {class(x)[[1]]}=="text"))))
  56. [1] 9
  57.  
  58. for (x in names(labels)) {
  59. v2[[ start + as.integer(gsub("a", "",x)) ]]$label = labels[[x]]
  60. }
  61.  
  62. grid.newpage()
  63. grid.draw(v2)
Add Comment
Please, Sign In to add comment