In [1]:

    

%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

import seaborn as sns
import tssb
reload(tssb)
from tssb.measure import create_gaussian_trm
from tssb.util import ScipySampler
import scipy.stats as stats


    

In [129]:

    

tssb = create_gaussian_trm(1, 
                           np.eye(2)*4, 
                           ScipySampler(stats.multivariate_normal([0, 0], np.eye(2)*0.0000001)),
                           25,
                           0.5,
                           1)


    

In [135]:

    

N = 1000
parameters = tssb.sample(N)
X = []
for i in xrange(N):
    X.append(stats.multivariate_normal(mean=parameters[i][1][0], cov=np.eye(2) / 1000.0).rvs())
X = np.array(X)
x, y = X[:, 0], X[:, 1]
scatter(x, y, alpha=0.5);
xlim([-10, 10]);
ylim([-10, 10]);