Another plot question I am afraid. Thanks for all the help with the previous one. I have the expression levels of a number of exons in a gene. I would look like to see how the expression levels vary across the gene. For each of the exons I have a start position and a stop position. It is easy enough to plot the expression level verses the start/stop position: > plot(exonAnnot[rownames(temp),names(exonAnnot)=="start"],temp[,16], main=names(temp)[16]) where temp is the expression matrix for the relevant exons. What I would like to do is for each exon draw a bar from the start position to the stop position at the level of expression. Anyway to do this? Many thanks Dan -- ************************************************************** Daniel Brewer, Ph.D. Institute of Cancer Research Molecular Carcinogenesis Email: daniel.brewer at icr.ac.uk ************************************************************** The Institute of Cancer Research: Royal Cancer Hospital, a charitable Company Limited by Guarantee, Registered in England under Company No. 534147 with its Registered Office at 123 Old Brompton Road, London SW7 3RP. This e-mail message is confidential and for use by the addre...{{dropped}}

On Monday 19 March 2007 11:18, Daniel Brewer wrote: > Another plot question I am afraid. Thanks for all the help with the > previous one. > I have the expression levels of a number of exons in a gene. I would > look like to see how the expression levels vary across the gene. For > each of the exons I have a start position and a stop position. > > It is easy enough to plot the expression level verses the start/stop > > position: > > plot(exonAnnot[rownames(temp),names(exonAnnot)=="start"],temp[,16], > > main=names(temp)[16]) > > where temp is the expression matrix for the relevant exons. > > What I would like to do is for each exon draw a bar from the start > position to the stop position at the level of expression. > > Anyway to do this? You might want to check the rect() function. Sean