12.7. Use Omnibem to filter annotation dataΒΆ

This example shows how to create a new genomic feature data set from an annotation file.

from gdsctools import *
data = gdsctools_data("test_omnibem_genomic_alterations.csv.gz")
bem = OmniBEMBuilder(data)

Finally, create a MoBEM dataframe

mobem = bem.get_mobem()

# features
mobem = bem.get_mobem()

# Then, let us create a dataframe that is compatible with
# GenomicFeature. We just need to make sure the columns are correct
mobem[[x for x in mobem.columns if x!="SAMPLE"]]
gf = GenomicFeatures(mobem[[x for x in mobem.columns if x!="SAMPLE"]])

The final volcano plot

an = ANOVA(ic50_test, gf, verbose=False)
results = an.anova_all(animate=False)

Total running time of the script: ( 0 minutes 2.144 seconds)

Gallery generated by Sphinx-Gallery