Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • DEXSeq output

    Hi!

    I've used the DESeq package and now I'm trying the DEXSeq one.

    In the DESeq package I can access the baseMean, baseMeanA, baseMeanB values after doing the analysis as a output result, but using the DEXSeq package, I only get the meanBase value.


    OUTPUT DEXSeq:

    > res1 <- DEUresultTable(ecs)
    > head(res1)
    geneID exonID dispersion pvalue padjust meanBase log2fold(all/sl)
    Smp_156060:E015 Smp_156060 E015 0.65669650 4.329870e-15 6.934596e-12 19.147740 -2.023423
    Smp_141410:E001 Smp_141410 E001 0.48777320 2.455751e-06 7.733547e-04 7.941987 -1.977033
    Smp_080520:E001 Smp_080520 E001 0.08365269 1.022540e-08 5.789747e-06 13.344642 -1.944765
    Smp_194240:E001 Smp_194240 E001 0.33620393 4.173149e-06 1.285307e-03 2.973169 -1.930388
    Smp_034190:E001 Smp_034190 E001 0.04309263 0.000000e+00 0.000000e+00 30.344916 -1.691994
    Smp_091390:E001 Smp_091390 E001 0.07264235 0.000000e+00 0.000000e+00 15.738090 -1.682837

    OUTPUT DESeq:

    > res1 <- nbinomTest( cds, "sl", "all" )
    > head(res1)
    id baseMean baseMeanA baseMeanB foldChange log2FoldChange pval padj
    1 Schisto_mansoni.mitochondria_rRNA1 1.081107e+04 2.668263 24321.567008 9115.130 13.15405 8.314941e-11 7.132301e-09
    2 Schisto_mansoni.mitochondria_rRNA2 7.196552e+03 1.820197 16189.965875 8894.625 13.11872 7.980118e-06 5.303116e-05
    3 Sm_mito_tRNA10.1 2.092785e+00 0.000000 4.708767 Inf Inf 3.392875e-01 5.555647e-01
    4 Sm_mito_tRNA11.1 7.391640e-01 0.000000 1.663119 Inf Inf 8.837895e-01 1.000000e+00
    5 Sm_mito_tRNA12.1 6.149918e-01 0.000000 1.383732 Inf Inf 1.000000e+00 1.000000e+00
    6 Sm_mito_tRNA13.1 9.076016e-01 0.000000 2.042104 Inf Inf 8.429498e-01 1.000000e+00
    How can I get the baseMeanA baseMeanB values using the DEXSeq package?

    Thanks!

    []s

    Mariana

  • #2
    Hi Mariana,

    If you are interested in the exon expression, you could just calculate the mean accross the raw counts. If you are interested in the fitted splicing coefficients (what is plotted in plotDEXSeq) you could have a look at the function estimatelog2FoldChanges to get this coefficients, especifically you have to set the parameter getOnlyEffects=TRUE.

    Alejandro

    Comment


    • #3
      DEXseq meanBaseA meanBaseB

      Could you please post some sample code how to extract the relevant coefficients in order to obtain proper meanBaseA / meanBaseB values in the DEUresultTable.
      I'm new to DEXseq and would highly appreciate your input.

      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, Yesterday, 06:37 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      9 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      50 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      67 views
      0 likes
      Last Post seqadmin  
      Working...
      X