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 <dwinsem...@comcast.net> > To: Debs Majumdar <debs_st...@yahoo.com> > 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 > > > David Winsemius Alameda, CA, USA ______________________________________________ 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.
