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.

2 participants