4858a3da74bc14eeaf724c9896bfe77df9451548,PyMC2/tests/test_norm_approx.py,test_norm_approx,check_draws,#test_norm_approx#,34

Before Change

        N = NormalApproximation(model)
        N.fit("fmin")
        draws = []
        for i in range(1000):
            N.draw()
            draws.append(hstack((N.alpha.value, N.beta.value)))
        draws = array(draws)
        plot(draws[:,0],draws[:,1],"k.")
        xlabel(r"$\alpha$")
        ylabel(r"$\beta$")

After Change

        N.fit("fmin")
        N.sample(1000)

        plot(N.alpha.trace(),N.beta.trace(),"k.")
        xlabel(r"$\alpha$")
        ylabel(r"$\beta$")
        
if __name__=="__main__":

Instances