我玩PYMC,在线性回归代码中,我通过随机装饰器生成了一个随机变量:
CTCallCenter
我想在MCMC迭代中绘制log_chisq的演变。但是,此变量的跟踪似乎不会存储...
因此我的问题是:
1)我在哪里可以找到pymc decorator关键字的定义或含义示例,例如plot或trace?
2)为什么stocastic(但观察到的)变量的痕迹不能保留为痕迹?
非常感谢
@pymc.stochastic(observed=True, trace=True)
def model(value = y_data, x_values = x_data, m = m_coef, n = n_coef, sigma = sigma):
value_theo = m*x_values + n
chi_sq = np.sum( np.square(value - value_theo) / np.square(sigma))
log_ChiSq = - chi_sq / 2.0
return log_ChiSq
答案 0 :(得分:1)
如果您直接实例化它,可以将任何传递给PyMC对象的参数传递给PyMC装饰器。例如,observed
,trace
和plot
都是Stochastic
可接受的所有参数。因此,要查看可以传递给@stochastic
的内容,请查看类docstring本身。