在我正在编写的模型的一个步骤中,我必须计算一个量的误差函数。我想要做的看起来像这样:
from math import erf
import numpy as np
import pymc as pm
sig = pm.Exponential('sig', beta=0.1, size=10)
x = erf(sig ** 2)
这失败了,因为erf
它不适用于数组。我试过了:
@pm.deterministic
def x(sig=sig):
return [erf(s) for s in sig]
但没有成功,我知道可以通过以下方式获得结果:
np_erf = np.vectorize(erf)
x = np_erf((sig ** 2).value)
但这似乎不是正确的方法,因为它不会产生 apm.Deterministic
而只是产生np.array
. 我该怎么做呢?(PyMC 是 2.3 版)
编辑:为清楚起见,对上述示例进行了简化,以下是相关段落在实际代码中的外观。理想情况下,我希望这个工作:
mu = pm.LinearCombination('mu', [...], [...])
sig2 = pm.exp(mu) ** 2
f = 1 / (pm.sqrt(np.pi * sig2 / 2.0) * erf(W / sig2))
但如果消息失败TypeError: only length-1 arrays can be converted to Python scalars
。走np.vectorize
路线
np_erf = np.vectorize(erf)
f = 1 / (pm.sqrt(np.pi * sig2 / 2.0) * np_erf(W / sig2))
崩溃并显示相同的错误消息。列表理解
@pm.deterministic
def f(sig2=sig2):
return [1 / (pm.sqrt(np.pi * s / 2.0) * erf(W / s)) for s in sig2]
可以这样工作,但稍后会在此位置的代码中导致错误:
@pm.observed(plot=True)
def y(value=df['dist'], sig2=sig2, f=f):
return (np.log(np.exp(-(value ** 2) / 2.0 / sig2) * f)).sum()
错误是AttributeError: log
。
我已经使用数值近似计算了误差函数,这应该意味着一般设置是正确的。直接使用该erf
功能会更好更清晰。