[R] finding most highly transcribed genes - ranking, sorting and subsets?

From: alison waller <alison.waller_at_utoronto.ca>
Date: Thu, 6 Dec 2007 14:03:28 -0500


I am not only interested in finding out which genes are the most highly upor  down-regulated (which I have done using the linear models and Bayesian statistics in Limma), but I also want to know which genes are consistently highly transcribed (ie. they have a high intensity in the channel of interest eg. Cy5 or Cy3 across the set of experiments). I might have missed a straight forward way to do this, or a valuable function, but I've been using my own methods and going around in circles.  

So far I've normalized within and between arrays, then returned the RG values using RG<-RG.MA, then I ranked each R and G values for each array as below.



rankRG$R[,2]<-rank(rankRG$R[,2]) .. and so on for 6 columns(ie. arrays, as well as the G's)  

then I thought I could pull out a subset of rankRG using something like;



However, this just returned me a matrix with one row of $R (the ranks were <500 for columns 1,2, and 5 and greater than 500 for 3,4,and 6). However, I can't believe that there is only one gene that is in the top 500 for R intensitiy among those three arrays.  

Am I doing something wrong? Can someone think of a better way of doing this?  



Alison S. Waller M.A.Sc.
Doctoral Candidate
416-978-4222 (lab)
Department of Chemical Engineering
Wallberg Building
200 College st.
Toronto, ON
M5S 3E5     

        [[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. Received on Thu 06 Dec 2007 - 20:11:24 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 06 Dec 2007 - 20:30:17 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.