[R] Proportional Venn diagram 6 sets with labels with pkg:venneuler
David Winsemius
dwinsemius at comcast.net
Sun Feb 23 17:57:31 CET 2014
On Feb 23, 2014, at 7:11 AM, Luigi Marongiu wrote:
> Dear all,
Dear Luigi;
Please stop posting triplicate postings and also learn to post in plain text.
>
> I would like to draw a Venn plot for data represented by 6 variables (or
> sets). I know how to do this using the package venneuler (which requires
> rJava). However this package does not report the numbers of elements within
> each class.
>
> Do you know an alternative to this package to draw this without having to
> calculate the area of all the possible combinations?
>
If you having difficulty calculating the counts in categories, then consider this:
> table(apply(data,1,function(x) paste(names(x[x==1]), collapse="_")))
a_e a_e_f a_f b_c_e b_d b_d_f b_e b_e_f b_f c_d
2 1 11 1 1 1 5 9 3 1
c_e c_e_f c_f d_e_f e e_f f
10 6 9 2 17 25 1
To answer your direct question, I do not know how to get venneuler to do that. But the maintainer of the package reads rhelp, so even though you should have posed the question to him first, you should eventually get his eyes on the question.
--
David.
> Best wishes,
>
> Luigi
>
>
>
> library(venneuler)
>
> library(rJava)
>
>
>
> a<-c( 1, 1, 1, 1, 1, 1,
> 1,
> 1, 1, 1, 1, 1, 1, 1,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0)
>
> b<-c( 0, 0, 0, 0, 0, 0,
> 0,
> 0, 0, 0, 0, 0, 0, 0,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0)
>
> c<-c( 0, 0, 0, 0, 0, 0,
> 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 1, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 1, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 1,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0)
>
> d<-c( 0, 0, 0, 0, 0, 0,
> 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 1, 1, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 1, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 1,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 1, 0,
> 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0)
>
> e<-c( 1, 1, 1, 0, 0, 0,
> 0,
> 0, 0, 0, 0, 0, 0, 0,
> 1, 1, 0, 1, 0, 0, 1,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 0, 0, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 0, 0,
> 0, 1, 0, 1, 1, 1, 1,
> 1, 1, 1, 1, 0, 0, 1,
> 1, 1, 1, 0, 0, 0, 0,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 0)
>
> f<-c( 0, 0, 1, 1, 1, 1,
> 1, 1, 1, 1, 1, 1, 1,
> 1, 1, 0, 1, 1, 0, 1,
> 1, 1, 0, 0, 1, 1, 0,
> 1, 0, 1, 1, 1, 1, 0,
> 0, 1, 0, 0, 0, 1, 1,
> 1, 1, 1, 0, 1, 0, 1,
> 0, 0, 1, 0, 1, 1, 1,
> 0, 0, 0, 0, 1, 1, 1,
> 1, 0, 1, 1, 1, 0, 1,
> 0, 0, 1, 1, 0, 0, 1,
> 1, 1, 0, 0, 0, 1, 1,
> 0, 1, 1, 1, 1, 1, 1,
> 0, 0, 1, 1, 1, 1, 1,
> 0, 0, 1, 0, 1, 1, 1,
> 1)
>
>
>
> data<-cbind(a,b,c,d,e,f)
>
>
>
> plot(venneuler(data), main="test")
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
David Winsemius
Alameda, CA, USA
More information about the R-help
mailing list