Pymc MAP - 如何提取最终值？

Question

我有一个相当大的模型适合pymc.MAP：

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

model 是一个包含许多 pymc 对象的大型列表，这些对象在 MCMC 采样器中运行良好。但是，我的问题是，当我尝试拟合上述 MAP 时，我不知道如何提取结果。模型很大，所以我不打算发布整个文件，但例如，它肯定包含：

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。