Bayesian Optimization of Trading Strategy Parameters

Developing a trading strategy without automatic parameter optimization is a lottery. Manual tuning takes weeks, grid search takes millions of runs, and the result is still suboptimal. Bayesian optimization (Wikipedia) with Optuna and TPE sampler solves the problem in 100 iterations instead of 10,000. We have implemented it in dozens of projects — from HFT bots to DeFi arbitrageurs. Let's break down a real case and show how it works.

Parameter optimization is a key stage of any trading system. Without it, the strategy remains raw. We solve this problem using Bayesian optimization on Optuna — the industry standard for hyperparameter search. Our Bayesian optimization system uses the TPE sampler for efficient parameter search. Our experience: over 10 years in trading and 50+ completed projects.

How does Bayesian parameter optimization work?

Bayesian optimization uses a surrogate model (Gaussian Process) to approximate the backtest function. After each evaluation, the model updates its representation of the dependence of parameters on the metric. The Acquisition Function (e.g., Expected Improvement) determines where to search next, balancing exploration and exploitation. This avoids useless runs and reduces search time by 5–10 times compared to Random Search. In a typical project, savings on computing resources exceed $10,000. The TPE sampler is robust to noisy metrics, which is critical for backtesting.

Implementation with Optuna

import optuna from optuna.samplers import TPESampler import pandas as pd optuna.logging.set_verbosity(optuna.logging.WARNING) class BayesianOptimizer: def __init__( self, backtest_fn: callable, n_trials: int = 100, n_startup_trials: int = 10, n_jobs: int = 1, metric: str = 'sharpe_ratio', direction: str = 'maximize', ): self.backtest_fn = backtest_fn self.n_trials = n_trials self.n_startup = n_startup_trials self.n_jobs = n_jobs self.metric = metric self.direction = direction self.results_log = [] def create_objective(self, param_space: dict): def objective(trial: optuna.Trial) -> float: params = {} for name, spec in param_space.items(): if spec['type'] == 'int': params[name] = trial.suggest_int(name, spec['low'], spec['high']) elif spec['type'] == 'float': params[name] = trial.suggest_float(name, spec['low'], spec['high']) elif spec['type'] == 'categorical': params[name] = trial.suggest_categorical(name, spec['choices']) elif spec['type'] == 'log': params[name] = trial.suggest_float(name, spec['low'], spec['high'], log=True) try: metrics = self.backtest_fn(params) value = metrics.get(self.metric, float('-inf')) n_trades = metrics.get('total_trades', 0) if n_trades < 15: value = value * n_trades / 15 max_dd = abs(metrics.get('max_drawdown_pct', 0)) if max_dd > 40: value = value * (40 / max_dd) ** 2 self.results_log.append({**params, self.metric: value, 'total_trades': n_trades}) return value except Exception as e: return float('-inf') return objective def run(self, param_space: dict) -> tuple[dict, pd.DataFrame]: sampler = TPESampler( n_startup_trials=self.n_startup, seed=42, ) study = optuna.create_study( direction=self.direction, sampler=sampler, ) study.optimize( self.create_objective(param_space), n_trials=self.n_trials, n_jobs=self.n_jobs, show_progress_bar=True, ) best_params = study.best_params results_df = pd.DataFrame(self.results_log).sort_values(self.metric, ascending=False) return best_params, results_df, study 

Example optimization run

optimizer = BayesianOptimizer( backtest_fn=lambda params: run_backtest(params, train_data), n_trials=150, n_startup_trials=15, n_jobs=4, ) param_space = { 'fast_period': {'type': 'int', 'low': 5, 'high': 30}, 'slow_period': {'type': 'int', 'low': 15, 'high': 100}, 'rsi_period': {'type': 'int', 'low': 7, 'high': 21}, 'rsi_oversold': {'type': 'int', 'low': 20, 'high': 40}, 'stop_loss_pct': {'type': 'float', 'low': 0.01, 'high': 0.10}, 'take_profit_pct': {'type': 'float', 'low': 0.02, 'high': 0.25}, 'commission': {'type': 'categorical', 'choices': ['market', 'limit']}, } best_params, results, study = optimizer.run(param_space) print("Best parameters:", best_params) print(f"Best {optimizer.metric}: {study.best_value:.3f}") 

Why TPE sampler?

TPE (Tree-structured Parzen Estimator) builds two histograms: for good and bad parameter values. The acquisition function picks points that are likely to belong to the good group. This provides robustness to noise and fast learning with few iterations. Unlike Gaussian Process, TPE scales better with 10+ parameters and does not require kernel tuning. In one project, we compared TPE with Random Search: TPE found the optimum in 80 iterations versus 600 for Random Search.

Result analysis

import optuna.visualization as vis fig = vis.plot_param_importances(study) fig.show() fig = vis.plot_contour(study, params=['fast_period', 'slow_period']) fig.show() fig = vis.plot_optimization_history(study) fig.show() 

How to protect against overfitting?

def time_series_cv_objective(params: dict, data: pd.DataFrame, n_splits: int = 5) -> dict: """K-fold cross-validation for time series""" fold_size = len(data) // (n_splits + 1) sharpe_scores = [] for fold in range(n_splits): train_start = 0 train_end = (fold + 1) * fold_size test_start = train_end test_end = train_end + fold_size train_data = data.iloc[train_start:train_end] test_data = data.iloc[test_start:test_end] metrics = run_backtest(params, test_data) sharpe_scores.append(metrics.get('sharpe_ratio', 0)) return { 'sharpe_ratio': np.mean(sharpe_scores), 'sharpe_std': np.std(sharpe_scores), 'min_sharpe': min(sharpe_scores), } 

Time series cross-validation is the only way to avoid false discoveries. We penalize for variance and minimum Sharpe: if the strategy fails in at least one period, it does not pass. Additionally, we impose penalties for too few trades (less than 15) and extreme drawdown (over 40%). In one project with crypto pairs, this cut off 70% of false patterns and saved tens of about $9–13 in savings in potential losses. Project cost starts from $5,000, offering a quick return on investment. Get a consultation on your strategy — we will evaluate it in 2 days and show how Bayesian optimization can reduce tuning time.

Comparison of optimization methods

For most tasks, Bayesian optimization with Optuna is the optimal choice: 2–3 times faster than Random Search with the same accuracy, and grid search loses by tens of times.

Work process

1. Strategy analysis — we review the code, define parameters and ranges.
2. Design — we choose the metric (Sharpe, Sortino, Profit Factor), set penalties.
3. Implementation — we write the backtest engine and integration with Optuna, implement cross-validation.
4. Testing — we run optimization on historical data (at least 3 years).
5. Deployment — we provide scripts for periodic re-optimization, documentation.

What's included in development

• Optimization module in Optuna with TPE sampler.
• Backtest engine accounting for commissions, slippage, and liquidity.
• Cross-validation with time blocks and penalties.
• Visualization of parameter importance and optimization history.
• Documentation and team training (2-hour webinar).
• 3 months of support after delivery.

Timelines and guarantees

A typical project takes from 2 to 4 weeks depending on strategy complexity. We guarantee stable code operation and assistance during deployment. Our experience — over 10 years in trading system development, 50+ completed optimization projects. Order optimization and reduce parameter tuning time by 10 times — contact us for an estimate of your project.