Pymc MAP -如何提取最终值？

Question

我有一个很大的模型，适合pymc.MAP

M = pymc.MAP(model)
M.fit(iterlim=1000)

其中model是一个包含许多pymc对象的大列表，它在MCMC采样器中工作得很好。然而，我的问题是，当我尝试拟合上面的地图时，我不知道如何提取结果。这个模型非常大，所以我不打算发布整个文件，但例如，它肯定包含：

arctanA = np.empty(numbins, dtype=object)
C = np.empty(numbins, dtype=object)
for i in range(numbins):
    arctanA[i] = pymc.Uniform('arctanA_%i' %i, -math.pi/2, math.pi/2)
    C[i] = pymc.Uniform('C_%i' %i, 0, 2)

但是，当我尝试使用M.C.value或M.C_0.value提取上面的值时，它会给出一个AttributeError: 'MAP' object has no attribute 'C_0'，依此类推。那么我到底做错了什么呢？是否有什么特殊的方式需要我定义对象，以便稍后提取它们？

Answer 1

最佳值位于您指定的变量集合的元素中。例如，尝试使用C[0].value。