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
                Recent Developments in Metagenomics
                by seqadmin





                Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
                09-23-2024, 06:35 AM
              • seqadmin
                Understanding Genetic Influence on Infectious Disease
                by seqadmin




                During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

                Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
                09-09-2024, 10:59 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 10-02-2024, 04:51 AM
              0 responses
              8 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 10-01-2024, 07:10 AM
              0 responses
              14 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 09-30-2024, 08:33 AM
              0 responses
              18 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 09-26-2024, 12:57 PM
              0 responses
              16 views
              0 likes
              Last Post seqadmin  
              Working...
              X