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

ATAC QC component #868

Merged
merged 32 commits into from
Dec 13, 2024
Merged

ATAC QC component #868

merged 32 commits into from
Dec 13, 2024

Conversation

VladimirShitov
Copy link
Collaborator

@VladimirShitov VladimirShitov commented Aug 27, 2024

Changelog

Add component for calculating QC metrics for ATAC data.

Issue ticket number and link

Contributes to #398

Checklist before requesting a review

  • I have performed a self-review of my code

  • Conforms to the Contributor's guide

  • Check the correct box. Does this PR contain:

    • Breaking changes
    • New functionality
    • Major changes
    • Minor changes
    • Documentation
    • Bug fixes
  • Proposed changes are described in the CHANGELOG.md

  • CI tests succeed!

Copy link
Member

@DriesSchaumont DriesSchaumont left a comment

Choose a reason for hiding this comment

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

Small comments, but looking good overall! 👍

src/qc/calculate_atac_qc_metrics/config.vsh.yaml Outdated Show resolved Hide resolved
src/qc/calculate_atac_qc_metrics/config.vsh.yaml Outdated Show resolved Hide resolved
mdata = mu.read(par["input"])

logger.debug("Making var names unique")
mdata.var_names_make_unique()
Copy link
Member

Choose a reason for hiding this comment

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

I think I want to avoid calling this here

atac.obs.drop(col, axis=1, inplace=True)

logger.info("Calculating QC metrics")
sc.pp.calculate_qc_metrics(atac, percent_top=None, log1p=False, inplace=True, layer=par["layer"])
Copy link
Member

@DriesSchaumont DriesSchaumont Oct 4, 2024

Choose a reason for hiding this comment

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

I want to avoid that things get overwritten by scanpy, could you pass a copy of the data?

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Did it manually

@VladimirShitov
Copy link
Collaborator Author

Tests are now failing due to incorrect scanpy version specification:
3.545 ERROR: No matching distribution found for scanpy~=1.10.4

Shouldn't be a problem on my side, I just use scanpy.yaml for installation now. Worth fixing elsewhere?

Copy link
Member

@DriesSchaumont DriesSchaumont left a comment

Choose a reason for hiding this comment

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

LGTM! Thanks @VladimirShitov

@DriesSchaumont DriesSchaumont merged commit 8669cb1 into main Dec 13, 2024
2 checks passed
@DriesSchaumont DriesSchaumont deleted the feature/ataq-qc branch December 13, 2024 10:06
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.

2 participants