Re: [R] read.frame

From: <DJNordlund_at_aol.com>
Date: Thu 08 Jul 2004 - 15:35:07 EST

 

In a message dated 7/7/2004 9:42:03 PM Pacific Daylight Time, biocperi@yahoo.com 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]]



R-help@stat.math.ethz.ch 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 Received on Thu Jul 08 15:39:30 2004

This archive was generated by hypermail 2.1.8 : Wed 03 Nov 2004 - 22:54:46 EST