Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Limma Voom Differential Splicing diffSplice() error

    Hi,

    I've been following the user's guide for Limma voom to perform differential splicing analysis and I have encountered a few problems.

    Firstly, In chapter 15.4 the following line of code produces and error:

    Code:
    dge <- dge[A>10,,keep.lib.sizes=FALSE]
    Error in object[[a]][i, , drop = FALSE] : incorrect number of dimensions
    If I perform filtering using the following code which I believe does something similar

    Code:
     
    isexpr <- rowSums(dge$counts) > 10
    x <- x[isexpr,]
    dge <- calcNormFactors(x)
    And then continue with the following as stated in the documentation I get an error:

    Code:
    v <- voom(dge,design,plot=TRUE)
    fit <- lmFit(v, design)
    
    ex <- diffSplice(fit, geneid="GeneID")
    Error: could not find function "diffSplice"
    Is there something (perhaps very stupid) that I'm doing wrong??

    Thanks in advance
    Last edited by Anomilie; 08-04-2014, 08:17 PM.

  • #2
    Please provide sessionInfo() so we can better assist you...

    Code:
    dge <- dge[A>10[COLOR="Red"],,[/COLOR]keep.lib.sizes=FALSE]
    Error in object[[a]][i, , drop = FALSE] : incorrect number of dimensions
    I see a double "," in there. That's causing your number of dimension error.

    Code:
    Error: could not find function "diffSplice"
    Update Limma if you have an old version.

    Comment


    • #3
      The double quote is directly copied from the documentation and is apparently requirement for subsetting the DGEList object.

      the following code sorted out that issue (in case someone else has similar problems)

      Code:
      dge <- DGEList(counts=raw_data, genes=geneID)
      In regards to the diffSplice error, the error message seems fairly obvious now, I had to update limma to version 3.20.

      Thanks for your reply Cytosine

      While on the topic of diffSplice, would you know the statistical difference between using that function and performing the DE with the standard limma voom pipeline using eBayes(fit)?
      Last edited by Anomilie; 08-04-2014, 11:42 PM.

      Comment


      • #4
        Originally posted by Anomilie View Post
        While on the topic of diffSplice, would you know the statistical difference between using that function and performing the DE with the standard limma voom pipeline using eBayes(fit)?
        It is exactly equivalent to a voom eBayes pipeline, but adapted to assess differential splicing instead of differential expression.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Latest Developments in Precision Medicine
          by seqadmin



          Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

          Somatic Genomics
          “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
          05-24-2024, 01:16 PM
        • seqadmin
          Recent Advances in Sequencing Analysis Tools
          by seqadmin


          The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
          05-06-2024, 07:48 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 05-24-2024, 07:15 AM
        0 responses
        15 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-23-2024, 10:28 AM
        0 responses
        18 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-23-2024, 07:35 AM
        0 responses
        21 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-22-2024, 02:06 PM
        0 responses
        10 views
        0 likes
        Last Post seqadmin  
        Working...
        X