[R] how to draw crossing circles to show relationship?

[Guang Dai]

Dear list,
To make it easier, let me start with an example. 
Say I have 20 different cookies. Dave likes 7 of them, Evan likes
14 of them and Torri likes 11 of them. The data is

cookies<-cbind(c(1:20), c(1,0,0,1,0,1,0,0,0,0,1,1,1,0,0,0,1,0,0,0),
c(1,1,1,0,0,1,1,1,1,1,0,1,1,0,1,0,1,0,1,1),
c(1,1,0,1,0,0,1,0,1,1,0,0,1,0,1,0,1,1,1,0))
colnames(cookies)<-c("cookie","Dave","Evan","Torri")

Obviously, cookie 1,13 and 17 are liked by all of the person. cookie 2,7,9,10,13,15,17 and 19 are liked by Evan and Torri, etc. 

I want a chart by using  crossing circles to show how many is liked to all of them, how many is like by Dave and Evan,..., and how many is like only by Dave,.... Just like the Olympic log but with numbers to show the relationship. Of course, each circle represents one person in our data.
Better if we can draw the circles according to the total cookies liked by its corresponding person. 

I'm very new to R and have no idea how to do it. Thanks for any advices. 
regards,
Guang