Maybe this example can help you. library(limma) hsb2<-read.table("http://www.ats.ucla.edu/stat/R/notes/hsb2.csv", sep=',', header=T) attach(hsb2) hw<-(write>=60) hm<-(math >=60) hr<-(read >=60) c3<-cbind(hw, hm, hr) a <- vennCounts(c3) vennDiagram(a) Best, Enrique Meza //