Hi,
can anybody help me with the BiomartGeneRegionTrack from the Gviz package of Bioconductor? My problem is the subset function. After creating a BiomartGeneRegionTrack for the chromosome 1 of hg19 from 179950000 to 180000000 this track contains 103 items.
In the following the code is depicted:
start <- 179950000;
end <- 180000000
track <- BiomartGeneRegionTrack(genome = 'hg19', chromosome = 'chr1', start = start, end = end)
length(track)
After creating a subset of this track from the same region there are only 51 items left.
subTrack <- subset(newTrack, from=start, to=end)
length(subTrack)
You can also plot these tracks with the command:
plotTracks(list(subTrack,newTrack), showId=TRUE)
Does anybody now how I can adjust this subset function that the result for this example are the same items then before subsetting?
Kind regards, Andreas
My sessionInfo()
shows the following information:
R version 3.0.1 (2013-05-16)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=de_AT.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=de_AT.UTF-8 LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=de_AT.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] grDevices datasets utils parallel stats graphics grid
[8] methods base
other attached packages:
[1] TxDb.Hsapiens.UCSC.hg19.knownGene_2.9.0
[2] org.Hs.eg.db_2.9.0
[3] RSQLite_0.11.3
[4] DBI_0.2-7
[5] BSgenome.Hsapiens.UCSC.hg19_1.3.19
[6] BSgenome_1.28.0
[7] Biostrings_2.28.0
[8] GenomicFeatures_1.12.1
[9] AnnotationDbi_1.22.5
[10] Biobase_2.20.0
[11] GenomicRanges_1.12.2
[12] IRanges_1.18.1
[13] BiocGenerics_0.6.0
[14] Gviz_1.4.1
loaded via a namespace (and not attached):
[1] biomaRt_2.16.0 biovizBase_1.8.0 bitops_1.0-5 cluster_1.14.4
[5] colorspace_1.2-2 dichromat_2.0-0 Hmisc_3.10-1.1 labeling_0.1
[9] lattice_0.20-15 munsell_0.4 plyr_1.8 RColorBrewer_1.0-5
[13] RCurl_1.95-4.1 Rsamtools_1.12.3 rtracklayer_1.20.2 scales_0.2.3
[17] stats4_3.0.1 stringr_0.6.2 tools_3.0.1 XML_3.96-1.1
[21] zlibbioc_1.6.0