Next Article in Journal
Robust Classification via Support Vector Machines
Next Article in Special Issue
Special Issue “Actuarial Mathematics and Risk Management”
Previous Article in Journal
Multiple Bonus–Malus Scale Models for Insureds of Different Sizes
Previous Article in Special Issue
Equivalent Risk Indicators: VaR, TCE, and Beyond
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparison of Macaulay Approximations

by
Stefanos C. Orfanos
Department of Risk Management & Insurance, J. Mack Robinson College of Business, Georgia State University, Atlanta, GA 30302, USA
Risks 2022, 10(8), 153; https://doi.org/10.3390/risks10080153
Submission received: 13 June 2022 / Revised: 8 July 2022 / Accepted: 15 July 2022 / Published: 29 July 2022
(This article belongs to the Special Issue Actuarial Mathematics and Risk Management)

Abstract

:
We discuss several known formulas that use the Macaulay duration and convexity of commonly used cash flow streams to approximate their net present value, and compare them with a new approximation formula that involves hyperbolic functions. Our objective is to assess the reliability of each approximation formula under different scenarios. The results in this note should be of interest to actuarial candidates and educators as well as analysts working in all areas of actuarial practice.

1. Introduction

Actuaries and actuarial science students at universities all over the world are familiar with approximation formulas for the present value of cash flow streams using some notion of cash flow duration or convexity. For example, the syllabus of Exam FM of the US-based Society of Actuaries includes the topic of approximations using the Macaulay and modified duration and convexity, while the UK-based Institute and Faculty of Actuaries in its material for exam CM1 mentions approximations derived from a Taylor series expansion.
Beside the academic and pedagogical interest in such approximation formulas, one may also consider the practical value in the management of interest rate risk. Although abundant computing power has enabled firms to implement elaborate immunization strategies that incorporate multi-factor stochastic interest rate models, non-parallel yield curve shifts, and complicated asset and liability characteristics, the restrictions posed by a simplistic valuation model are not unreasonable if rates remain historically low, yield curves stay relatively flat, and we can control the potential errors. Indeed, it may be helpful to know which approximation formula proves to be the most reliable, and to use it as a quick validation tool when time constraints preclude the use of a more sophisticated approach.
Alps (2017) describes a realistic scenario involving an investment actuary and her CEO, where the use of an approximation formula would be warranted or even necessitated. This is especially true in today’s world of fast-changing rates, when companies have to react almost instantly to benchmark fund rates and quantitative tightening decisions by the Federal Reserve or other central banks.
In this note, we discuss several known formulas that use the Macaulay duration and convexity of commonly used cash flow streams to approximate their net present value, and compare them with a new approximation formula that involves hyperbolic functions. In addition to annuities, dividend stocks, and bonds, we also consider the cases of negative payments and embedded options to perform a deeper assessment. The notions of effective duration and convexity are defined in the next section and used to price the embedded options. Our objective is to measure the reliability of each approximation formula under different scenarios. As alluded to earlier, we only consider parallel interest rate shocks to flat yield curves.

1.1. Literature Review

The idea of using a bond’s duration to approximate changes to its price goes back to Macaulay (1938). Some authors credit Fischer and Weil (1971) with the publication of the first duration–convexity approximation formula. Enhancements of that formula by controlling the missing higher-order terms and incorporating passage of time were announced in Jarjir and Rakotondratsimba (2008, 2012), though the resulting formulas contain parameters that are unintuitive and hard to calibrate. A conceptually simpler formula was given in Barber (1995), and independently in Livingston and Zhou (2005) for the modified duration, which was subsequently generalized to a duration–convexity model in Tchuindjo (2008). Further work in Barber and Dandapani (2017) considered negative-yielding bonds, and Johansson (2012) added passage of time. A separate duration–convexity formula appears in Alps (2017) and is applied to an empirical study of basic immunization strategies in Nie et al. (2021), while a very recent paper by Barber (2022) further generalizes a duration–convexity approximation by introducing an additional ‘compounding’ parameter. Finally, traditional approximations have been implemented in statistical analysis packages; see Lee (2021) for R code.

1.2. Notation

