figure() err = y.get_piecewise_pdf() - M3.as1DDistr().get_piecewise_pdf() err.plot() show() stop print "---", [K] + Y print M M.varschange(A, K) print M for i in xrange(n): M.varschange(Y[i], Y[i+1]) print M M.varschange(K, A) M.plot() #M.condition(Y[n], 2) print M M.eliminate(K) print M for i in xrange(n-1,-1,-1): M.eliminate(Y[i]) print M M.eliminate(A) print M figure() M.plot() print M.nddistr.pdf(linspace(0,2.5,100)) X0 = BetaDistr(1, 1) y = X0 * exp(A) y.summary()
figure() err = y.get_piecewise_pdf() - M3.as1DDistr().get_piecewise_pdf() err.plot() show() stop print "---", [K] + Y print M M.varschange(A, K) print M for i in xrange(n): M.varschange(Y[i], Y[i + 1]) print M M.varschange(K, A) M.plot() #M.condition(Y[n], 2) print M M.eliminate(K) print M for i in xrange(n - 1, -1, -1): M.eliminate(Y[i]) print M M.eliminate(A) print M figure() M.plot() print M.nddistr.pdf(linspace(0, 2.5, 100)) X0 = BetaDistr(1, 1) y = X0 * exp(A) y.summary()