Does this do what you want:


x <- "Exon    start           end
5'UTR   22540060        22540121
1       22540122        22540140
2       22540303        22540493
3       22541552        22541565
4       22542373        22542519
5       22544265        22544432
3'UTR   22544433        22544856"

y <- read.table(textConnection(x), quote='', header=TRUE)

# create the plot area
plot(range(y$start, y$end), c(1,7), type='n', xaxt='n', yaxt='n',
    ylab='', xlab='')

# create the colored rectangles
rect(y$start, 1, y$end, 2, col=c('red', rep('green', nrow(y)-2), 'red'))

# put labels at midpoints of rectangles
axis(1, at=(y$start + y$end)/2, labels=as.character(y$Exon))



On 7/10/06, kannaiah@bsd.uchicago.edu <kannaiah@bsd.uchicago.edu> wrote:
>
> Hello,
>
> I am new to R and still feeling my way thru it.
>
> I am trying to plot the values from this file below on the X-axis of a
> plot. I
> have attached the graph to the email...the one i am trying to recreate.
>
> Exon    start           end
> 5'UTR   22540060        22540121
> 1       22540122        22540140
> 2       22540303        22540493
> 3       22541552        22541565
> 4       22542373        22542519
> 5       22544265        22544432
> 3'UTR   22544433        22544856
>
>
> I would like to create small rectangles on the x-axis as colored boxes
> from
> start position to end position of each exon...with the label showing 1, 2
> etc
> under each respective exon.
>
> The 5' and 3' UTR would be colored different from the exons 1-5.
>
>
> Any suggestions and ideas would greatly appreciated.
>
> Thank you
> Kiran
>
>
> -------------------------------------------------
> This email is intended only for the use of the individual or entity to
> which
> it is addressed and may contain information that is privileged and
> confidential.  If the reader of this email message is not the intended
> recipient, you are hereby notified that any dissemination, distribution,
> or
> copying of this communication is prohibited.  If you have received this
> email
> in error, please notify the sender and destroy/delete all copies of the
> transmittal.  Thank you.
> -------------------------------------------------
>
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>
>


-- 
Jim Holtman
Cincinnati, OH
+1 513 646 9390 (Cell)
+1 513 247 0281 (Home)

What is the problem you are trying to solve?

	[[alternative HTML version deleted]]

