In a message dated 7/7/2004 9:42:03 PM Pacific Daylight Time,
[EMAIL PROTECTED] writes:
Dear Alec,
Thank you for your response and it worked. Further, I
have another mathematical problem. I apologize ahead
as this question is not apt for this list. I am a
biologist working at Johns Hopkins School of Medicine.
As I listed in my previous e-mail (attached below)I
mentioned I have 4 protein sets. I am now trying to
calculate the combinatorics of these sets. My
ultimate aim is to draw a venn diagram and find out
the proteins that are unique to set C,S,E and H.
I drew a venn diagram graph and I am banging my head
to deduce the combinations. It is easy for me to
deduce the intersections - that means the protein
entries that are present in common. However, it proved
very difficult to deduce the following:
I could calculate the following:
NP_*** present in both C and E (C intersection E)
NP_*** present in both C and H (C intersection H)
NP_*** present in both C and S (C ^ S)
NP_*** present in both E and H (E ^ H)
NP_*** present in both E and S (E ^ S)
NP_*** present in both H and S (H ^ S)
NP_*** present in C, E and H (C^E^H)
NP_*** present in C, H and S (C^H^S)
NP_*** present in E, H and S (E^H^S)
NP_*** present in E, S and C (E^S^C)
It is very difficult for me to deduce the following:
NP_**** entries specific to E
NP_**** entries specific to H
NP_**** entries specific to S
I am not sure about your notation, but in R terms, if you want the items
specific to E, H, or S, how about something like:
specific to E -- setdiff(e, union(c,union(h,s)))
specific to H -- setdiff(h, union(c,union(e,s)))
specific to S -- setdiff(s, union(c,union(e,h)))
Dan Nordlund
[[alternative HTML version deleted]]
______________________________________________
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html