0.11.1

Bayesian Optimization in PyTorch

Introduction
Get Started
Tutorials

Key Features

Modular

Plug in new models, acquisition functions, and optimizers.

Built on PyTorch

Easily integrate neural network modules. Native GPU & autograd support.

Scalable

Support for scalable GPs via GPyTorch. Run code on multiple devices.

References

BoTorch: A Framework for Efficient Monte-Carlo Bayesian Optimization
@inproceedings{balandat2020botorch,
  title = {{BoTorch: A Framework for Efficient Monte-Carlo Bayesian Optimization}},
  author = {Balandat, Maximilian and Karrer, Brian and Jiang, Daniel R. and Daulton, Samuel and Letham, Benjamin and Wilson, Andrew Gordon and Bakshy, Eytan},
  booktitle = {Advances in Neural Information Processing Systems 33},
  year = 2020,
  url = {http://arxiv.org/abs/1910.06403}
}
Check out some other papers using BoTorch.

Get Started

  1. Install BoTorch:

    via Conda (strongly recommended for OSX):
    conda install botorch -c pytorch -c gpytorch -c conda-forge
    via pip:
    pip install botorch

  2. Fit a model:

    import torch
from botorch.models import SingleTaskGP
from botorch.fit import fit_gpytorch_mll
from botorch.utils import standardize
from gpytorch.mlls import ExactMarginalLogLikelihood

train_X = torch.rand(10, 2)
Y = 1 - torch.linalg.norm(train_X - 0.5, dim=-1, keepdim=True)
Y = Y + 0.1 * torch.randn_like(Y)  # add some noise
train_Y = standardize(Y)

gp = SingleTaskGP(train_X, train_Y)
mll = ExactMarginalLogLikelihood(gp.likelihood, gp)
fit_gpytorch_mll(mll)

  3. Construct an acquisition function:

    from botorch.acquisition import UpperConfidenceBound

UCB = UpperConfidenceBound(gp, beta=0.1)

  4. Optimize the acquisition function:

    from botorch.optim import optimize_acqf

bounds = torch.stack([torch.zeros(2), torch.ones(2)])
candidate, acq_value = optimize_acqf(
    UCB, bounds=bounds, q=1, num_restarts=5, raw_samples=20,
)
candidate  # tensor([0.4887, 0.5063])