bystrogenomics / bystro

Natural Language Search and Analysis of High Dimensional Genomic Data
Mozilla Public License 2.0
44 stars 14 forks source link

[prs] Distance-based PRS with corrected betas #565

Closed akotlar closed 2 months ago

akotlar commented 2 months ago

I've attached dave's output vs what I believe is ours. I will double check by re-running again, but this is roughly the difference I recall. Our PRS values are within 10% of Dave's at most, and the differences could be down to precision potentially. I spent a good deal of time tracking down the differences, and it seemed they likely happened in the final step, when adjusting dosages (the ancestry-adjusted allele frequencies were identical, as were the genotypes). I will look into this more in the future, but this PR is a strict improvement, and likely right.

Archive 2.zip

Summary by CodeRabbit

coderabbitai[bot] commented 2 months ago

Walkthrough

The recent changes significantly enhance the proteomics data analysis codebase, introducing new functionality for polygenic risk score calculations and improved data handling through streaming. Comprehensive documentation and logging capabilities have been added, promoting clarity and monitoring. Additionally, the structure of data classes has evolved to better manage covariates, ensuring robustness and accuracy in genetic data processing.

Changes

Files Change Summary
python/python/bystro/examples/...Query.ipynb Updates to execution counts, introduction of streaming functionality, added user credentials, and identified a variable initialization bug.
python/python/bystro/prs/handler.py, python/python/bystro/prs/messages.py, python/python/bystro/prs/prs_types.py Enhanced calculate_prs_scores function with new parameters, and expanded PRSJobData attributes to improve genetic data processing.
python/python/bystro/prs/listener.py Added logging functionality for improved tracking of PRS job operations.
python/python/bystro/prs/preprocess_for_prs.py Refactored code for clarity, added documentation, introduced constants for data handling, and improved PRS calculation logic.
python/python/bystro/utils/covariates.py New file defining structured data classes for covariates and experiment mappings to enhance data integrity and management.

Sequence Diagram(s)

sequenceDiagram
    participant User
    participant API
    participant PRSHandler
    participant PRSJobData
    participant Logger

    User->>API: Initiate PRS Job
    API->>PRSHandler: Call calculate_prs_scores()
    PRSHandler->>PRSJobData: Extract parameters
    PRSJobData->>PRSHandler: Return parameter data
    PRSHandler->>Logger: Log job initiation
    PRSHandler->>PRSHandler: Process data
    PRSHandler->>Logger: Log processing details
    PRSHandler->>API: Return results
    User->>API: Receive PRS results

πŸ‡ In the fields where we hop and play,
New features bloom, bright as day.
Data streaming flows like a stream,
Together we build, we share a dream.
With logs to guide us, we leap and bound,
In coding joy, our hearts are found! 🌼


Thank you for using CodeRabbit. We offer it for free to the OSS community and would appreciate your support in helping us grow. If you find it useful, would you consider giving us a shout-out on your favorite social media?

Share - [X](https://twitter.com/intent/tweet?text=I%20just%20used%20%40coderabbitai%20for%20my%20code%20review%2C%20and%20it%27s%20fantastic%21%20It%27s%20free%20for%20OSS%20and%20offers%20a%20free%20trial%20for%20the%20proprietary%20code.%20Check%20it%20out%3A&url=https%3A//coderabbit.ai) - [Mastodon](https://mastodon.social/share?text=I%20just%20used%20%40coderabbitai%20for%20my%20code%20review%2C%20and%20it%27s%20fantastic%21%20It%27s%20free%20for%20OSS%20and%20offers%20a%20free%20trial%20for%20the%20proprietary%20code.%20Check%20it%20out%3A%20https%3A%2F%2Fcoderabbit.ai) - [Reddit](https://www.reddit.com/submit?title=Great%20tool%20for%20code%20review%20-%20CodeRabbit&text=I%20just%20used%20CodeRabbit%20for%20my%20code%20review%2C%20and%20it%27s%20fantastic%21%20It%27s%20free%20for%20OSS%20and%20offers%20a%20free%20trial%20for%20proprietary%20code.%20Check%20it%20out%3A%20https%3A//coderabbit.ai) - [LinkedIn](https://www.linkedin.com/sharing/share-offsite/?url=https%3A%2F%2Fcoderabbit.ai&mini=true&title=Great%20tool%20for%20code%20review%20-%20CodeRabbit&summary=I%20just%20used%20CodeRabbit%20for%20my%20code%20review%2C%20and%20it%27s%20fantastic%21%20It%27s%20free%20for%20OSS%20and%20offers%20a%20free%20trial%20for%20proprietary%20code)
Tips ### Chat There are 3 ways to chat with [CodeRabbit](https://coderabbit.ai): - Review comments: Directly reply to a review comment made by CodeRabbit. Example: - `I pushed a fix in commit .` - `Generate unit testing code for this file.` - `Open a follow-up GitHub issue for this discussion.` - Files and specific lines of code (under the "Files changed" tab): Tag `@coderabbitai` in a new review comment at the desired location with your query. Examples: - `@coderabbitai generate unit testing code for this file.` - `@coderabbitai modularize this function.` - PR comments: Tag `@coderabbitai` in a new PR comment to ask questions about the PR branch. For the best results, please provide a very specific query, as very limited context is provided in this mode. Examples: - `@coderabbitai generate interesting stats about this repository and render them as a table.` - `@coderabbitai show all the console.log statements in this repository.` - `@coderabbitai read src/utils.ts and generate unit testing code.` - `@coderabbitai read the files in the src/scheduler package and generate a class diagram using mermaid and a README in the markdown format.` - `@coderabbitai help me debug CodeRabbit configuration file.` Note: Be mindful of the bot's finite context window. It's strongly recommended to break down tasks such as reading entire modules into smaller chunks. For a focused discussion, use review comments to chat about specific files and their changes, instead of using the PR comments. ### CodeRabbit Commands (invoked as PR comments) - `@coderabbitai pause` to pause the reviews on a PR. - `@coderabbitai resume` to resume the paused reviews. - `@coderabbitai review` to trigger an incremental review. This is useful when automatic reviews are disabled for the repository. - `@coderabbitai full review` to do a full review from scratch and review all the files again. - `@coderabbitai summary` to regenerate the summary of the PR. - `@coderabbitai resolve` resolve all the CodeRabbit review comments. - `@coderabbitai configuration` to show the current CodeRabbit configuration for the repository. - `@coderabbitai help` to get help. Additionally, you can add `@coderabbitai ignore` anywhere in the PR description to prevent this PR from being reviewed. ### CodeRabbit Configuration File (`.coderabbit.yaml`) - You can programmatically configure CodeRabbit by adding a `.coderabbit.yaml` file to the root of your repository. - Please see the [configuration documentation](https://docs.coderabbit.ai/guides/configure-coderabbit) for more information. - If your editor has YAML language server enabled, you can add the path at the top of this file to enable auto-completion and validation: `# yaml-language-server: $schema=https://coderabbit.ai/integrations/schema.v2.json` ### Documentation and Community - Visit our [Documentation](https://coderabbit.ai/docs) for detailed information on how to use CodeRabbit. - Join our [Discord Community](https://discord.com/invite/GsXnASn26c) to get help, request features, and share feedback. - Follow us on [X/Twitter](https://twitter.com/coderabbitai) for updates and announcements.