Sample code to call dice.stan with input data

dice_sample.py

#!/usr/bin/env python                                                                                                                                                                                              

import pystan
from numpy import random
import matplotlib.pyplot as plt

m = 6
N = 1000

prob = [1/float(m) for i in range(m)]

x = random.multinomial(N,prob)

dice_data = {'m': m, 'x': x}

fit = pystan.stan(file='dice.stan', data=dice_data, iter=100000, chains=1)

print(fit)
fit.plot()
plt.show()