Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • zinky
    replied
    Originally posted by fanli View Post
    Are you getting 300-500 DE genes from a pooled variance model (e.g. DESeq or edgeR) with no replicates? Or is that based solely on the fold change?

    I'd be really surprised if you managed to get 300 hits with no replicates even between two very different cell types.
    yes, FoldChange is taken in consideration mainly. I also performed a test which was proposed by Audic S et al in 1997(ref:Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7:986–995. ). This statistic
    model combine the bayes and poisson distribution, it finally gives a p value to help screening DE features(fdr was also calculated). So based on this method, i got 300-500 DE genes while the correlation of two samples transcript profile is bigger than 0.97.

    Leave a comment:


  • SNPsaurus
    replied
    Originally posted by amitm View Post
    hi zinky,
    You say based on fold change. Did you have replicates? If yes then maybe you can apply a standard statistical testing like a t-test.
    If you are not expecting DE and seeing some (by fold change) then those might be due to sample to sample variation.
    A t-test would consider standard dev. and then give p-values
    Care should be taken not to interpret a t-test p-value that is above some threshold like 0.05 to mean that the samples are the same. The samples could, for instance, be very different but the experiment lacks the statistical power to detect the differences, or the variance is too high (as you say). But the opposite of differently expressed with statistical significance is not similarly expressed with statistical significance, it is "not differently expressed with statistical significance". There are tests that try to measure statistically if two distributions are the same, and this was the basis of a flurry of reports on problems with how statistics are used in many clinical studies, but the exact tests escape me!

    Leave a comment:


  • fanli
    replied
    Are you getting 300-500 DE genes from a pooled variance model (e.g. DESeq or edgeR) with no replicates? Or is that based solely on the fold change?

    I'd be really surprised if you managed to get 300 hits with no replicates even between two very different cell types.

    Leave a comment:


  • zinky
    replied
    hi amitm, i have no replicates but only two samples here. I know the way to detect DE features while my purpose is to demonstrate that there is no difference between them. You advice is useful for understanding but may not work for making a conclusion.
    thanks

    Leave a comment:


  • amitm
    replied
    hi zinky,
    You say based on fold change. Did you have replicates? If yes then maybe you can apply a standard statistical testing like a t-test.
    If you are not expecting DE and seeing some (by fold change) then those might be due to sample to sample variation.
    A t-test would consider standard dev. and then give p-values

    Leave a comment:


  • zinky
    replied
    Originally posted by fanli View Post
    Do you mean absolute or relative transcription levels? If absolute, then you'd need an external spike-in like the ERCCs. Relative, I'd say the lack of differentially expressed genes as even in a normal two-condition comparison the correlation values tend to be very high (>0.95).
    thanks for your suggestion. yes i mean relative transcripts level, because the library size may result in differ counts number( PS: there are alseo some other variables). After normalization, my data show correlation test << 0.05, and the correlation values is 0.97, but there are still some(300-500) DE genes based on Foldchange (abs(threshold) >2). How to explan this ?
    thanks!

    Leave a comment:


  • fanli
    replied
    Do you mean absolute or relative transcription levels? If absolute, then you'd need an external spike-in like the ERCCs. Relative, I'd say the lack of differentially expressed genes as even in a normal two-condition comparison the correlation values tend to be very high (>0.95).

    Leave a comment:


  • How to decide the two sample are the same in trancription level?

    Hi, guys.
    Is there any criterion to help us to draw a conclusion that two samples are the same or unchanged in trancription level based on RNA-seq data? Correlation test ? or No differential gene indentified ?

Latest Articles

Collapse

  • seqadmin
    Best Practices for Single-Cell Sequencing Analysis
    by seqadmin



    While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
    06-06-2024, 07:15 AM
  • 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

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 06:54 AM
0 responses
10 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-14-2024, 07:24 AM
0 responses
17 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-13-2024, 08:58 AM
0 responses
16 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-12-2024, 02:20 PM
0 responses
17 views
0 likes
Last Post seqadmin  
Working...
X