[BioC] Gviz UcscTrack error
Anna Sawicka
anna.sawicka at univie.ac.at
Thu Jan 9 10:43:50 CET 2014
Dear Florian,
thank you for the response. Everything works today. This must have been a server problem. And thanks a lot for you suggestion concerning RefSeq gene models!
Best wishes,
Ania
> Hi Ania,
> I can't seem to reproduce your problem here:
>> library(Gviz)
> Loading required package: grid
>> genetrack=UcscTrack(track="RefSeq Genes", table="refGene",
>
> + trackType="GeneRegionTrack",chromosome="chr5",
> genome="mm9",
> + rstart="exonStarts", rends="exonEnds", gene="name",
> symbol="name2",
> + transcript="name", strand="strand", name="RefSeq
> Genes", feature="name2",
> + showId=T, from=30334169, to=30349988, fill = "green")
>> sessionInfo()
> R version 3.0.2 Patched (2013-10-27 r64116)
> Platform: i386-apple-darwin12.5.0/i386 (32-bit)
>
> locale:
> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
>
> attached base packages:
> [1] parallel grid stats graphics grDevices utils datasets
> [8] methods base
>
> other attached packages:
> [1] GenomicFeatures_1.14.2 AnnotationDbi_1.24.0 Biobase_2.22.0
> [4] GenomicRanges_1.14.4 XVector_0.2.0 IRanges_1.20.6
> [7] BiocGenerics_0.8.0 Gviz_1.6.0 BiocInstaller_1.12.0
>
> loaded via a namespace (and not attached):
> [1] biomaRt_2.18.0 Biostrings_2.30.1 biovizBase_1.10.7
> [4] bitops_1.0-6 BSgenome_1.30.0 cluster_1.14.4
> [7] colorspace_1.2-4 DBI_0.2-7 dichromat_2.0-0
> [10] Formula_1.1-1 Hmisc_3.13-0 labeling_0.2
> [13] lattice_0.20-24 latticeExtra_0.6-26 munsell_0.4.2
> [16] plyr_1.8 RColorBrewer_1.0-5 RCurl_1.95-4.1
> [19] Rsamtools_1.14.2 RSQLite_0.11.4 rtracklayer_1.22.0
> [22] scales_0.2.3 splines_3.0.2 stats4_3.0.2
> [25] stringr_0.6.2 survival_2.37-4 tools_3.0.2
> [28] XML_3.98-1.1 zlibbioc_1.8.0
>
> Maybe the UCSC server had some issue on the day you tried this? Could you
> please check again today whether the problem has gone away? Both Gviz and
> rtracklayer versions are the same on my machine and yours, so package
> versions should not be an issue. The only difference I see in the
> sessionInfo is a slightly newer R version in my end (3.0.2 Patched)
>
> If all that you are interested in are RefSeq gene models from UCSC you
> might as well consider building a GeneRegionTrack directly from a
> TranscriptDB object. This should be much more efficient and you will also
> get UTR information for free. Please see the package vignette for details.
>
> Florian
>
>
>
>
>
>
> On 1/7/14 10:44 PM, "Anna Sawicka [guest]" <guest at bioconductor.org> wrote:
>
>> genetrack=UcscTrack(track="RefSeq Genes", table="refGene",
>> trackType="GeneRegionTrack",chromosome="chr5",
>> genome="mm9",
>> rstart="exonStarts", rends="exonEnds", gene="name",
>> symbol="name2",
>> transcript="name", strand="strand", name="RefSeq
>> Genes", feature="name2",
>> showId=T, from=30334169, to=30349988, fill = "green")
>
More information about the Bioconductor
mailing list