如何在 PyGPGO 中使用数组？

问题描述

我正在通过包 PyGPGO 研究贝叶斯优化，它在底层使用 PyMC3，将高斯过程拟合到底层函数并提出运动以最大化奖励。

到目前为止，根据标量值定义函数似乎很容易，尤其是在参数有限的情况下。但是，我不确定是否有办法提供向量/数组参数。

from pyGPGO.covfunc import squaredExponential
from pyGPGO.acquisition import Acquisition
from pyGPGO.surrogates.GaussianProcess import GaussianProcess
from pyGPGO.GPGO import GPGO

def f(x):
    return (np.sin(x))


sexp = squaredExponential()
gp = GaussianProcess(sexp)
acq = Acquisition(mode='ExpectedImprovement')
param = {'x': ('cont',[0,2 * np.pi])}

np.random.seed(23)
gpgo = GPGO(gp,acq,f,param)
gpgo.run(max_iter=20)

我想要的是：

#        var:  type,range
param = {'x': (np.array(3,10),2 * np.pi])}

如果这不可能，我可能会以 {x1: ...,x2: ...,} 的形式设置多个变量，然后通过函数定义本身相应地访问它们，尽管当参数空间很大时这似乎不是最理想的。

编辑：这是我的方法，但不起作用...

import numpy as np
from pyGPGO.covfunc import squaredExponential
from pyGPGO.acquisition import Acquisition
from pyGPGO.surrogates.GaussianProcess import GaussianProcess
from pyGPGO.GPGO import GPGO

param = {f'x_{i}': ('cont',1]) for i in range(3)}
param.update({f'y_{i}': ('cont',1]) for i in range(3)})

def f(p):
    v1 = np.array([p[f'x_{i}'] for i in range(3)])
    v2 = np.array([p[f'y_{i}'] for i in range(3)])
    return np.dot(v1,v2)


sexp = squaredExponential()
gp = GaussianProcess(sexp)
acq = Acquisition(mode='ExpectedImprovement')

np.random.seed(23)
gpgo = GPGO(gp,param)
gpgo.run(max_iter=20)

和错误：

Evaluation   Proposed point       Current eval.      Best eval.
---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-49-f58f602dca8c> in <module>()
     26 np.random.seed(23)
     27 gpgo = GPGO(gp,param)
---> 28 gpgo.run(max_iter=20)

1 frames
/usr/local/lib/python3.7/dist-packages/pyGPGO/GPGO.py in _firstRun(self,n_eval)
     86             s_param_val = list(s_param.values())
     87             self.X[i] = s_param_val
---> 88             self.y[i] = self.f(**s_param)
     89         self.GP.fit(self.X,self.y)
     90         self.tau = np.max(self.y)

TypeError: f() got an unexpected keyword argument 'x_1'

因此，将参数包装到字典中似乎比我预期的要棘手一些，因为对象 param 中的参数必须由函数 f 显式调用，而不是作为字典，p。

以前有人解决过这个问题吗？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）