Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Gordon Smyth
    replied
    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.

    Leave a comment:


  • Anomilie
    replied
    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.

    Leave a comment:


  • Cytosine
    replied
    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.

    Leave a comment:


  • 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.

Latest Articles

Collapse

  • 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, Today, 07:35 AM
0 responses
3 views
0 likes
Last Post seqadmin  
Started by seqadmin, Yesterday, 02:06 PM
0 responses
8 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-14-2024, 07:03 AM
0 responses
27 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-10-2024, 06:35 AM
0 responses
47 views
0 likes
Last Post seqadmin  
Working...
X