We denote the net present value of a cash flow stream by P. The interest rate r is annualized and continuously compounded (i.e., force of interest). Δ r is the change in interest rates from the initial value r 0 to r. Finally, the annual discount factor v is by definition equal to e r .
Throughout the remainder of this paper and for convenience, assume r 0 = 1.6 % , which is approximately the yield on the 10-year T-bond at the beginning of this year.

2. Materials and Methods

Recall that the Macaulay duration of a stream of cash flows { C F t j } j = 1 n being paid at future times { t j } j = 1 n is defined by
d = j = 1 n C F t j v t j t j j = 1 n C F t j v t j = d P / d r P ,
while its Macaulay convexity is given by
c = j = 1 n C F t j v t j t j 2 j = 1 n C F t j v t j = d 2 P / d r 2 P .
We do not consider the modified duration here because the Macaulay duration has a more intuitive interpretation (being the ‘average’ timing of the cash flows) and tends to result in tighter approximations for non-negative rates.
In the case of bonds with embedded options, it will be necessary to price the value of the option using a simple Black model. Recall that the pricing formula for, say, a European call option with expiration at time t and strike K is
V = v t P Φ ( d 1 ) K Φ ( d 2 )
where Φ represents the standard normal c.d.f. and the quantities d 1 , 2 are given by
d 1 , 2 = ln ( P / K ) σ t ± σ t 2
with σ the bond price volatility. We also need more flexible measures of bond duration and convexity. To that end, define the effective duration by means of
d e = P ( r 0 + Δ r ) P ( r 0 Δ r ) 2 P 0 Δ r
and the effective convexity as
c e = P ( r 0 + Δ r ) 2 P 0 + P ( r 0 Δ r ) P 0 ( Δ r ) 2 .

2.1. Fischer–Weil’s Approximation

This follows immediately from Calculus and the definitions above.
Δ P P 0 d 0 Δ r + c 0 2 ( Δ r ) 2 .
It is assumed that the Macaulay duration and convexity are computed at rate r 0 , hence the subscripts.

2.2. Barber’ 1995 Approximation

Instead of the second-order Taylor polynomial of P, we consider the first-order Taylor polynomial of ln P , thus obtaining
ln P ln P 0 d 0 Δ r ,
thus
P P 0 e d 0 Δ r .
Unlike the first-order Taylor polynomial in P that has no convexity, the functional form of Barber’s approximation bequeaths it with a certain degree of positive curvature. This leads to good approximation results whenever c 0 d 0 2 and poor performance for c 0 < 0 .

2.3. Tchuindjo’ Approximation

Similar to Barber’s approximation, but involving the second-order Taylor polynomial of ln P
ln P ln P 0 d 0 Δ r + c 0 d 0 2 2 ( Δ r ) 2 ,
from which one solves for P. The added quadratic term gives better results in cases where c 0 d 0 2 is non-trivial, but may still introduce large errors whenever c 0 < 0 .

2.4. Alps’ Approximation

The approximation formula and its derivation can be found in Alps (2017). The central idea in the derivation of this approximation is to compute a Taylor polynomial for the current value of the cash flow stream at time t = d 0 . This choice results in high accuracy in situations where d 0 > 0 , and less so for d 0 < 0 .
We have rewritten it below in terms of continuously compounded interest rates.
P P 0 e d 0 Δ r 1 + c 0 d 0 2 2 e Δ r 1 2 .

2.5. Hyperbolic Approximation

