Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

[Question]: Should one Re-run DESeq after removing genes with low values #15

Open
tamuanand opened this issue Nov 14, 2024 · 1 comment

Comments

@tamuanand
Copy link

Hi @federicomarini

I happened to find this behavior on the "Counts Overview" tab. Let's say I use "airway demo" data and "filter DDS by threshold on row sums of the counts" and choose 10 the first time

  • the resulting table shrinks to 22,008 entries
  • If I now change the threshold to 8 (or anything <10), the number of rows in the table does not increase (it remains the same at 22,008). Is this expected?

Does one have to re-run Step 3 in the Data Setup tab (run DESeq) after doing the above filtering (say using 10 as threshold with row sums)? I ask so because the ideal userguide has this - https://federicomarini.github.io/ideal/articles/ideal-usersguide.html
Additionally, [ideal](https://bioconductor.org/packages/3.20/ideal) has an option to include a filter step at the gene level by removing genes with low absolute or averages low values. After this, it might be possible to have to re-run the analysis in step 3 from the Data Setup panel

Thanks in advance.

@federicomarini
Copy link
Owner

After changing the threshold, you should still "re-create" the dds object in use for the app.

In theory, you should run DESeq() on the features you keep in your object.
I did not run any benchmark on what you obtain for the retained features if you keep the original values from the "full set" vs the re-done one, but I am expecting them to not vary too massively.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants