{"id":208425,"date":"2024-06-28T10:23:01","date_gmt":"2024-06-28T14:23:01","guid":{"rendered":"https:\/\/ibkrcampus.com\/?p=208425"},"modified":"2024-08-12T15:28:50","modified_gmt":"2024-08-12T19:28:50","slug":"coupon-bearing-bond-pricing-using-r-code","status":"publish","type":"post","link":"https:\/\/www.interactivebrokers.com\/campus\/ibkr-quant-news\/coupon-bearing-bond-pricing-using-r-code\/","title":{"rendered":"Coupon Bearing Bond Pricing Using R Code"},"content":{"rendered":"\n
This post explains how to calculate the price of some complicated coupon-bearing bonds using R code.<\/p>\n\n\n\n
Introduction<\/h3>\n\n\n\n
There are pricing formula for coupon bond as well as discount bond which are used among practioners under the market convention.<\/p>\n\n\n\n
Bond Pricing Formula<\/h3>\n\n\n\n<\/figure>\n\n\n\n<\/figure>\n\n\n\n
From the above formula, denominator \nindicates that if valuation date tval<\/sub>\nfall on between interest payment dates \nthe simple interest rate calculation is applied to the period tval<\/sub> \u2212 tni<\/sub>, not continously compounding.<\/p>\n\n\n\n
Bond Pricing Example<\/h3>\n\n\n\n
The following table shows the cash flow schedule of coupon bond with in-advance<\/strong> interest payments. From the fact that the number of cash flow starts from 8 to 14, we can find that the recent settled cash flow payment is 6th not 7th. This contract is of in-advance interest payment features. See the accrual begin date and payment date.<\/p>\n\n\n\n<\/figure>\n\n\n\n
Consider two cases :<\/p>\n\n\n\n
\n
valuation date = 7th interest payment date<\/li>\n\n\n\n
7th interet payment date < valuation date < 8th interest payment date<\/li>\n<\/ol>\n\n\n\n
The maturity date is 2019-07-29 which is 4 more days later than 2019-07-25 (last day in a quarterly basis). Taking this into account, we can draw cash flows time lines for two cases.<\/p>\n\n\n\n<\/figure>\n\n\n\n
In other words, since maturity date is 2019-07-29, interet for 4 remaining days is paid at 2019-07-25 in advance. Of couse, principal redemption is made at 2019-07-29.<\/p>\n\n\n\n
But for calculation is done for the second case since the second case is more general than the first and the first case is simple to calculate. With these in-advance payment features, the pricing equation for this coupon bond (clean price) at 2018-03-22 is as follows.<\/p>\n\n\n\n<\/figure>\n\n\n\n
The explanation of pricing process with cash flow timeline and equation above is<\/p>\n\n\n\n
1) Amount in [ ] is the discounted value of all cashflows since 2018\/04\/26 evaluating at 2018\/04\/26 (first payment date) not 2018\/03\/22 (valuation date).<\/p>\n\n\n\n
2) Discount amount in [] with 35\/90 quarter simply compouding to 2018\/03\/22 value (dirty price).<\/p>\n\n\n\n
3) The portion of 35\/90 in interest previously paid at 2018\/01\/26 should be returned to buyer. This is called as accrued interest and reflected to the last term.<\/p>\n\n\n\n
4) Therefore, dirty price + accrued interest is the clean price which is equal to P????.<\/p>\n\n\n\n
It is a market convention that dirty and clean price of in-advance interest payments is the other way around of regular coupon bond with in-arrears interest payment.<\/p>\n\n\n\n
R code for Bond Pricing<\/h3>\n\n\n\n
The following R code is implemented using the above formula and specific considerations such as in-advance features, in-between settlement, … etc.<\/p>\n\n\n\n
#=========================================================================#\n# Financial Econometrics & Derivatives, ML\/DL using R, Python, Tensorflow \n# by Sang-Heon Lee \n#\n# https:\/\/shleemodel.blogspot.com \n#-------------------------------------------------------------------------#\n# Coupon Bond Valuation with in advance interest payment\n#=========================================================================#\n \ngraphics.off() # clear all graphs\nrm(list=ls()) # clear all datasets\n \n#-----------------------------------------------------------------#\n# Specification\n#-----------------------------------------------------------------#\n# valuation date : 2018-03-22\n# maturity date : 2019-07-29\n# coupon rate : annual 6.2%, quarterly coupon with in advance\n# principal amount : Bullet repayment at maturity \n#-----------------------------------------------------------------#\n \n cpn.rate <- 0.062 # coupon rate\n int.rate <- 0.062 # interest rate\n mat.qrt <- 6 # # of interest payments\n \n # initialization\n cpn.amt <- 0 # coupon payments\n prn.amt <- 0 # principal amount\n price <- 0 # bond price\n \n # 2018.04.26's discounted value of coupon payments\n for (t in 0:(mat.qrt-2)) { \n cpn.amt <- cpn.amt + (cpn.rate\/4)*10000\/\n (1+int.rate\/4)^t\n }\n \n # 2018.04.26's discounted value of \n # 4-day coupon payment at 2019.07.25\n cpn.amt <- cpn.amt + (cpn.rate\/4)*(4\/93)*10000\/\n (1+int.rate\/4)^(mat.qrt-1)\n \n # 2018.04.26's discounted value of \n #principal amount at 2019.07.29\n prn.amt <- 10000\/\n ((1+(int.rate\/4)*(4\/93))*(1+int.rate\/4)^(mat.qrt-1))\n \n # 2018.04.26's bond price\n price <- cpn.amt + prn.amt\n \n print(paste(\" bond price before discounting at pricing date:\", price))\n \n #-----------------------------------------------------------------#\n # 2018.03.22's discounted value of all cash flow\n # = present value of all cash flow\n # = bond price at valuation date\n #-----------------------------------------------------------------#\n # present value discounted at pricing date using simple interest \n #-----------------------------------------------------------------#\n price = price\/(1+(int.rate\/4)*(35\/90))\n \n print(paste(\"bond price discounting at pricing date using simple interest :\", price))\n \n dirty.price <- round(price,2)\n accrued.interest <- round((10000*int.rate\/4)*(35\/90),2)\n clean.price <- round(dirty.price + accrued.interest,2)\n \n print(\"clean price = dirty price + accrued interest\") \n print(paste(clean.price, \"=\", price, \"+\", accrued.interest))<\/pre>\n\n\n\n
The result is as follows.<\/p>\n\n\n\n
[1] \"price before discounting at pricing date: 10011.47\"\n[1] \"price discounting at pricing date using simple interest: 9951.49\"\n[1] \"clean price = dirty price + accrued interest\"\n[1] \"10011.77 = 9951.4922404715 + 60.28\"<\/pre>\n\n\n\n
This post has showed the pricing of some complicated coupon bond, not simple regular coupon bond. But the bottom line is to apply simple interest to incomplete accrual period and compound interest to complete accrual period.<\/p>\n\n\n\n
Next time, let’s dive into YTM(Yield to Maturity) and reinvestment risk with EXCEL illustrations and R code.<\/p>\n\n\n\n
![\"\"](\"\/campus\/wp-content\/uploads\/sites\/2\/2024\/06\/shlee-discount-bond-formula.png\")
<\/figure>\n\n\n\n![\"\"](\"\/campus\/wp-content\/uploads\/sites\/2\/2024\/06\/shlee-coupon-bond-formula.png\")
<\/figure>\n\n\n\n![\"\"](\"\/campus\/wp-content\/uploads\/sites\/2\/2024\/06\/shlee-bond-formula-1.png\")
\nindicates that if valuation date tval<\/sub>\nfall on between interest payment dates ![\"\"](\"\/campus\/wp-content\/uploads\/sites\/2\/2024\/06\/shlee-bond-formula-2.png\")
\nthe simple interest rate calculation is applied to the period tval<\/sub> \u2212 tni<\/sub>, not continously compounding.<\/p>\n\n\n\n![\"\"](\"\/campus\/wp-content\/uploads\/sites\/2\/2024\/06\/CF_schdule-shlee.png\")
<\/figure>\n\n\n\n![\"\"](\"\/campus\/wp-content\/uploads\/sites\/2\/2024\/06\/bond_cash_flow_time_line-shlee.png\")
<\/figure>\n\n\n\n![\"\"](\"\/campus\/wp-content\/uploads\/sites\/2\/2024\/06\/bond_price_formula_application-shlee.png\")
<\/figure>\n\n\n\n