Originally posted by jwfoley
View Post
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
-
It's interesting that your number of up and down genes is so different with the two approaches. I wonder if it is actually a "normalization" difference. If I were you I would look at some MA plots or smear plots coloured for each of the methods and see if you can see anything the might push you to choose one over the other. Nothing better than visualising your data in my opinion
Comment
-
Originally posted by A Oshlack View PostIt's interesting that your number of up and down genes is so different with the two approaches. I wonder if it is actually a "normalization" difference. If I were you I would look at some MA plots or smear plots coloured for each of the methods and see if you can see anything the might push you to choose one over the other. Nothing better than visualising your data in my opinion
i did have MA plots and smear plots generated by those two packages, but i do not ubderstand what you have inferred.
would you please give me more details.
thank you
Comment
-
As I'm sure you know, an MA plot is just a plot of the raw data. In the normalization assumptions, at least for TMM in edgeR, we assume that the "cloud" of points in an MA plot should be centred at M=0. You have used a two fold cut-off as well as the statistical test so I wonder if the difference between methods that you are using is actually more dependent on the normalization that you have chosen. Are the DE gene in the MA plot equal distance up and down from where you think that the normalization line should sit. Often your eyes are very good at assessing this if you generate the right plots. Hope that helps.
Comment
-
I believe that replicates are very important to have good quality results. RNA-seq is becoming cheaper and cheaper but still quite expensive for small labs. In this case I also believe that RNAseq without replicates could be used as screening and then confirm by replicating qRT-PCR and based you conclusion on these results.
Comment
-
I think both EdgeR and DESeq are both pretty bad. Not sure how this community came to believe that "borrowing variance" from genes with similar average expression makes sense. It does not make biological sense and it is mathematically provably wrong. Therefore, both EdgeR, DESeq, Limma produce more false positives and more false negatives than a more robust and reliable statistical test. If it made good sense, respected statistical leaders like SAS and Partek would have adopted this methodology, and neither company has.
Comment
-
Originally posted by rfilbert View PostI think both EdgeR and DESeq are both pretty bad. Not sure how this community came to believe that "borrowing variance" from genes with similar average expression makes sense. It does not make biological sense and it is mathematically provably wrong. Therefore, both EdgeR, DESeq, Limma produce more false positives and more false negatives than a more robust and reliable statistical test. If it made good sense, respected statistical leaders like SAS and Partek would have adopted this methodology, and neither company has.
Metapress is a fast growing digital platform that helps visitors to answer questions, solve problems, learn new skills, find inspiration and provide the latest Technology news.
http://www.ncbi.nlm.nih.gov/pmc/arti...rtype=abstract.
etc
I don't buy this argument of "SAS and Partek don't do it, therefore its wrong".
Comment
-
Also, I think your claims about how DESeq works are somewhat inaccurate.
The default for DESeq is to calculate both the variance for each gene and the variance by borrowing from other genes of similar expression and then take the maximum of the two. Since one is taking the maximum of the two values, this would reduce the number of false positives, but increase the number of false negatives. Its easy enough to alter DESeq so that it uses one or the other so that if you really feel that passionate about not borrowing variance, than you don't have to.Last edited by chadn737; 12-21-2012, 11:54 AM.
Comment
-
The same goes for edgeR - you don't need to use the moderated dispersions if you don't want to.
I usually prefer SAMSeq (a non-parametric method) for DE analysis when there are enough replicates, but the nice thing with DESeq and edgeR is that you can consider complex designs in your analysis.
Comment
-
I would definitely recommend Partek Flow for the following reasons:
1. They fit 5 different distribution assumptions and use the best fit for each gene or transcript. This gives more statistical power and more biological meaning.
2. There are no limit to the number of factors (handles any type of experiment design)
3. It has a really easy to use point & click web-based GUI.
4. They have excellent technical support.
Comment
-
There are multiple options for normalization, but I believe the default option is to simply normalize to the total number of reads for each sample. I don't think any normalization based on the length of the transcript (like RPKM) matters as for this analysis you are comparing the same transcript in different groups of samples.
Comment
-
1. They fit 5 different distribution assumptions and use the best fit for each gene or transcript. This gives more statistical power and more biological meaning.
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 05-02-2024, 08:06 AM
|
0 responses
17 views
0 likes
|
Last Post
by seqadmin
05-02-2024, 08:06 AM
|
||
Started by seqadmin, 04-30-2024, 12:17 PM
|
0 responses
20 views
0 likes
|
Last Post
by seqadmin
04-30-2024, 12:17 PM
|
||
Started by seqadmin, 04-29-2024, 10:49 AM
|
0 responses
27 views
0 likes
|
Last Post
by seqadmin
04-29-2024, 10:49 AM
|
||
Started by seqadmin, 04-25-2024, 11:49 AM
|
0 responses
28 views
0 likes
|
Last Post
by seqadmin
04-25-2024, 11:49 AM
|
Comment