[R] Adding markers and text for some data points after drawing a plot

David Winsemius dwinsemius at comcast.net
Tue Feb 26 23:32:36 CET 2013


On Feb 26, 2013, at 2:28 PM, Debs Majumdar wrote:

> Hi David, 
> 
> Thanks for your suggestions.
> 
> I am not sure how to plot the points by using gene names, but I will give it a try.

The x and y arguments would be the same as were used in the text() call.

(You are _still_ sending formatted text.)
http://help.yahoo.com/kb/index?locale=en_US&y=PROD_ACCT&page=content&id=SLN3514


-- 
David.
> 
> Thanks,
> 
> Debs
> 
> From: David Winsemius <dwinsemius at comcast.net>
> To: Debs Majumdar <debs_stata at yahoo.com> 
> Cc: "r-help at r-project.org" <r-help at r-project.org> 
> Sent: Tuesday, February 26, 2013 1:54 PM
> Subject: Re: [R] Adding markers and text for some data points after drawing a plot
> 
> 
> On Feb 26, 2013, at 1:04 PM, Debs Majumdar wrote:
> 
> > 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)
> 
> Using 'attach' is a bad idea. In particular we have trouble keeping straight what is a column name or a first class object in your workspace.
> 
> > 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
> 
> Use points()
> 
> > 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.
> 
> 
> text(x= d[ which( d[['gene']]=="AADACL4"), 2], 
>     y= kblength.ecdf( d[ which( dat[['gene']]=="AADACL4"), 2]), 
>     "AADACL4")
> > 
> > I am fighting with this but wasn't able to get a solution.
> 
> > 
> > Any help is appreciated.
> > 
> > Thanks,
> >  Debs
> > 
> >     [[alternative HTML version deleted]]
> 
> Please learn to post in plain text. And please, please do read the fine Posting Guide.
> 
> -- 
> David Winsemius
> Alameda, CA, USA
> 
> 
> 

David Winsemius
Alameda, CA, USA



More information about the R-help mailing list