Skip to content

Fix zero in each gene#98

Merged
BorisMuzellec merged 13 commits intomainfrom
fix_zero_in_each_gene
Mar 24, 2023
Merged

Fix zero in each gene#98
BorisMuzellec merged 13 commits intomainfrom
fix_zero_in_each_gene

Conversation

@BorisMuzellec
Copy link
Copy Markdown
Collaborator

@BorisMuzellec BorisMuzellec commented Mar 15, 2023

Reference Issue or PRs

Fixes #95

What does your PR implement? Be specific.

When each gene has at least one zero, it is not possible to compute size factors using the median of ratios method.
This PR implements an iterative MLE scheme to fit size factors when this happens (see #95).

@BorisMuzellec BorisMuzellec marked this pull request as ready for review March 15, 2023 16:03
@BorisMuzellec BorisMuzellec requested a review from a user March 15, 2023 16:03
Copy link
Copy Markdown
Collaborator

@maikia maikia left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Nice work. some minor comments before merge. Especially the test might need some improvement.

@BorisMuzellec BorisMuzellec requested a review from maikia March 21, 2023 14:44
Copy link
Copy Markdown
Collaborator

@maikia maikia left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@BorisMuzellec BorisMuzellec merged commit 67a9f02 into main Mar 24, 2023
@BorisMuzellec BorisMuzellec deleted the fix_zero_in_each_gene branch March 24, 2023 14:14
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

Successfully merging this pull request may close these issues.

[BUG] DeseqDataSet.deseq2() fails when "all genes contain a sample with a zero"

2 participants