/usr/include/ql/termstructures/yield/nonlinearfittingmethods.hpp is in libquantlib0-dev 1.7.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 | /* -*- mode: c++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
/*
Copyright (C) 2007 Allen Kuo
Copyright (C) 2010 Alessandro Roveda
This file is part of QuantLib, a free-software/open-source library
for financial quantitative analysts and developers - http://quantlib.org/
QuantLib is free software: you can redistribute it and/or modify it
under the terms of the QuantLib license. You should have received a
copy of the license along with this program; if not, please email
<quantlib-dev@lists.sf.net>. The license is also available online at
<http://quantlib.org/license.shtml>.
This program is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
FOR A PARTICULAR PURPOSE. See the license for more details.
*/
/*! \file nonlinearfittingmethods.hpp
\brief nonlinear methods to fit a bond discount function
*/
#ifndef quantlib_nonlinear_fitting_methods_hpp
#define quantlib_nonlinear_fitting_methods_hpp
#include <ql/termstructures/yield/fittedbonddiscountcurve.hpp>
#include <ql/math/bspline.hpp>
#include <boost/shared_ptr.hpp>
namespace QuantLib {
//! Exponential-splines fitting method
/*! Fits a discount function to the exponential form
\f[
d(t) = \sum_{i=1}^9 c_i \exp^{-kappa i t}
\f]
where the constants \f$ c_i \f$ and \f$ \kappa \f$ are to be
determined. See:Li, B., E. DeWetering, G. Lucas, R. Brenner
and A. Shapiro (2001): "Merrill Lynch Exponential Spline
Model." Merrill Lynch Working Paper
\warning convergence may be slow
*/
class ExponentialSplinesFitting
: public FittedBondDiscountCurve::FittingMethod {
public:
ExponentialSplinesFitting(bool constrainAtZero = true,
const Array& weights = Array(),
boost::shared_ptr<OptimizationMethod> optimizationMethod
= boost::shared_ptr<OptimizationMethod>());
std::auto_ptr<FittedBondDiscountCurve::FittingMethod> clone() const;
private:
Size size() const;
DiscountFactor discountFunction(const Array& x, Time t) const;
};
//! Nelson-Siegel fitting method
/*! Fits a discount function to the form
\f$ d(t) = \exp^{-r t}, \f$ where the zero rate \f$r\f$ is defined as
\f[
r \equiv c_0 + (c_0 + c_1)*(1 - exp^{-\kappa*t}/(\kappa t) -
c_2 exp^{ - \kappa t}.
\f]
See: Nelson, C. and A. Siegel (1985): "Parsimonious modeling of yield
curves for US Treasury bills." NBER Working Paper Series, no 1594.
*/
class NelsonSiegelFitting
: public FittedBondDiscountCurve::FittingMethod {
public:
NelsonSiegelFitting(const Array& weights = Array(),
boost::shared_ptr<OptimizationMethod> optimizationMethod
= boost::shared_ptr<OptimizationMethod>());
std::auto_ptr<FittedBondDiscountCurve::FittingMethod> clone() const;
private:
Size size() const;
DiscountFactor discountFunction(const Array& x, Time t) const;
};
//! Svensson Fitting method
/*! Fits a discount function to the form
\f$ d(t) = \exp^{-r t}, \f$ where the zero rate \f$r\f$ is defined as
\f[
r \equiv c_0 + (c_0 + c_1)(\frac {1 - exp^{-\kappa t}}{\kappa t})
- c_2exp^{ - \kappa t}
+ c_3{(\frac{1 - exp^{-\kappa_1 t}}{\kappa_1 t} -exp^{-\kappa_1 t})}.
\f]
See: Svensson, L. (1994). Estimating and interpreting forward
interest rates: Sweden 1992-4.
Discussion paper, Centre for Economic Policy Research(1051).
*/
class SvenssonFitting
: public FittedBondDiscountCurve::FittingMethod {
public:
SvenssonFitting(const Array& weights = Array(),
boost::shared_ptr<OptimizationMethod> optimizationMethod
= boost::shared_ptr<OptimizationMethod>());
std::auto_ptr<FittedBondDiscountCurve::FittingMethod> clone() const;
private:
Size size() const;
DiscountFactor discountFunction(const Array& x, Time t) const;
};
//! CubicSpline B-splines fitting method
/*! Fits a discount function to a set of cubic B-splines
\f$ N_{i,3}(t) \f$, i.e.,
\f[
d(t) = \sum_{i=0}^{n} c_i * N_{i,3}(t)
\f]
See: McCulloch, J. 1971, "Measuring the Term Structure of
Interest Rates." Journal of Business, 44: 19-31
McCulloch, J. 1975, "The tax adjusted yield curve."
Journal of Finance, XXX811-30
\warning "The results are extremely sensitive to the number
and location of the knot points, and there is no
optimal way of selecting them." James, J. and
N. Webber, "Interest Rate Modelling" John Wiley,
2000, pp. 440.
*/
class CubicBSplinesFitting
: public FittedBondDiscountCurve::FittingMethod {
public:
CubicBSplinesFitting(const std::vector<Time>& knotVector,
bool constrainAtZero = true,
const Array& weights = Array(),
boost::shared_ptr<OptimizationMethod> optimizationMethod
= boost::shared_ptr<OptimizationMethod>());
//! cubic B-spline basis functions
Real basisFunction(Integer i, Time t) const;
std::auto_ptr<FittedBondDiscountCurve::FittingMethod> clone() const;
private:
Size size() const;
DiscountFactor discountFunction(const Array& x, Time t) const;
BSpline splines_;
Size size_;
//! N_th basis function coefficient to solve for when d(0)=1
Natural N_;
};
//! Simple polynomial fitting method
/* Fits a discount function to the simple polynomial form:
\f[
d(t) = \sum_{i=0}^{degree} c_i * t^{i}
\f]
where the constants \f$ c_i \f$ are to be determined.
This is a simple/crude, but fast and robust, means of fitting
a yield curve.
*/
class SimplePolynomialFitting
: public FittedBondDiscountCurve::FittingMethod {
public:
SimplePolynomialFitting(Natural degree,
bool constrainAtZero = true,
const Array& weights = Array(),
boost::shared_ptr<OptimizationMethod> optimizationMethod
= boost::shared_ptr<OptimizationMethod>());
std::auto_ptr<FittedBondDiscountCurve::FittingMethod> clone() const;
private:
Size size() const;
DiscountFactor discountFunction(const Array& x, Time t) const;
Size size_;
};
//! Spread fitting method helper
/* Fits a spread curve on top of a discount function according to given parametric method
*/
class SpreadFittingMethod
: public FittedBondDiscountCurve::FittingMethod {
public:
SpreadFittingMethod(boost::shared_ptr<FittingMethod> method,
Handle<YieldTermStructure> discountCurve);
std::auto_ptr<FittedBondDiscountCurve::FittingMethod> clone() const;
protected:
void init();
private:
Size size() const;
DiscountFactor discountFunction(const Array& x, Time t) const;
// underlying parametric method
boost::shared_ptr<FittingMethod> method_;
// adjustment in case underlying discount curve has different reference date
DiscountFactor rebase_;
// discount curve from on top of which the spread will be calculated
Handle<YieldTermStructure> discountingCurve_;
};
}
#endif
|