Gviz plotting gene region
1
0
Entering edit mode
swebb1 • 0
@swebb1-11141
Last seen 5.7 years ago

I am using Gviz to plot a region given by a gene name, for instance:

biomTrack <- BiomartGeneRegionTrack(genome ="sacCer3",name="ENSEMBL",gene="YBL087C")
plotTracks(list(btrack))

This works fine. I also want to add another annotation track:

btrack<-AnnotationTrack(start = c(100,200,300),width=100,chromosome = "chrII",name="test",showFeatureId=T)
plotTracks(list(biomTrack,btrack))

Now the region displayed appears to be all of chrII. I know I can use the from/to parameters in plotTracks but I want to select regions by gene name. Is this possible?

 
gviz • 874 views
ADD COMMENT
0
Entering edit mode
@florianhahnenovartiscom-3784
Last seen 5.6 years ago
Switzerland

This is not possible. Unless explicitly provided via the from and to arguments, Gviz will derive the plotting region from the extremes of all provided tracks. So in your case it will plot everything from 100 to the rightmost exon of gene yBL08. You could of course use the ranges in biomTrack as limits:

plotTracks(list(biomTrack,btrack), from=min(biomTrack), to=max(biomTrack))

However I am not quite sure how that is supposed to make any sense, since the ranges you provide in btrack do not overlap with YBL08....

Florian

ADD COMMENT

Login before adding your answer.

Traffic: 967 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6