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

Lineage specific plot trouble shooting #220

Open
Aishah86 opened this issue Mar 12, 2024 · 3 comments
Open

Lineage specific plot trouble shooting #220

Aishah86 opened this issue Mar 12, 2024 · 3 comments
Labels
bug Something isn't working

Comments

@Aishah86
Copy link

Dear,

We are looking into how freyja performs plotting on specific lineage (--lineages) on the visualization. We have plot our samples and the visualization display majorities of the lineages were group under "other". But if we look at the aggregate data "other" were not listed. Can you clarify what is happening, or any other config that we need to do.

Thank you

Specific lineage SARS-CoV-2

Hereby i attached aggregate data for Feb, for your reference

aggregate data February 2023

-AIshah-

@Aishah86 Aishah86 added the bug Something isn't working label Mar 12, 2024
@joshuailevy
Copy link
Collaborator

Hi Aishah!

Given the outputs provided, it seems you may just have a bunch of low level variants. When running demix, I'd recommend setting --eps to a small values (say, 0.00001) and for plotting you can lower the --thresh value as well. It may also be helpful to use the --config option to specify lineage groupings/colors, otherwise you can run out of lineage colors quite quickly and you end up with a lot of "Other".

Best,
Josh

@Aishah86
Copy link
Author

Do you mean the --thresh used during freyja plot? most of our data the plotting will be for data with < 0.5 abundance. if let say we put --thresh 0.1 , does it mean that lineages with >0.1 abundance will be plotted?

@joshuailevy
Copy link
Collaborator

Apologies- this seems to have evaded me. Yes, exactly. Under the hood this gets implemented via the prepLineageDict() function.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Development

No branches or pull requests

2 participants