#!/usr/bin/env python
# coding: utf-8

# # One Drug One Feature Analysis (ODOF)

# Similarly to notebook 1 (introduction) and 3 (one drug analysis), we use the ANOVA class but call another function called anova_one_drug_one_feature

# In[1]:


get_ipython().run_line_magic('pylab', 'inline')
matplotlib.rcParams['figure.figsize'] = (10,6)


# In[2]:


from gdsctools import ANOVA, ic50_test


# In[3]:


an = ANOVA(ic50_test)


# To get the drug idenfiers or genomic features, use the following code

# In[4]:


an.drugIds[0:5]


# In[5]:


an.feature_names[0:6]


# **note:** the first 3 feature names are not real features and should be ignored for now.

# In[6]:


results = an.anova_one_drug_one_feature(1047, 'BRAF_mut', 
                              show=True)


# Note that here this is a PANCAN analysis so all tissues are used. The MSI is also
# a feature used in the regression.