Search
Question: DESeq2: How to find out how well the model fits the data?
0
gravatar for lina.faller
12 months ago by
lina.faller0 wrote:

Hi all,

I am using DESeq2 to calculate differentially expressed genes. Is there a good way to determine how well the model fits the data? I'd appreciate it if you could share any insight or resources.

Thanks!
~Lina

ADD COMMENTlink modified 12 months ago by Bio_Ram0 • written 12 months ago by lina.faller0

Hi!!

I have the same question. I've got several factors and I would like to know if using all of them creates a better model. Did you find anything about it? I've been looking for an answer, but still not luck.... Thanks!!

ADD REPLYlink written 23 days ago by ceboral0

You seem to be thinking that you are fitting a single model. You aren't. Instead you are simultaneously fitting thousands of models, so asking 'how well a model fits' is in some sense nonsensical. Which model might you be asking about?

In general people just fit a model containing any nuisance variables that may affect the gene expression, and call it good. If you really care, you can test all of the nuisance variables you are including in your model and drop those that aren't 'significant enough', where by that I mean those variables that are either not significant for any gene, or for only a small subset of genes, when testing using a likelihood ratio test.

The downside of including too many nuisance variables is that you may be wasting degrees of freedom and reducing power, or that you may be including variables that are not orthogonal, which can be problematic. Otherwise an overspecified model isn't that big of a deal.

 

ADD REPLYlink written 23 days ago by James W. MacDonald46k

Yes, I know that a different model is tested for each gene considered, but I am testing two different factors and I would like to know which one fits the expression better and I cannot find any Akaike implementation or bayesian information criterion (I cannot rely only in the wald test of the factors because the models aren't nested). 

ADD REPLYlink written 16 days ago by ceboral0
0
gravatar for Bio_Ram
12 months ago by
Bio_Ram0
Bio_Ram0 wrote:

The workflows below will give you better insights 

https://www.bioconductor.org/help/workflows/rnaseqGene/

https://scilifelab.github.io/courses/rnaseq/labs/DEseq2

ADD COMMENTlink written 12 months ago by Bio_Ram0

This answer is not helpful at all. There is no where in these manuals that tells us how well the model fits the data

ADD REPLYlink written 5 months ago by rrcutler50
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 320 users visited in the last hour