We have not encountered this approximation formula in the literature and we assume its derivation is presented here for the first time. Consider the homogeneous differential equation P c 0 P = 0 that mimics the definition of Macauley convexity given earlier in this note. Its general solution takes the form
P = a e c 0 r + b e c 0 r
(do not worry for the time being about the case c 0 < 0 .) Setting P ( r 0 ) = P 0 and P ( r 0 ) = d 0 P 0 , which is a reformulation of the definition of Macauley duration, one obtains the approximation
P P 0 1 2 1 d 0 c 0 e c 0 Δ r + 1 2 1 + d 0 c 0 e c 0 Δ r
which can be rewritten as
P P 0 cosh c 0 Δ r d 0 c 0 sinh c 0 Δ r .
The well-known trig identities
cosh ( i θ ) = e i θ + e i θ 2 = cos θ , sinh ( i θ ) = e i θ e i θ 2 = i sin θ
can be used in the case c 0 < 0 to obtain
P P 0 cos | c 0 | Δ r d 0 | c 0 | sin | c 0 | Δ r ,
which is useful whenever there is a computational issue with imaginary numbers.
In the next section, we demonstrate that the hyperbolic approximation is less prone to errors than other well-known approximations in situations where the duration and/or convexity are negative. Recall that negative convexity cash flow streams can be easily constructed with the addition of negative cash flows to a stream of positive payments, when considering callable bonds, or with mortgage-backed securities due to the prepayment option in conventional residential mortgages.

3. Results

Approximation formulas such as Equations (1)–(5) should ideally be intuitive and behave well in special cases.
(i)
The simplest cash flow is cash, which has trivial duration and convexity and is unaffected by interest rate changes. By substituting d 0 = c 0 = 0 or taking the corresponding limit in the case of (5) and using the fact that
lim θ 0 sinh θ θ = lim θ 0 sin θ θ = 1 ,
we obtain P = P 0 as expected.
(ii)
Next, take a zero-coupon bond, for which c 0 = d 0 2 : Except for Fischer–Weil’s approximation, the rest reduce to Barber’s approximation, which is perfectly accurate in this case. On the other hand, the error associated with Fischer-Weil’s approximation increases with the bond duration and it can be as high as 0.56% for a 30-year zero-coupon bond after a 100 bp increase in rates.
(iii)
For a convexity-hedged ( c 0 0 ) portfolio, Fischer–Weil’s and the hyperbolic approximations reduce to the first-order approximation Δ P d 0 P 0 Δ r . The corresponding results for the other approximation formulas are not as intuitive and their accuracy relative to the above approximation cannot be determined without additional details about the cash flow characteristics.
We supplement the theoretical tests above with some concrete examples.
(iv)
Consider a 10-year annuity-immediate with annual payments of 10. Recall that our assumption is r 0 = 1.6 % and compute the present value P 0 = 10 a 10 = 91.6728 . Another easy calculation gives the Macaulay duration and convexity as d 0 = 5.3681 and c 0 = 37.0554 , respectively.
In Table 1, the exact value of P is computed using the same formula as for P 0 but at the new continuously compounded rate.
We can observe that Tchuindjo’s approximation outperforms the rest, while Barber’s lags behind for sizable rate changes. This was to be expected, since Barber’s approximation lacks a convexity term and will not do well in cases when c 0 d 0 2 is non-trivial. On the other hand, Alps’ and the hyperbolic approximations are roughly equally accurate behind Tchuindjo’s.
(v)
Next, add a negative cash flow at time 20. We have chosen C F 20 = 120 in the example below; the net present value is P 0 = 4.5349 and the Macaulay duration and convexity are d 0 = 275.7817 and c 0 = 6936.8498 , respectively.
Looking at Table 2 below, it may come as a surprise that the approximations by Tchuindjo and Alps blow up completely. However, we can provide a simple mathematical explanation for the bizarre behavior. Whenever c 0 < 0 , the quadratic term of these two approximations that includes the expression c 0 d 0 2 has the potential to be extremely influential. As Δ r increases, said term can overwhelm the baseline value P 0 and the linear term, resulting in large errors. Barber’s approximation exhibits the opposite weakness: missing a quadratic term implies that the negative convexity is not accounted for at all. In fact, for suitable C F 20 , we can obtain d 0 = 0 , in which case Barber’s approximation fails to yield any results.
We conclude this example by mentioning that the top-performing approximation is Fischer–Weil’s, while the hyperbolic approximation is second-best.
(vi)
Let us now consider a dividend stock, whose theoretical price is computed using Gordon’s dividend discount model
P = D r g
with D representing next year’s dividend and g its constant continuously compounded growth rate in perpetuity. A quick calculation gives d = ( r g ) 1 and c = 2 ( r g ) 2 ; for g = 0.6 % we obtain d 0 = 100 and c 0 = 20 , 000 . Assume D = 1 .
Some of the results in Table 3 may appear counterintuitive at first sight. Gordon’s model suggests that P has an inverse relationship to r; however, all approximations except for Alps’ and Barber’s eventually produce a divergent estimate for P as Δ r increases. However, this is explained by the fact that we are attempting to trace a hyperbola using quadratic curves. Moreover, all approximations struggle to keep up with P for large negative values of Δ r .
Overall, Alps’ approximation proves to be the most dependable for moderate changes in the interest rates.
(vii) 
Next, consider a 10-year bond with a coupon rate of r 0 and face value of 100. A quick calculation yields d 0 = 9.3151 and c 0 = 90.6932 .
It turns out that the last three approximation formulas clearly outperform the rest, with Tchuindjo’s having a slight advantage over Alps’ and the hyperbolic approximation, as evidenced from Table 4. The subpar performance of Fischer–Weil on bonds is one of the reasons why this approximation is not widely utilized, despite its robustness in cases such as (v).
It has been shown empirically that although investment-grade bonds fall in price when interest rates rise, that is not necessarily the case with high-yield bonds whose duration can be negative due to default risk; see Melentyev and Yu (2020). For such bonds, care should be exercised when using the approximations by Tchuindjo or Alps.
(viii)
Finally, assume the bond is callable, with the European call strike set at K = 101.0000 and bond price volatility σ = 8 % . The call is exercised a year ahead of the bond’s maturity and has price V = 9.9431 , which is subtracted from the price of a conventional bond to arrive at the callable bond price. Using Δ r = 20 bp in the calculation of the effective duration and convexity, we obtain d e = 5.3333 and c e = 50.3993 . The positive convexity may surprise some readers, but note that the convexity turns negative when the interest rate gets closer to 0 and the bond price approaches the strike.
It is important to observe in Table 5 that none of the approximation formulas can consistently outperform the rest, if our objective is to estimate the full range of prices for such a bond. In effect, we are trying to approximate a function with an inflection point using quadratic curves, and thus significant approximation errors are inevitable.
The only useful conclusion is that the hyperbolic approximation is never the worst one, since it tends to be “sandwiched” between other approximations.

