I would like to run the "BiomartGeneRegionTrack" in package Gviz to get the annotation file. But it seems I can not use the range function of it. Every time I get a start and end to "BiomartGeneRegionTrack" it will report an error. Moreover, "plotTracks" also shows nothing. The code is below:
> afrom <- 2960000 > ato <- 3160000 > bmt <- BiomartGeneRegionTrack(genome="hg19", chromosome="chr12", start=afrom, end=ato) Warning message: In .Seqinfo.mergexy(x, y) : The 2 combined objects have no sequence levels in common. (Use suppressWarnings() to suppress this warning.) >plotTracks(bmt, from=afrom, to=ato, chromosome="chr12")
Could anyone give me a tip about how to solve this problem? Thank you in advance.