Skip to content

Conversation

@BorisMuzellec
Copy link
Collaborator

What does your PR implement? Be specific.

This PR implements a method in DeseqStats to make volcano plots.

Example:

counts = load_example_data()
metadata = load_example_data('metadata')

dds = DeseqDataSet(counts=counts, metadata=metadata)
dds.deseq2()

ds = DeseqStats(dds)
ds.summary()

ds.plot_volcano(LFC_threshold = .25, pval_threshold = 0.05)

volcano_plot_example

Copy link
Collaborator

@umarteauowkin umarteauowkin left a comment

Choose a reason for hiding this comment

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

I approve as it is a good first draft ! In an ideal world, all the fixed kwargs should be modifiable, but in that case I think people would just recopy the code !

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.

3 participants