4. Discussion

We have established through a number of theoretical considerations and concrete examples that the accuracy of various Macaulay approximations can vary widely. Approximations that outperform in one case turn out to be unreliable in another case. The hyperbolic approximation, introduced in this paper, exhibited modest errors in most cases and thus the most reliability among the five approximations studied.
We can envision a variety of uses for the results presented here:
  • To perform expeditious interest risk calculations by practitioners;
  • As a study note to gain insight into risk management concepts that are tested in the actuarial examinations in the US and Europe;
  • As potential areas of student research or as assigned projects that utilize real financial data in actuarial science classes taught by academics.
There is also potential to expand the scope of this study by incorporating non-flat yield curves, key rate durations, passage of time, and more complex financial instruments.

Funding

This research received no external funding.

Acknowledgments

The author wishes to thank the anonymous referees for their constructive comments. Partial results were presented at the 56th Actuarial Research Conference (see Orfanos 2021) and this work benefitted from the attendees’ questions.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Alps, Robert. 2017. Using Duration and Convexity to Approximate Change in Present Value. SOA Financial Mathematics Study Note FM-24-17. pp. 1–18. Available online: https://www.soa.org/globalassets/assets/Files/Edu/2017/fm-duration-convexity-present-value.pdf (accessed on 22 July 2022).
  2. Barber, Joel R. 1995. A Note on Approximating Bond Price Sensitivity Using Duration and Convexity. The Journal of Fixed Income 4: 95–98. [Google Scholar] [CrossRef]
  3. Barber, Joel R. 2022. Power Law Bond Price and Yield Approximation. The Journal of Risk Finance 23: 14–31. [Google Scholar] [CrossRef]
  4. Barber, Joel R., and Krishnan Dandapani. 2017. Interest Rate Risk in a Negative Yielding World. Frontiers in Finance and Economics 14: 1–19. [Google Scholar]
  5. Fischer, Lawrence, and Roman L. Weil. 1971. Coping with the Risk of Interest-Rate Fluctuations: Returns to Bondholders from Naive and Optimal Strategies. The Journal of Business 44: 406–31. [Google Scholar] [CrossRef]
  6. Livingston, Miles, and Lei Zhou. 2005. Exponential Duration: A More Accurate Estimate of Interest Rate Risk. The Journal of Financial Research 28: 343–61. [Google Scholar] [CrossRef]
  7. Jarjir, Souad L., and Yves Rakotondratsimba. 2008. Revisiting the Bond Duration-Convexity Approximation. SSRN Electronic Journal. [Google Scholar] [CrossRef]
  8. Jarjir, Souad L., and Yves Rakotondratsimba. 2012. Enhancement of the Bond Duration-Convexity Approximation. International Journal of Economics and Finance 4: 115–25. [Google Scholar]
  9. Johansson, Bo. 2012. A Note on Approximating Bond Returns Allowing for Both Yield Change and Time Passage. Munich Personal RePEc Archive Paper No. 92607. Available online: https://mpra.ub.uni-muenchen.de/92607/1/MPRA_paper_92607.pdf (accessed on 22 July 2022).
  10. Lee, Sang-Heon. 2021. Percent Change of Bond Price Using Duration and Convexity in R. R bloggers. Available online: https://www.r-bloggers.com/2021/09/percent-change-of-bond-price-using-duration-and-convexity-in-r/ (accessed on 22 July 2022).
  11. Macaulay, Fblackerick R. 1938. Some Theoretical Problems Suggested by Movement of Interest Rates, Bonds, Yields, and Stock Prices in the United States Since 1856. New York: National Bureau of Economic Research, pp. 45–53. [Google Scholar]
  12. Melentyev, Oleg, and Eric Yu. 2020. Bank of America/Merrill Lynch High Yield Strategy. New York: Bank of America Merrill Lynch, September 4. [Google Scholar]
  13. Nie, Junming, Zhihao Wu, Shiyao Wang, and Ye Chen. 2021. Duration Strategy for Bond Investment Based on an Empirical Study. Advances in Economics, Business and Management Research 203: 1726–30. [Google Scholar]
  14. Orfanos, Stefanos. 2021. A Comparison of Macaulay approximations. Research talk at the 56th Actuarial Research Conference. Available online: https://csh.depaul.edu/academics/mathematical-sciences/about/arc2021/Documents/Parallel-Session-Abstracts.pdf (accessed on 22 July 2022).
  15. Tchuindjo, Leonard. 2008. An Accurate Formula for Bond-Portfolio Stress Testing. The Journal of Risk Finance 9: 262–77. [Google Scholar] [CrossRef]
