''' Created on Nov 24, 2009 @author: johnsalvatier ''' import model2 import pymc import multichain_mcmc as mc from pylab import * import pydevd pydevd.set_pm_excepthook() sampler = mc.AmalaSampler(model2.model) history, time = sampler.sample(nChains = 5, ndraw = 3000, maxGradient = 100) print time mc.show_samples(plot, history, ('mean', 'sd'))
''' Created on Nov 24, 2009 @author: johnsalvatier ''' import model1 import pymc import multichain_mcmc as mc from pylab import * import numpy import pydevd pydevd.set_pm_excepthook() sampler = mc.AmalaSampler(model1.model) history, time = sampler.sample(nChains = 5, ndraw = 3000, maxGradient = 100) print time mc.show_samples(plot, history, ('a','b', 'sd'))