Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • DESeq2 error

    Hi all,

    New to DESeq2 and getting error consistently while running this line:

    dds <- DESeq(dds)

    The error message:

    estimating size factors
    estimating dispersions
    same number of samples and coefficients to fit, estimating dispersion by treating samples as replicates
    gene-wise dispersion estimates
    mean-dispersion relationship
    Error in lfproc(x, y, weights = weights, cens = cens, base = base, geth = geth, :
    newsplit: out of vertex space
    Calls: DESeq ... estimateDispersionsFit -> localDispersionFit -> locfit -> lfproc -> .C
    In addition: There were 18 warnings (use warnings() to see them)
    Execution halted


    Any ideas? Thanks!
    Last edited by stormin; 09-06-2014, 11:38 AM.

  • #2
    Try manually specifying the fitType as parametric (dds <- DESeq(dds, fitType="parametric")). I'm actually not sure why it's using local fitting to begin with (last I looked, parametric was the default).

    Comment


    • #3
      Hmm, that didn't work unfortunately. Here is my R script with the same error message. I'm new to R, so sorry in advance if there is a very obvious mistake. Thanks!


      #!/usr/bin/Rscript

      library('DESeq2')

      directory <- getwd()
      sink("DESeq_output.txt")
      sampleTable <- read.csv("sampleTable.csv")

      dds <- DESeqDataSetFromHTSeqCount(sampleTable=sampleTable, directory=directory, design = ~ condition)
      colData(dds)$condition <- relevel(colData(dds)$condition, "noCoverage")
      dds <- DESeq(dds, fitType="parametric")
      results <- results(dds)
      results <- results[order(results$padj), ]
      write.csv(results,file="DESeq_results.csv")
      sink()

      detach("packageESeq2")

      quit()
      Same error message plus addition 18 warnings:
      1: In estimateDispersionsFit(object, fitType = fitType, quiet = quiet) :
      the parametric fit of dispersion estimates over the mean of counts
      failed, which occurs when the trend is not well captured by the
      function y = a/x + b. A local regression fit is automatically performed,
      and the analysis can continue. You can specify fitType='local' or 'mean'
      to avoid this message if re-running the same data.
      When using local regression fit, the user should examine plotDispEsts(dds)
      to make sure the fitted line is not sharply curving up or down based on
      the position of individual points.
      2: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      3: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      4: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      5: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      6: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      7: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      8: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      9: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      10: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      11: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      12: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      13: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      14: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      15: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      16: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      17: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      18: In lfproc(x, y, weights = weights, cens = cens, base = base, ... :
      procv: no points with non-zero weight
      Content of the sampleTable.csv
      sampleName,fileName,condition
      2hnoCoverage,counts_Coverage.txt,Coverage
      2hCoverage,counts_noCoverage.txt,noCoverage
      Last edited by stormin; 09-07-2014, 09:45 AM.

      Comment


      • #4
        That's actually a different error. So it seems that a parametric fit doesn't work well for this dataset. One possibility would be to use local fitting and increase the number of vertexes. I don't think there's a way to do that with DESeq2 (there was with DESeq), though you could modify the code a bit.

        Your best bet would be to post this to the bioconductor email list and CC Michael Love, the DESeq2 maintainer. Please include the output of sessionInfo() in your email. It may turn that you'll need to send one of the DESeq2 authors your dataset so they can have a look.

        As an aside, I probably wouldn't put too much effort into an experiment where you're just comparing two samples against each other (as appears to be the case from your sampleTable).

        Comment


        • #5
          Can you include the sessionInfo()?

          As the local fit isn't working, you can just use:

          DESeq(dds, fitType="mean")

          How many rows/genes do you have?

          Comment


          • #6
            Hi Devon, so what I am trying to do is to use DESeq2 to compare two data manipulation methods. The same RNAseq fastq dataset was used in tophat. Condition1 used coverage search (which takes a lot longer to run) while Condition2 used no coverage search. I suspect the additional coverage based junction search will not significantly impact reads and thus DE results. Although I still want to check.

            Perhaps the reason why DESeq is having so much trouble with this comparison is because the counts are mostly identical between the two conditions?

            To Michael, the mean fit worked and I will delve more into the manual to see why. However, the log2foldchange column doesn't make sense. The value is way too small for almost all of the genes with some being 0!

            Zach
            Last edited by stormin; 09-07-2014, 10:00 AM.

            Comment


            • #7
              hi Zach,

              Check out the DESeq2 vignette. We explain the new methods there (e.g. Fig 1). For even more details, check out the DESeq2 manuscript, which is linked on the front page of the vignette. Unreliable log fold changes (not supported by information) are shrunk towards 0.

              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, 03-27-2024, 06:37 PM
              0 responses
              12 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-27-2024, 06:07 PM
              0 responses
              11 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              52 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