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

# In[2]:


import pymc3 as pm
import arviz as az
import numpy as np


# In[3]:


np.random.seed(42)
x = np.random.randn(100)


# In[4]:


with pm.Model() as m1:
    m = pm.Normal('m', 0, 1)
    l = pm.Normal('l', m, observed=x)


# In[9]:


with pm.Model() as m2:
    m = pm.Normal('m', 0, 1.1)
    l = pm.Normal('l', m, observed=x)


# In[6]:


with m1:
    tr1 = pm.sample()


# In[10]:


with m2:
    tr2 = pm.sample()


# In[11]:


az.compare({'model1': tr1, 'model2': tr2})