2013-04-01 35 views
1

是否有经验的人使用bioconductor程序包:Gviz知道如何通过DataTrack直接添加AnnotationTrackGviz中的注释轨迹

例如,在ggplot2我可以使用+ geom_text添加到prexitsting情节,但我一直没能找到一个类似的功能为Gviz

谢谢!

+0

您是否阅读了小插曲,例如'vignette(“Gviz”)'或http://bioconductor.org/packages/release/bioc/html/Gviz.html,特别是第2节快速浏览?我想你想使用'plotTracks()',提供一个轨道列表作为第一个参数。 –

+0

是的!我一直在梳理文档...'plotTracks()'会在同一个plot/pdf中同时绘制'AnnotationTrack'和'DataTrack' - 但是一个接一个,我希望使用Annotation Track功能添加一个“DataTrack”数据直接位于特定染色体区域的几个关键元素(例如,着丝粒落在与数据中观察到的模式相关的位置等) – mfk534

+0

我意识到可以在“IdeogramTrack”中观察着丝粒'但是我再次感兴趣的是覆盖两条曲目的可能性,以便更直观地表现出特定的关系。 – mfk534

回答

0

虽然这不完全是你想要的,但一种可能的解决方案是添加一个覆盖您感兴趣的区域的HighlightTrack。虽然这不会在您的DataTrack上特别标注/添加关键元素,但它有助于突出显示不同DataTracks之间的对齐。

实施例:

library(Gviz) 
library(GenomicRanges) 
data(geneModels) 
data(cpgIslands) 
gen <- genome(cpgIslands) 
chr <- 1 
start <- 120005434 
end <- 129695434 
itrack <- IdeogramTrack(genome = gen, chromosome = chr) 
gtrack <- GenomeAxisTrack() 
grtrack <- GeneRegionTrack(geneModels, genome = gen, chromosome = chr, name = "foo") 
htrack <- HighlightTrack(trackList = list(gtrack, grtrack), start = 121535434, end = 124535434, chromosome = chr) 
plotTracks(list(itrack, htrack), from = start, to = end) 

图形输出

enter image description here

虽然grtrack是空的,它说明了如何将HighlightTrack将跨越指定DataTracks(在这种情况下,grtrackgtrack )。

有关详细信息,请参阅GViz documentation