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

# In[1]:


get_ipython().run_line_magic('matplotlib', 'inline')
import numpy as np
import pandas as pd
import pymc3 as pm
import seaborn as sns


# In[6]:


fake_data = np.random.random((100, 301))
X = pd.DataFrame(fake_data, columns=['y']+['x%i' % i for i in range(300)])


# In[8]:


formula = 'y ~ ' + ' + '.join(X.columns.values[1:])
formula


# In[9]:


with pm.Model() as model:
    pm.glm.glm(formula, X)