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
          Current Approaches to Protein Sequencing
          by seqadmin


          Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
          04-04-2024, 04:25 PM
        • 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

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 04-11-2024, 12:08 PM
        0 responses
        31 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 10:19 PM
        0 responses
        33 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 09:21 AM
        0 responses
        28 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-04-2024, 09:00 AM
        0 responses
        53 views
        0 likes
        Last Post seqadmin  
        Working...
        X