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

Somalier module: the label on y-axis is incorrect #1546

Closed
leonorpalmeira opened this issue Sep 9, 2021 · 4 comments
Closed

Somalier module: the label on y-axis is incorrect #1546

leonorpalmeira opened this issue Sep 9, 2021 · 4 comments
Labels
bug: module Bug in a MultiQC module

Comments

@leonorpalmeira
Copy link

Description of bug

In the Somalier module (MultiQC v1.9), the label on the y-axis reads: "Scaled mean depth on X" (which should be around 2 for females and 1 for males given that this is computed by Somalier as "x_depth_mean": 2 * stats[i].x_dp.mean / stats[i].gtdp.mean). However, the value which is plotted is X_mean_depth, see below:

Screenshot from 2021-09-09 14-36-34

I believe the y-label should be modified.

File that triggers the error

No response

MultiQC Error log

No response

@leonorpalmeira
Copy link
Author

From what I see in the current code, the label is still incorrect in the current code base:

Screenshot from 2021-09-09 14-47-16

@ewels ewels added the bug: module Bug in a MultiQC module label Jan 27, 2022
@ewels
Copy link
Member

ewels commented Jan 27, 2022

oof, so this is a difference between x_depth_mean and x_mean_depth? I can see why that was confused! 😅

Fixing the plot axis is easy, but it might be nicer to have the scaled version which is not sequencing-depth dependent. It should be possible to replicate the formula that you quoted from Somalier above.

With the example data I then get this:
somalier_sex_check_plot

Instead of the original:
somalier_sex_check_plot (1)

I guess that this was the original intention, so let's go with this..

@ewels ewels closed this as completed in 6374939 Jan 27, 2022
@ewels
Copy link
Member

ewels commented Jan 27, 2022

Ok, changes made in 6374939 - please let me know if you spot any issues with the code.

Thanks for the detailed report!

@leonorpalmeira
Copy link
Author

The code looks good to me. Thanks for fixing this and for reverting back to what was probably intended!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug: module Bug in a MultiQC module
Projects
None yet
Development

No branches or pull requests

2 participants