In past several weeks, I spent a tremendous amount of time on reading literature about automatic parameter tuning in the context of Machine Learning (ML), most of which can be classified into two major categories, e.g. search and optimization. Searching mechanisms, such as grid search, random search, and Sobol sequence, can be somewhat computationally expensive. However, they are extremely easy to implement and parallelize on a multi-core PC, as shown in https://statcompute.wordpress.com/2019/02/03/sobol-sequence-vs-uniform-random-in-hyper-parameter-optimization. On the other hand, optimization algorithms, especially gradient-free optimizers such as Nelder–Mead simplex and particle swarm, are often able to quickly locate close-to-optimal solutions in cases that the global optimal is neither feasible nor necessary, as shown in https://statcompute.wordpress.com/2019/02/10/direct-optimization-of-hyper-parameter and https://statcompute.wordpress.com/2019/02/23/gradient-free-optimization-for-glmnet-parameters.
In the example below, another interesting approach, namely Bayesian optimization (https://arxiv.org/abs/1206.2944), is demonstrated and compared with CMA-ES (https://www.researchgate.net/publication/227050324_The_CMA_Evolution_Strategy_A_Comparing_Review), which is also a popular gradient-free optimizer based on the evolution strategy. As shown in the result, the output from Bayesian optimization is closer to the one from Nelder–Mead simplex and particle swarm. What’s more, Bayesian optimization is more consistent than CMA-ES among multiple trials in the experiment.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| cma_out <- cmaes::cma_es( | |
| par = 0.5, | |
| fn = function(x) grnn.optim(x, net, 4, 2019), | |
| lower = 0.1, upper = 1, | |
| control = list(fnscale = -1, mu = 20, lambda = 50)) | |
| #$par | |
| #[1] 0.5766267 | |
| #$value | |
| #[1] 0.8018076 | |
| bay_out <- rBayesianOptimization::BayesianOptimization( | |
| FUN = function(x) list(Score = grnn.optim(x, net, 4, 2019), Pred = 0), | |
| bounds = list(x = c(0.1, 1)), | |
| init_points = 5, n_iter = 20, | |
| acq = "ucb", verbose = F) | |
| # Best Parameters Found: | |
| #Round = 20 x = 0.5583 Value = 0.8019 |
Yet Another Blog in Statistical Computing 