Dear r-helpers,
Here I would like to have your kind helps on generating Venn diagram.
There are some packages within R on this task, like venneuler, VennDiagram,
vennerable. But, vennerable can not be installed on my Mac book. It seems
VennDiagram can not work on my data. And, venneuler may have generated a
wrong Venn diagram to me.
Do you have any experience/expertise on those Venn diagram? Could you please
give me any directions on that?
Thanks in advance.
Best wishes,
Jian-Feng,
######################################################################
# (1) my code for venneuler
vd <-
venneuler(c("A&B&C&D&E&F"69604,"A&B&C&D&E"=426120
,"A&B&C&D&F"=20297,"A&B&C&D"=123063,"A&B&C&E&F"=12695,"A&B&C&E"=115100,"A&B&C&F"=11667,"A&B&C"=95656,"A&B&D&E&F"=1755,"A&B&D&E"=20113,"A&B&D&F"=1903,"A&B&D"=19218,"A&B&E&F"=2831,"A&B&E"=38362,"A&B&F"=4950,"A&B"=68289,"A&C&D&E&F"=11657,"A&C&D&E"=107235,"A&C&D&F"=14883,"A&C&D"=193338,"A&C&E&F"=6284,"A&C&E"=79985,"A&C&F"=14710,"A&C"271416
,"A&D&E&F"=1069,"A&D&E"=17628,"A&D&F"=3152,"A&D"=71573,"A&E&F"=2786,"A&E"=57511,"A&F"=13684,A475970
,"B&C&D&E&F"=2722,"B&C&D&E"=30528,"B&C&D&F"=2740,"B&C&D"=30986,"B&C&E&F"=3579,"B&C&E"=55443,"B&C&F"=7789,"B&C"=101005,"B&D&E&F"=917,"B&D&E"=14894,"B&D&F"=1436,"B&D"=24972,"B&E&F"=3975,"B&E"=105527,"B&F"16877,B=718570
,"C&D&E&F"=1587,"C&D&E"=26289,"C&D&F"=4902,"C&D"=101947,"C&E&F"=3326,"C&E"=77289,"C&F"=20125,C689330,"D&E&F"=892,"D&E"=22666,"D&F"=4661,D=200020,"E&F"=8518,E=521290,F401622))
pdf("myvenn.pdf")
plot(vd)
dev.off()
#####################################################################
# (2) the problem of the plot venneuler generated me is sets (A,B,C,D,E,F)
should shared 69604 elements.
# But, it illustrated nothing for me for this 6 sets sharing.
#####################################################################
# (3) I prepared my code for vennerable package, but it can not be installed
now.
myVenn <- Venn(SetNames =
c("Norway","Russia","Iceland",
"Scotland","Austria","North American"), Weight =
c('111111'=69604,'111110'426120
,'111101'=20297','111100'=123063,'111011'=12695,'111010'=115100,'111001'=11667,'111000'=95656,'110111'=1755,'110110'=20113,'110101'=1903,'110100'=19218,'110011'=2831,'110010'=38362,'110001'=4950,'110000'=68289,'101111'=11657,'101110'=107235,'101101'=14883,'101100'=193338,'101011'=6284,'101010'=79985,'101001'=14710,'101000'271416
,'100111'=1069,'100110'=17628,'100101'=3152,'100100'=71573,'100011'=2786,'100010'=57511,'100001'=13684,'100000'475970
,'011111'=2722,'011110'=30528,'011101'=2740,'011100'=30986,'011011'=3579,'011010'=55443,'011001'=7789,'011000'=101005,'010111'=917,'010110'=14894,'010101'=1436,'010100'=24972,'010011'=3975,'010010'=105527,'010001'16877,'010000'=718570
,'001111'=1587,'001110'=26289,'001101'=4902,'001100'=101947,'001011'=3326,'001010'=77289,'001001'=20125,'001000'=689330,'000111'=892,'000110'=22666,'000101'=4661,'000100'=200020,'000011'=8518,'000010'=521290,'000001'=401622)
pdf("myVenn")
plot(myVenn, doWeight = T, type = "circles")
dev.off()
[[alternative HTML version deleted]]