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.
Thanks,
Debs
________________________________
From: David Winsemius <dwinsem...@comcast.net>
Cc: "r-help@r-project.org" <r-help@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
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
