如何 Save/Load 优化 GPy 回归模型

Question

我正在尝试保存优化的高斯过程模型以用于不同的脚本。我目前的想法是利用 GPy 的内置 to_dict 和 from_dict 函数将模型信息存储在 json 文件中。大致如下：

import GPy
import numpy as np
import json

X = np.random.uniform(-3.,3.,(20,1))
Y = np.sin(X) + np.random.randn(20,1)*0.05
kernel = GPy.kern.RBF(input_dim=1, variance=1., lengthscale=1.)

m = GPy.models.GPRegression(X, Y, kernel)

m.optimize(messages=True)
m.optimize_restarts(num_restarts = 10)

jt = json.dumps(m.to_dict(save_data=False), indent=4)
with open("j-test.json", 'w') as file:
    file.write(jt)

这一步没有问题，但是当我尝试使用加载模型信息时运行遇到了问题：

with open("j-test.json", 'r') as file:
    d = json.load(file)  # d is a dictionary

m2 = GPy.models.GPClassification.from_dict(d, data=None)

这给了我一个断言错误，因为“数据不是 None”，它是——或者至少我是这么认为的。

我对 GPy 和使用 jsons 真的很陌生，所以我真的不确定我哪里误入歧途了。我试着查看文档，但文档有点含糊，我找不到它的使用示例。我错过了 step/concept 吗？另外，这是存储和重新加载我的模型的最佳方式吗？对此的任何帮助将不胜感激！谢谢！

Answer 1

pickle 模块是你的好朋友！

import pickle
with open('save.pkl', 'wb') as file:
    pickle.dump(m, file)

您可以在以后的脚本中调用它：

with open('save.pkl', 'rb') as file:
    loaded_model = pickle.load(file)

Answer 2

没有建议将 Pickle 作为执行此操作的推荐方法。请参阅结尾部分的 here。以下是相同的示例。

# let X, Y be data loaded above
# Model creation:
m = GPy.models.GPRegression(X, Y)
m.optimize()
# 1: Saving a model:
np.save('model_save.npy', m.param_array)
# 2: loading a model
# Model creation, without initialization:
m_load = GPy.models.GPRegression(X, Y, initialize=False)
m_load.update_model(False) # do not call the underlying expensive algebra on load
m_load.initialize_parameter() # Initialize the parameters (connect the parameters up)
m_load[:] = np.load('model_save.npy') # Load the parameters
m_load.update_model(True) # Call the algebra only once
print(m_load)

如何 Save/Load 优化 GPy 回归模型

How to Save/Load Optimized GPy Regression Model

python

json

python-3.x

gaussian-process

gpy