Table 1. PV of 10-year annuity with annual payments of 10.
Table 1. PV of 10-year annuity with annual payments of 10.
Δ r ExactFischer-WeilBarberTchuindjoAlpsHyperbolic
100 bp96.768296.763796.728396.768296.766996.7668
80 bp95.720795.718495.695495.720795.719995.7199
60 bp94.687694.686694.673594.687694.687194.6873
40 bp93.668793.668493.662593.668793.668593.6686
20 bp92.663992.663892.662392.663992.663892.6639
0 bp91.672891.672891.672891.672891.672891.6728
20 bp90.695490.695490.693990.695490.695390.6954
40 bp89.731389.731689.725489.731389.731189.7314
60 bp88.780488.781388.767288.780488.780088.7807
80 bp87.842587.844787.819387.842587.841787.8432
100 bp86.917386.921686.881586.917386.916286.9186
Table 2. NPV of 10-year annuity with annual payments of 10 and a payment of 120 at time 20.
Table 2. NPV of 10-year annuity with annual payments of 10 and a payment of 120 at time 20.
Δ r ExactFischer-WeilBarberTchuindjoAlpsHyperbolic
100 bp 9.6622 9.5445 0.28770.0045 0.8684 8.0590
80 bp 6.5366 6.4769 0.49940.0351 0.7850 5.7162
60 bp 3.5601 3.5352 0.86690.1946 0.3873 3.2151
40 bp 0.7266 0.7193 1.50480.77470.5372 0.6250
20 bp1.96981.97072.61232.21282.19241.9824
0 bp4.53494.53494.53494.53494.53494.5349
20 bp6.97426.97337.87256.66856.60696.9619
40 bp9.29299.285913.66647.03574.87859.1962
60 bp11.496011.472623.72435.3262 10.6004 11.1758
80 bp13.588513.533541.18462.8930 64.7471 12.8461
100 bp15.574815.468671.49501.1275 215.8388 14.1608
Table 3. Price of a dividend stock with D = 1 and g = 0.6 % .
Table 3. Price of a dividend stock with D = 1 and g = 0.6 % .
Δ r ExactFischer-WeilBarberTchuindjoAlpsHyperbolic
100 bpn/a300.0000271.8282448.1689403.4619354.6482
80 bp500.0000244.0000222.5541306.4854291.5285269.3175
60 bp250.0000196.0000182.2119218.1472213.9771205.6762
40 bp166.6667156.0000149.1825161.6074160.7412158.5990
20 bp125.0000124.0000122.1403101.8813101.8813101.8813
0 bp100.0000100.0000100.0000100.0000100.0000100.0000
20 bp83.333384.000081.873183.527083.459083.7590
40 bp71.428676.000067.032072.614972.225774.2635
60 bp62.500076.000054.881265.704764.448770.7488
80 bp55.555684.000044.932961.878358.858672.9319
100 bp50.0000100.000036.787960.653154.602680.9885
Table 4. PV of 10-year par bond with coupon rate r 0 .
Table 4. PV of 10-year par bond with coupon rate r 0 .
Δ r ExactFischer-WeilBarberTchuindjoAlpsHyperbolic
100 bp109.7839109.7686109.7628109.7843109.7836109.7830
80 bp107.7501107.7423107.7368107.7503107.7499107.7497
60 bp105.7556105.7523105.7482105.7557105.7555105.7554
40 bp103.7996103.7986103.7963103.7996103.7995103.7995
20 bp101.8813101.8812101.8805101.8813101.8813101.8813
0 bp100.0000100.0000100.0000100.0000100.0000100.0000
20 bp98.155098.155198.154298.155098.155098.1550
40 bp96.345696.346596.342596.345596.345496.3456
60 bp94.571094.574294.564294.570994.570794.5712
80 bp92.830692.838192.818892.830492.830192.8310
100 bp91.123891.138391.105691.123491.122991.1246
Table 5. NPV of 10-year par bond with coupon rate r 0 , callable for 101 at t = 9 .
Table 5. NPV of 10-year par bond with coupon rate r 0 , callable for 101 at t = 9 .
Δ r ExactFischer-WeilBarberTchuindjoAlpsHyperbolic
100 bp95.059495.086894.990395.094695.091395.0910
80 bp94.036794.044593.982494.048594.046494.0466
60 bp93.019793.020492.985393.022093.020993.0213
40 bp92.014892.014491.998792.014992.014492.0146
20 bp91.026591.026591.022691.026691.026591.0266
0 bp90.056990.056990.056990.056990.056990.0569
20 bp89.105389.105389.101489.105389.105289.1053
40 bp88.169288.172088.156088.171588.171088.1717
60 bp87.243787.256887.220787.255287.254187.2559
80 bp86.322786.359786.295386.355986.354086.3577
100 bp85.399185.480885.379785.473585.470585.4768
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Orfanos, S.C. A Comparison of Macaulay Approximations. Risks 2022, 10, 153. https://doi.org/10.3390/risks10080153

AMA Style

Orfanos SC. A Comparison of Macaulay Approximations. Risks. 2022; 10(8):153. https://doi.org/10.3390/risks10080153

Chicago/Turabian Style

Orfanos, Stefanos C. 2022. "A Comparison of Macaulay Approximations" Risks 10, no. 8: 153. https://doi.org/10.3390/risks10080153

APA Style

Orfanos, S. C. (2022). A Comparison of Macaulay Approximations. Risks, 10(8), 153. https://doi.org/10.3390/risks10080153

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop