Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • emolinari
    Member
    • May 2013
    • 47

    #1

    GoSeq -unable to calculate pwf!

    hello All!
    I am desperately trying to run this GO analysis on my RNAseq data(human samples).
    I got stuck in creating a pwf file: I have -I guess- successfully created a DE matrix, with ENS id and 0/1 as identifier of diff-non diff expressed genes.
    here is the error message I encounter
    genes
    0 1
    19004 6715
    > pwf <- nullp(genes, "hg19", "ensGene")
    Loading hg19 length data...
    Error in qr.R(qrx) :
    could not find symbol "..." in environment of the generic function
    I noticed that the ENS ids in my matrix are not into double quotes...I am not sure if it really has an effect...I have just followed the manual step by step, but I am not able to figure out why is not going through (and I am newbie in all this).
    Anybody can help?!!

    Manu
    > sessionInfo()
    R version 2.15.2 (2012-10-26)
    Platform: x86_64-apple-darwin9.8.0/x86_64 (64-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] stats graphics grDevices utils datasets methods base

    other attached packages:
    [1] goseq_1.10.0 geneLenDataBase_0.99.10 BiasedUrn_1.05 edgeR_3.0.8 limma_3.14.4 BiocInstaller_1.8.3

    loaded via a namespace (and not attached):
    [1] AnnotationDbi_1.20.7 Biobase_2.18.0 BiocGenerics_0.4.0 biomaRt_2.14.0 Biostrings_2.26.3 bitops_1.0-4.2 BSgenome_1.26.1
    [8] DBI_0.2-5 GenomicFeatures_1.10.2 GenomicRanges_1.10.7 grid_2.15.2 IRanges_1.16.6 lattice_0.20-15 Matrix_1.0-12
    [15] mgcv_1.7-22 nlme_3.1-108 parallel_2.15.2 RCurl_1.95-4.1 Rsamtools_1.10.2 RSQLite_0.11.2 rtracklayer_1.18.2
    [22] stats4_2.15.2 tools_2.15.2 XML_3.96-1.1 zlibbioc_1.4.0
  • Michael_RNAseq
    Junior Member
    • Aug 2012
    • 3

    #2
    Hi,

    The description is a little vague but I think I can help troubleshoot the problem.

    1. Are you sure you are working with 'ensGene' symbols? Could your genes be better classified as something else?

    2. Did you manipulate your ensGene ID's at all? Are you supplying two different types of ids in your complete genes list and DE genes list?

    3. Make sure you are running the correct R version 3.0 and bioconductor version 2.12. I was having issues with older versions and this helped.

    Michael
    Last edited by Michael_RNAseq; 05-30-2013, 09:29 AM.

    Comment

    • emolinari
      Member
      • May 2013
      • 47

      #3
      Hi Michael!!!

      Thanks for the hints... Actually I have already updated my R version and now everything works beautifully!!!

      So if it can be helpful for anyone else...!

      Thanks guys!
      Manu

      Comment

      Latest Articles

      Collapse

      • SEQadmin2
        16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
        by SEQadmin2




        NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...
        10-30-2026, 06:46 AM
      • SEQadmin2
        Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
        by SEQadmin2



        Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

        Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
        ...
        10-14-2026, 09:55 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 10-28-2026, 09:57 AM
      0 responses
      33 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-28-2026, 09:51 AM
      0 responses
      21 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-22-2026, 11:13 AM
      0 responses
      72 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-22-2026, 10:48 AM
      0 responses
      48 views
      0 reactions
      Last Post SEQadmin2  
      Working...