Entering edit mode
Hello,
I get an HTTP connection error when using Gviz::IdeogramTrack():
library(Gviz)
IdeogramTrack(genome="hg19", chromosome="chr7")
Error: failed to load HTTP resource
traceback()
18: stop(e)
17: (function (msg, code, domain, line, col, level, filename, class = "XMLError")
{
e = makeXMLError(msg, code, domain, line, col, level, filename,
class)
dom = names(e$domain)
class(e) = c(names(e$code), sprintf("%s_Error", gsub("_FROM_",
"_", dom)), class(e))
if (e$code == xmlParserErrors["XML_IO_LOAD_ERROR"])
stop(e)
})("failed to load HTTP resource\n", 1549L, 8L, 0L, 0L, 1L, character(0))
16: .Call("RS_XML_ParseTree", as.character(file), handlers, as.logical(ignoreBlanks),
as.logical(replaceEntities), as.logical(asText), as.logical(trim),
as.logical(validate), as.logical(getDTD), as.logical(isURL),
as.logical(addAttributeNamespaces), as.logical(useInternalNodes),
as.logical(isHTML), as.logical(isSchema), as.logical(fullNamespaceInfo),
as.character(encoding), as.logical(useDotNames), xinclude,
error, addFinalizer, as.integer(options), as.logical(parentFirst),
PACKAGE = "XML")
15: htmlParse(doc, encoding = encoding)
14: .local(doc, header, colClasses, skip.rows, trim, elFun, as.data.frame,
which, ...)
13: readHTMLTable(url)
12: readHTMLTable(url)
11: ucscGenomes()
10: eval(expression, envir = callEnv)
9: eval(expression, envir = callEnv)
8: eval(expression, envir = callEnv)
7: .doCache(availToken, expression(ucscGenomes()), env)
6: .cacheGenomes(genome = genome)
5: .local(.Object, ...)
4: initialize(value, ...)
3: initialize(value, ...)
2: new("IdeogramTrack", chromosome = chromosome, genome = genome,
name = name, bands = bands, ...)
1: IdeogramTrack(genome = "hg19", chromosome = "chr7")
The connection between rtracklayer and UCSC is working:
library(rtracklayer)
session <- browseSession()
track(session, "knownGene")
UCSC track 'knownGene'
UCSCData object with 197782 ranges and 5 metadata columns:
seqnames ranges strand | name score itemRgb thick blocks
<Rle> <IRanges> <Rle> | <character> <numeric> <character> <IRanges> <IRangesList>
[1] chr1 [17369, 17436] - | uc031tla.1 0 <NA> [17369, 17368] [1, 68]
[2] chr1 [29554, 31097] + | uc057aty.1 0 <NA> [29554, 29553] [ 1, 486] [1011, 1114] [1423, 1544]
[3] chr1 [30267, 31109] + | uc057atz.1 0 <NA> [30267, 30266] [ 1, 401] [710, 843]
[4] chr1 [30366, 30503] + | uc031tlb.1 0 <NA> [30366, 30365] [1, 138]
[5] chr1 [34554, 36081] - | uc001aak.4 0 <NA> [34554, 34553] [ 1, 621] [ 724, 928] [1168, 1528]
... ... ... ... . ... ... ... ... ...
[197778] chr22_KI270733v1_random [179370, 179418] + | uc062fri.1 0 <NA> [179370, 179369] [1, 49]
[197779] chr22_KI270734v1_random [ 72411, 74814] + | uc062frj.1 0 <NA> [ 72838, 74124] [1, 2404]
[197780] chr22_KI270734v1_random [131494, 137392] + | uc032quk.2 0 <NA> [131646, 136994] [ 1, 262] [ 343, 503] [3950, 4050] ...
[197781] chr22_KI270734v1_random [138082, 161750] - | uc032qun.2 0 <NA> [138480, 156446] [ 1, 586] [ 662, 750] [4113, 4211] ...
[197782] chr22_KI270734v1_random [138082, 161852] - | uc032qum.2 0 <NA> [138480, 161586] [ 1, 586] [ 662, 750] [4113, 4211] ...
-------
seqinfo: 455 sequences from hg38 genome
I updated both rtracklayer and Gviz today. Any pointers on how to fix this would be much appreciated.
Thanks!
Enrico
sessionInfo() R version 3.4.2 (2017-09-28) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 7 x64 (build 7601) Service Pack 1 Matrix products: default locale: [1] LC_COLLATE=English_United Kingdom.1252 LC_CTYPE=English_United Kingdom.1252 LC_MONETARY=English_United Kingdom.1252 LC_NUMERIC=C LC_TIME=English_United Kingdom.1252 attached base packages: [1] grid parallel stats4 stats graphics grDevices utils datasets methods base other attached packages: [1] rtracklayer_1.38.3 Gviz_1.22.2 GenomicRanges_1.30.0 GenomeInfoDb_1.14.0 IRanges_2.12.0 S4Vectors_0.16.0 BiocGenerics_0.24.0 loaded via a namespace (and not attached): [1] Biobase_2.38.0 httr_1.3.1 RMySQL_0.10.13 AnnotationHub_2.10.1 bit64_0.9-7 splines_3.4.2 shiny_1.0.5 [8] Formula_1.2-2 assertthat_0.2.0 interactiveDisplayBase_1.16.0 latticeExtra_0.6-28 blob_1.1.0 BSgenome_1.46.0 GenomeInfoDbData_1.0.0 [15] Rsamtools_1.30.0 yaml_2.1.16 progress_1.1.2 pillar_1.1.0 RSQLite_2.0 backports_1.1.2 lattice_0.20-35 [22] biovizBase_1.26.0 digest_0.6.15 RColorBrewer_1.1-2 XVector_0.18.0 checkmate_1.8.5 colorspace_1.3-2 httpuv_1.3.5 [29] htmltools_0.3.6 Matrix_1.2-12 plyr_1.8.4 XML_3.98-1.9 biomaRt_2.34.2 zlibbioc_1.24.0 xtable_1.8-2 [36] scales_0.5.0 BiocParallel_1.12.0 htmlTable_1.11.2 tibble_1.4.2 AnnotationFilter_1.2.0 ggplot2_2.2.1 SummarizedExperiment_1.8.1 [43] GenomicFeatures_1.30.1 nnet_7.3-12 lazyeval_0.2.1 mime_0.5 survival_2.41-3 magrittr_1.5 memoise_1.1.0 [50] foreign_0.8-69 BiocInstaller_1.28.0 tools_3.4.2 data.table_1.10.4-3 prettyunits_1.0.2 matrixStats_0.53.0 stringr_1.2.0 [57] munsell_0.4.3 cluster_2.0.6 DelayedArray_0.4.1 ensembldb_2.2.0 AnnotationDbi_1.40.0 Biostrings_2.46.0 compiler_3.4.2 [64] rlang_0.1.6 RCurl_1.95-4.10 dichromat_2.0-0 rstudioapi_0.7 VariantAnnotation_1.24.5 htmlwidgets_1.0 bitops_1.0-6 [71] base64enc_0.1-3 gtable_0.2.0 curl_3.1 DBI_0.7 R6_2.2.2 GenomicAlignments_1.14.1 gridExtra_2.3 [78] knitr_1.19 bit_1.1-12 Hmisc_4.1-1 ProtGenerics_1.10.0 stringi_1.1.6 Rcpp_0.12.15 rpart_4.1-11 [85] acepack_1.4.1

This might be related to gviz ideogramTrack connection error though the error message appears to be different.
Hopefully just a transient issue. This is a simple download of a web page, via the XML package and its embedded nanohttp library (part of libxml2).
will you explain more about it
I too face the similar issue
itrack <- IdeogramTrack(genome = gen, chromosome = chr) Error: failed to load HTTP resource