Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Marianna85
    replied
    Hi Puggie,
    thanks for your reply.

    Originally posted by puggie View Post
    E.g. select 5-10.000 genes that are expressed in all samples and then normalize by these
    what do you mean normalize by 5-10.000 gene, how? you mean RPKM normalization?

    Originally posted by puggie View Post
    If samples are of different origin this approach may prove insufficient.
    the two samples I have come from the same tissue but different maturation stage. But several genes seem to be extremely DE (0 vs 100.000 reads!)

    Anyway RPKM normalization makes sense only if you don't need to compare different samples.

    Marianna

    Leave a comment:


  • puggie
    replied
    I work on whole-transcriptome also. I believe there is a lot of alternatives for normalizing this kind of data. I have tried edgeR but it did not produce meaningfull normalization by human eyes. What I do instead is based on FPKM/RPKM counts e.g. via Cufflinks. E.g. select 5-10.000 genes that are expressed in all samples and then normalize by these. If samples are of different origin this approach may prove insufficient.

    Leave a comment:


  • Marianna85
    replied
    Originally posted by Simon Anders View Post
    No matter how you normalize, all genes will show huge differences. And you will have no way to figure out whether you might have seen as huge differences has you sequenced two eggs from the same developmental stage (this is well possible: maybe mollusc eggs differ a lot from each other, or maybe your wet-lab protocol is unstable). This is, of course, just another example while doing an experiment without at least duplicates is simply bad science and a waste of time and money (and by now, no decent journal will accept such studies any more, I hope), but I am repeating, what I already said often in previous threads.

    Sorry for the pessimism, but quite frankly, I doubt that it's worth putting much more effort into this analysis. What exactly had you hoped to find?
    The experiment is mainly a transcriptome characterization...this is the reason why we didn't plan any replicates.In order to study the whole transcriptome, in my opinion, a simple normalization intra-library is enough.
    But I also hoped to find few genes differentially expressed, without a p-values and without statistical power, just to find transcripts (maybe 5-10) which show a huge difference between the two conditions. Just a first step to be validated, no more. To do this I need an inter-library normalization.
    Considering the data, do you think that this objective is too ambitious?
    Marianna

    Leave a comment:


  • Simon Anders
    replied
    No matter how you normalize, all genes will show huge differences. And you will have no way to figure out whether you might have seen as huge differences has you sequenced two eggs from the same developmental stage (this is well possible: maybe mollusc eggs differ a lot from each other, or maybe your wet-lab protocol is unstable). This is, of course, just another example while doing an experiment without at least duplicates is simply bad science and a waste of time and money (and by now, no decent journal will accept such studies any more, I hope), but I am repeating, what I already said often in previous threads.

    Sorry for the pessimism, but quite frankly, I doubt that it's worth putting much more effort into this analysis. What exactly had you hoped to find?

    Leave a comment:


  • Marianna85
    replied
    But...the two samples correspond each to a different stage...Maybe the expression profile between the two conditions is huge...and in this case, which is the best solution in order to normalize the raw count?

    Leave a comment:


  • Marianna85
    replied
    this is a RNAseq experiment in mollusc eggs. Raw reads after orthologs clustering. Yes, there are some points with a huge different among the two libraries. Is this the reason why I obtained so different normalization factors?

    Leave a comment:


  • Simon Anders
    replied
    Wow, this is quite a mess. Huge amounts of points with 0 reads in one and more than 10,000 reads in the other library. What have you sequenced there? If this is ordinary RNA-Seq data, things went quite wrong.

    Leave a comment:


  • Marianna85
    replied

    Leave a comment:


  • Marianna85
    replied
    sorry again...i think it's not ok

    Leave a comment:


  • Marianna85
    replied
    sorry

    Leave a comment:


  • Simon Anders
    replied
    ".emf"? That's Windows extended metafile, right? Haven't seen this graphics file format in ten years, and frankly, I have no idea how to open it. Could you use something more common, please, maybe png?

    Leave a comment:


  • Marianna85
    replied
    of course! I defined the cds rows, not the columns.
    So this is the plot...



    something strange in your opinion??

    Leave a comment:


  • Simon Anders
    replied
    Sorry, I made a type. It's

    Code:
    plot( log10( 1 + counts(cds)[,1] ), log10( 1 + counts(cds)[,2] ), pch="." )

    Leave a comment:


  • Simon Anders
    replied
    This will be hard to debug via the forum. You may need to get some local help.

    To try one thing: If you simply type "counts(cds)", you get your table of raw counts (or, if you just want the first 100 lines, try "head( counts(cds), 100 )". Check whether they make sense.

    Leave a comment:


  • Marianna85
    replied
    Hi Simon,
    the plot seems empty...
    may I change the axis scale?

    Leave a comment:

Latest Articles

Collapse

  • seqadmin
    An Introduction to the Technologies Transforming Precision Medicine
    by seqadmin


    In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
    01-27-2025, 07:46 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 09:30 AM
0 responses
15 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-05-2025, 10:34 AM
0 responses
23 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-03-2025, 09:07 AM
0 responses
26 views
0 likes
Last Post seqadmin  
Started by seqadmin, 01-31-2025, 08:31 AM
0 responses
34 views
0 likes
Last Post seqadmin  
Working...
X