Hi All,

I have a data set of around 17,000 gene names and their lengths. E.g.

          gene      kblength    
 A3GALT2       14.333    
AADACL3       12.609    
AADACL4       22.532    
     ABCA4     128.312    
    ABCB10      42.114    
      ABCD3    100.287    
............
.........

and I was able to draw a reverse cumulative frequency plot using the following code:
d <- read.table("gene_length_data.txt", header=T)
attach(d)
kblength.ecdf <-ecdf(kblength)
r <-range(kblength)
curve(1-kblength.ecdf(x), from=r[1], to=r[2], col="red", xlim=r)
detach(d)

Here's the link to the output: http://www.tiikoni.com/tis/view/?id=5c41810

I want to put markers and the corresponding gene names in the graph for a small list of 6 genes. i.e. say a hollow circle on the curve followed by the corresponding gene name.

I am fighting with this but wasn't able to get a solution.

Any help is appreciated.

Thanks,
 Debs

	[[alternative HTML version deleted]]

