Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Parsing issue in pathview package

    I got a parsing error when using pathview package to visualize pathways in my RNAseq data. This error still showed up even I tried to run the example (quick start with demo data) in the vignette. I found that XML was loaded successfully as shown in sessionInfo. Could someone help me to see what's wrong? Thank you very much!

    > library(pathview)
    Loading required package: org.Hs.eg.db
    Loading required package: AnnotationDbi
    Loading required package: stats4
    Loading required package: BiocGenerics
    Loading required package: parallel

    Attaching package: ‘BiocGenerics’

    The following objects are masked from ‘packagearallel’:

    clusterApply, clusterApplyLB, clusterCall, clusterEvalQ, clusterExport,
    clusterMap, parApply, parCapply, parLapply, parLapplyLB, parRapply,
    parSapply, parSapplyLB

    The following objects are masked from ‘package:stats’:

    IQR, mad, sd, var, xtabs

    The following objects are masked from ‘package:base’:

    anyDuplicated, append, as.data.frame, cbind, colMeans, colnames,
    colSums, do.call, duplicated, eval, evalq, Filter, Find, get, grep,
    grepl, intersect, is.unsorted, lapply, lengths, Map, mapply, match,
    mget, order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank,
    rbind, Reduce, rowMeans, rownames, rowSums, sapply, setdiff, sort,
    table, tapply, union, unique, unsplit, which, which.max, which.min

    Loading required package: Biobase
    Welcome to Bioconductor

    Vignettes contain introductory material; view with 'browseVignettes()'.
    To cite Bioconductor, see 'citation("Biobase")', and for packages
    'citation("pkgname")'.

    Loading required package: IRanges
    Loading required package: S4Vectors

    Attaching package: ‘S4Vectors’

    The following object is masked from ‘package:base’:

    expand.grid


    ##############################################################################
    Pathview is an open source software package distributed under GNU General
    Public License version 3 (GPLv3). Details of GPLv3 is available at
    http://www.gnu.org/licenses/gpl-3.0.html. Particullary, users are required to
    formally cite the original Pathview paper (not just mention it) in publications
    or products. For details, do citation("pathview") within R.

    The pathview downloads and uses KEGG data. Non-academic uses may require a KEGG
    license agreement (details at http://www.kegg.jp/kegg/legal.html).
    ##############################################################################
    > data(gse16873.d)
    > pv.out <- pathview(gene.data = gse16873.d[, 1], pathway.id = "04110", species = "hsa", out.suffix = "gse16873")
    Info: Downloading xml files for hsa04110, 1/1 pathways..
    Info: Downloading png files for hsa04110, 1/1 pathways..
    Warning: Parsing ./hsa04110.xml file failed, please check the file!
    > sessionInfo()
    R version 3.4.1 (2017-06-30)
    Platform: x86_64-w64-mingw32/x64 (64-bit)
    Running under: Windows 10 x64 (build 15063)

    Matrix products: default

    locale:
    [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252
    [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C
    [5] LC_TIME=English_United States.1252

    attached base packages:
    [1] parallel stats4 stats graphics grDevices utils datasets methods
    [9] base

    other attached packages:
    [1] pathview_1.16.1 org.Hs.eg.db_3.4.1 AnnotationDbi_1.38.1
    [4] IRanges_2.10.2 S4Vectors_0.14.3 Biobase_2.36.2
    [7] BiocGenerics_0.22.0

    loaded via a namespace (and not attached):
    [1] graph_1.54.0 Rcpp_0.12.12 KEGGgraph_1.38.0 XVector_0.16.0
    [5] zlibbioc_1.22.0 bit_1.1-12 R6_2.2.2 rlang_0.1.1
    [9] blob_1.1.0 httr_1.2.1 tools_3.4.1 grid_3.4.1
    [13] png_0.1-7 DBI_0.7 bit64_0.9-7 digest_0.6.12
    [17] tibble_1.3.3 Rgraphviz_2.20.0 KEGGREST_1.16.0 memoise_1.1.0
    [21] RSQLite_2.0 compiler_3.4.1 Biostrings_2.44.1 XML_3.98-1.9
    [25] pkgconfig_2.0.1

  • #2
    I happen to keep another version of R (3.3.1) on my PC. I tried the same code using R_3.3.1 and it worked. I'm not sure whether it is a version issue. But at least it works. Just let anyone who's interested in known.

    > library("pathview")
    > data(gse16873.d)
    > pv.out <- pathview(gene.data = gse16873.d[, 1], pathway.id = "04110", species = "hsa", out.suffix = "gse16873")
    Info: Downloading xml files for hsa04110, 1/1 pathways..
    Info: Downloading png files for hsa04110, 1/1 pathways..
    'select()' returned 1:1 mapping between keys and columns
    Info: Working in directory C:/Users/mahon/Desktop
    Info: Writing image file hsa04110.gse16873.png
    > list.files(pattern = "hsa04110*")
    [1] "hsa04110.gse16873.png" "hsa04110.png" "hsa04110.xml"

    Comment


    • #3
      Thank you so much for reporting your results! I was having the same issue and it was driving me nuts.

      I used Shiny's debugging breakpoints to drill into the pathview function, and traced the error to what seems to be a namespace problem between `pathview` and `KEGGgraph`. If you print the error message passed to tryCatch, you get:

      ```
      Error in parsePathwayInfo(r): could not find function "parsePathwayInfo"
      ```

      If you simply redefine `parsePathwayInfo()` in the global namespace, you get the same error with `xmlAttrs()`.

      Fortunately pathview functions normally in R 3.3.x, but hopefully the package developer can fix this for R 3.4!

      Comment


      • #4
        Weijun is very nice to help to figure out the issue and give a solution posted in Bioconductor website (see below). Thank you, Weijun!

        Comment


        • #5
          Dear all,
          Thank you for reporting this issue. Indeed, I was able to replicate the same problem with the latest R/BioC versions.
          I just posted the cause and solution for this problem on bioc support forum:


          you may always use pathview web server as an backup in case something like this occurred again:

          Pathview Web server provides a intuitive GUI and API to Pathview, a full pathway analysis workflow supporting multiple omics data and integrated analysis, and much more.

          please do not hesitate to contact us for any future issues. Thank you for your patience and support to pathview.

          Comment


          • #6
            pathview package is running error-free now with the latest R/BioC and KEGGgraph:

            Pathview is a tool set for pathway based data integration and visualization. It maps and renders a wide variety of biological data on relevant pathway graphs. All users need is to supply their data and specify the target pathway. Pathview automatically downloads the pathway graph data, parses the data file, maps user data to the pathway, and render pathway graph with the mapped data. In addition, Pathview also seamlessly integrates with pathway and gene set (enrichment) analysis tools for large-scale and fully automated analysis.


            Thanks to everyone!

            Comment


            • #7
              Thank you so much!

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Strategies for Sequencing Challenging Samples
                by seqadmin


                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                03-22-2024, 06:39 AM
              • seqadmin
                Techniques and Challenges in Conservation Genomics
                by seqadmin



                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                Avian Conservation
                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                03-08-2024, 10:41 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, Yesterday, 06:37 PM
              0 responses
              11 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              10 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              51 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              68 views
              0 likes
              Last Post seqadmin  
              Working...
              X