@
user1572508提供的功能现在是PyMC的一部分，名称为
stochastic_from_data()或
Histogram()。然后，该线程的解决方案变为：
from pymc import *import matplotlib.pyplot as pltxtrue = 2 # unknown in the real applicationprior = rnormal(0,1,10000) # initial guess is inaccuratefor i in range(5): x = stochastic_from_data('x', prior) y = x*x obs = Normal('obs', y, 0.1, xtrue*xtrue + rnormal(0,1), observed=True) model = Model([x,y,obs]) mcmc = MCMC(model) mcmc.sample(10000) Matplot.plot(mcmc.trace('x')) plt.show() prior = mcmc.trace('x')[:]
