Question: Question about BiomartGeneRegionTrack
gravatar for
5.3 years ago by
What you are trying to do is not that straight forward. What do you want to do with overlapping exons from different transcript models? Do you just want to merge them? And how treat UTRs that are overlapping with coding exons in other transcripts, as it is the case for your THADA gene? In principle the tracks in the Gviz package are just GRanges objects, so you could modify them using existing tools. You may want to treat the separate genes individually, so splitting the GRanges object into a GRangesList should help. I don't think that there is anything useful that could be done with the UTR information, so I'd skip it completely and just treat all exons as coding exons. The following untested code might do what you want to archive: gen="hg19" chr="2" start=43625705 end= 43826133 biomTrack <- BiomartGeneRegionTrack(genome = gen, chromosome = chr, start = start, end = end, name = "ENSEMBL") plotTracks(biomTrack, showId=TRUE) feature(biomTrack) <- "protein_coding" r <- split(ranges(biomTrack), gene(biomTrack)) rNew <- endoapply(r, function(x){ rx <- reduce(x, with.mapping=TRUE) mcols(rx) <- mcols(x)[sapply(rx$mapping, head, 1),] rx$transcript <- rx$gene rx }) ranges(biomTrack) <- unlist(rNew) plotTracks(biomTrack, showId=TRUE) Florian From: Tiphaine Martin <<>> Date: Friday, February 21, 2014 5:19 PM To: Florian Hahne <<>> Cc: "<>" <<>> Subject: Question about BiomartGeneRegionTrack end = end, name = "ENSEMBL",stacking="squish") [[alternative HTML version deleted]]
gviz • 734 views
ADD COMMENTlink written 5.3 years ago by florian.hahne@novartis.com1.6k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 261 users visited in the last hour