Files
T
Jhonathan Abreu 9eb484625d Default option pricing models for options indicators (#8229)
* Set best recommended model based on option style.

The best option model will be used in options indicators based on option style if not specified

* Reduced Greeks/IV numerical models steps to 100

* Minor fixes

* Minor changes and improvements

* Changed default pricing model for IV

* Minor performance improvements

* Minor performance improvements and fixes

* Minor change

* Minor regression algorithm fix

* Update greek indicators internal IV indicator only when not user provided

* Revert unnecessary changes

* Minor changes

* Cleanup
2024-07-30 10:15:21 -04:00

261 lines
15 KiB
C#

/*
* QUANTCONNECT.COM - Democratizing Finance, Empowering Individuals.
* Lean Algorithmic Trading Engine v2.0. Copyright 2014 QuantConnect Corporation.
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
using System;
using MathNet.Numerics.Distributions;
using Python.Runtime;
using QuantConnect.Data;
namespace QuantConnect.Indicators
{
/// <summary>
/// Option Delta indicator that calculate the delta of an option
/// </summary>
/// <remarks>sensitivity of option price relative to $1 of underlying change</remarks>
public class Delta : OptionGreeksIndicatorBase
{
/// <summary>
/// Initializes a new instance of the Delta class
/// </summary>
/// <param name="name">The name of this indicator</param>
/// <param name="option">The option to be tracked</param>
/// <param name="riskFreeRateModel">Risk-free rate model</param>
/// <param name="dividendYieldModel">Dividend yield model</param>
/// <param name="mirrorOption">The mirror option for parity calculation</param>
/// <param name="optionModel">The option pricing model used to estimate Delta</param>
/// <param name="ivModel">The option pricing model used to estimate IV</param>
public Delta(string name, Symbol option, IRiskFreeInterestRateModel riskFreeRateModel, IDividendYieldModel dividendYieldModel, Symbol mirrorOption = null,
OptionPricingModelType? optionModel = null, OptionPricingModelType? ivModel = null)
: base(name, option, riskFreeRateModel, dividendYieldModel, mirrorOption, optionModel, ivModel)
{
}
/// <summary>
/// Initializes a new instance of the Delta class
/// </summary>
/// <param name="option">The option to be tracked</param>
/// <param name="riskFreeRateModel">Risk-free rate model</param>
/// <param name="dividendYieldModel">Dividend yield model</param>
/// <param name="mirrorOption">The mirror option for parity calculation</param>
/// <param name="optionModel">The option pricing model used to estimate Delta</param>
/// <param name="ivModel">The option pricing model used to estimate IV</param>
public Delta(Symbol option, IRiskFreeInterestRateModel riskFreeRateModel, IDividendYieldModel dividendYieldModel, Symbol mirrorOption = null,
OptionPricingModelType? optionModel = null, OptionPricingModelType? ivModel = null)
: this($"Delta({option},{mirrorOption},{GetOptionModel(optionModel, option.ID.OptionStyle)})", option, riskFreeRateModel,
dividendYieldModel, mirrorOption, optionModel, ivModel)
{
}
/// <summary>
/// Initializes a new instance of the Delta class
/// </summary>
/// <param name="name">The name of this indicator</param>
/// <param name="option">The option to be tracked</param>
/// <param name="riskFreeRateModel">Risk-free rate model</param>
/// <param name="dividendYieldModel">Dividend yield model</param>
/// <param name="mirrorOption">The mirror option for parity calculation</param>
/// <param name="optionModel">The option pricing model used to estimate Delta</param>
/// <param name="ivModel">The option pricing model used to estimate IV</param>
public Delta(string name, Symbol option, PyObject riskFreeRateModel, PyObject dividendYieldModel, Symbol mirrorOption = null,
OptionPricingModelType? optionModel = null, OptionPricingModelType? ivModel = null)
: base(name, option, riskFreeRateModel, dividendYieldModel, mirrorOption, optionModel, ivModel)
{
}
/// <summary>
/// Initializes a new instance of the Delta class
/// </summary>
/// <param name="option">The option to be tracked</param>
/// <param name="riskFreeRateModel">Risk-free rate model</param>
/// <param name="dividendYieldModel">Dividend yield model</param>
/// <param name="mirrorOption">The mirror option for parity calculation</param>
/// <param name="optionModel">The option pricing model used to estimate Delta</param>
/// <param name="ivModel">The option pricing model used to estimate IV</param>
public Delta(Symbol option, PyObject riskFreeRateModel, PyObject dividendYieldModel, Symbol mirrorOption = null,
OptionPricingModelType? optionModel = null, OptionPricingModelType? ivModel = null)
: this($"Delta({option},{mirrorOption},{GetOptionModel(optionModel, option.ID.OptionStyle)})", option, riskFreeRateModel,
dividendYieldModel, mirrorOption, optionModel, ivModel)
{
}
/// <summary>
/// Initializes a new instance of the Delta class
/// </summary>
/// <param name="name">The name of this indicator</param>
/// <param name="option">The option to be tracked</param>
/// <param name="riskFreeRateModel">Risk-free rate model</param>
/// <param name="dividendYield">Dividend yield, as a constant</param>
/// <param name="mirrorOption">The mirror option for parity calculation</param>
/// <param name="optionModel">The option pricing model used to estimate Delta</param>
/// <param name="ivModel">The option pricing model used to estimate IV</param>
public Delta(string name, Symbol option, IRiskFreeInterestRateModel riskFreeRateModel, decimal dividendYield = 0.0m, Symbol mirrorOption = null,
OptionPricingModelType? optionModel = null, OptionPricingModelType? ivModel = null)
: base(name, option, riskFreeRateModel, dividendYield, mirrorOption, optionModel, ivModel)
{
}
/// <summary>
/// Initializes a new instance of the Delta class
/// </summary>
/// <param name="option">The option to be tracked</param>
/// <param name="riskFreeRateModel">Risk-free rate model</param>
/// <param name="dividendYield">Dividend yield, as a constant</param>
/// <param name="mirrorOption">The mirror option for parity calculation</param>
/// <param name="optionModel">The option pricing model used to estimate Delta</param>
/// <param name="ivModel">The option pricing model used to estimate IV</param>
public Delta(Symbol option, IRiskFreeInterestRateModel riskFreeRateModel, decimal dividendYield = 0.0m, Symbol mirrorOption = null,
OptionPricingModelType? optionModel = null, OptionPricingModelType? ivModel = null)
: this($"Delta({option},{mirrorOption},{GetOptionModel(optionModel, option.ID.OptionStyle)})", option, riskFreeRateModel,
dividendYield, mirrorOption, optionModel, ivModel)
{
}
/// <summary>
/// Initializes a new instance of the Delta class
/// </summary>
/// <param name="name">The name of this indicator</param>
/// <param name="option">The option to be tracked</param>
/// <param name="riskFreeRateModel">Risk-free rate model</param>
/// <param name="dividendYield">Dividend yield, as a constant</param>
/// <param name="mirrorOption">The mirror option for parity calculation</param>
/// <param name="optionModel">The option pricing model used to estimate Delta</param>
/// <param name="ivModel">The option pricing model used to estimate IV</param>
public Delta(string name, Symbol option, PyObject riskFreeRateModel, decimal dividendYield = 0.0m, Symbol mirrorOption = null,
OptionPricingModelType? optionModel = null, OptionPricingModelType? ivModel = null)
: base(name, option, riskFreeRateModel, dividendYield, mirrorOption, optionModel, ivModel)
{
}
/// <summary>
/// Initializes a new instance of the Delta class
/// </summary>
/// <param name="option">The option to be tracked</param>
/// <param name="riskFreeRateModel">Risk-free rate model</param>
/// <param name="dividendYield">Dividend yield, as a constant</param>
/// <param name="mirrorOption">The mirror option for parity calculation</param>
/// <param name="optionModel">The option pricing model used to estimate Delta</param>
/// <param name="ivModel">The option pricing model used to estimate IV</param>
public Delta(Symbol option, PyObject riskFreeRateModel, decimal dividendYield = 0.0m, Symbol mirrorOption = null,
OptionPricingModelType? optionModel = null, OptionPricingModelType? ivModel = null)
: this($"Delta({option},{mirrorOption},{GetOptionModel(optionModel, option.ID.OptionStyle)})", option, riskFreeRateModel,
dividendYield, mirrorOption, optionModel, ivModel)
{
}
/// <summary>
/// Initializes a new instance of the Delta class
/// </summary>
/// <param name="name">The name of this indicator</param>
/// <param name="option">The option to be tracked</param>am>
/// <param name="riskFreeRate">Risk-free rate, as a constant</param>
/// <param name="dividendYield">Dividend yield, as a constant</param>
/// <param name="mirrorOption">The mirror option for parity calculation</param>
/// <param name="optionModel">The option pricing model used to estimate Delta</param>
/// <param name="ivModel">The option pricing model used to estimate IV</param>
public Delta(string name, Symbol option, decimal riskFreeRate = 0.05m, decimal dividendYield = 0.0m, Symbol mirrorOption = null,
OptionPricingModelType? optionModel = null, OptionPricingModelType? ivModel = null)
: base(name, option, riskFreeRate, dividendYield, mirrorOption, optionModel, ivModel)
{
}
/// <summary>
/// Initializes a new instance of the Delta class
/// </summary>
/// <param name="option">The option to be tracked</param>
/// <param name="riskFreeRate">Risk-free rate, as a constant</param>
/// <param name="dividendYield">Dividend yield, as a constant</param>
/// <param name="mirrorOption">The mirror option for parity calculation</param>
/// <param name="optionModel">The option pricing model used to estimate Delta</param>
/// <param name="ivModel">The option pricing model used to estimate IV</param>
public Delta(Symbol option, decimal riskFreeRate = 0.05m, decimal dividendYield = 0.0m, Symbol mirrorOption = null,
OptionPricingModelType? optionModel = null, OptionPricingModelType? ivModel = null)
: this($"Delta({option},{mirrorOption},{GetOptionModel(optionModel, option.ID.OptionStyle)})", option, riskFreeRate, dividendYield,
mirrorOption, optionModel, ivModel)
{
}
/// <summary>
/// Calculate the Delta of the option
/// </summary>
protected override decimal CalculateGreek(decimal timeTillExpiry)
{
var math = OptionGreekIndicatorsHelper.DecimalMath;
switch (_optionModel)
{
case OptionPricingModelType.BinomialCoxRossRubinstein:
var upFactor = math(Math.Exp, ImpliedVolatility * math(Math.Sqrt, timeTillExpiry / OptionGreekIndicatorsHelper.Steps));
if (upFactor == 1)
{
// provide a small step to estimate delta
upFactor = 1.00001m;
}
var sU = UnderlyingPrice * upFactor;
var sD = UnderlyingPrice / upFactor;
var fU = OptionGreekIndicatorsHelper.CRRTheoreticalPrice(
ImpliedVolatility, sU, Strike, timeTillExpiry, RiskFreeRate, DividendYield, Right);
var fD = OptionGreekIndicatorsHelper.CRRTheoreticalPrice(
ImpliedVolatility, sD, Strike, timeTillExpiry, RiskFreeRate, DividendYield, Right);
return OptionGreekIndicatorsHelper.Divide(fU - fD, sU - sD);
case OptionPricingModelType.ForwardTree:
var discount = math(Math.Exp, (RiskFreeRate - DividendYield) * timeTillExpiry / OptionGreekIndicatorsHelper.Steps);
upFactor = math(Math.Exp, ImpliedVolatility * math(Math.Sqrt, timeTillExpiry / OptionGreekIndicatorsHelper.Steps)) * discount;
if (upFactor == 1)
{
// provide a small step to estimate delta
upFactor = 1.00001m;
}
var downFactor = math(Math.Exp, -ImpliedVolatility * math(Math.Sqrt, timeTillExpiry / OptionGreekIndicatorsHelper.Steps)) * discount;
if (downFactor == 1)
{
// provide a small step to estimate delta
downFactor = 0.99999m;
}
sU = UnderlyingPrice * upFactor;
sD = UnderlyingPrice * downFactor;
fU = OptionGreekIndicatorsHelper.ForwardTreeTheoreticalPrice(
ImpliedVolatility, sU, Strike, timeTillExpiry, RiskFreeRate, DividendYield, Right);
fD = OptionGreekIndicatorsHelper.ForwardTreeTheoreticalPrice(
ImpliedVolatility, sD, Strike, timeTillExpiry, RiskFreeRate, DividendYield, Right);
return OptionGreekIndicatorsHelper.Divide(fU - fD, sU - sD);
case OptionPricingModelType.BlackScholes:
default:
var norm = new Normal();
var d1 = OptionGreekIndicatorsHelper.CalculateD1(UnderlyingPrice, Strike, timeTillExpiry, RiskFreeRate, DividendYield, ImpliedVolatility);
decimal wholeShareDelta;
if (Right == OptionRight.Call)
{
wholeShareDelta = math(norm.CumulativeDistribution, d1);
}
else
{
wholeShareDelta = -math(norm.CumulativeDistribution, -d1);
}
return wholeShareDelta * math(Math.Exp, -DividendYield * timeTillExpiry);
}
}
}
}