SwapPricer is on Github
In the previous posts we have seen how easy it is to price interest rate swaps using R. I am honoured to announce that I have decided to put all the functions I have described together into a package that is called…SwapPricer!
Ok, the name is not super original, but it should at least be easy to remember.
You can install it as follows:
# library(devtools) devtools::install_github("DavideMagno/SwapPricer")
The package is still unfortunately not on CRAN but it has an official hexagon. Here it is:
Let me know if you like it in the Disqus form below in the post.
In order to price a swap you just need to run the following code.
library(SwapPricer) SwapPortfolioPricing(SwapPricer::swap.basket, lubridate::ymd(20190414), SwapPricer::df.table)
## # A tibble: 5 x 7 ## swap.id clean.mv dirty.mv accrual.pay accrual.receive par pv01 ## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> ## 1 Swap 25y -8.82e5 -8.75e5 5441. 1379. 0.00771 -12394. ## 2 Swap 30y 2.34e5 1.24e5 -97222. -12470 0.0111 20867. ## 3 Swap 10y 2.22e5 2.36e5 6702. 7361. -0.00138 -5724. ## 4 Swap 2y16y 3.60e5 3.60e5 0 0 0.0118 -11163. ## 5 Swap non … -2.59e6 -2.87e6 -263836. -14681. 0.0107 27914.
You can see that I have used two objects that are delivered with the package:
swap.basket which consists in a 5 swaps portfolio that can be referenced as blueprint for your swap portfolio
df.table this is the discount curve downloaded from Bloomberg as at the 14th of April 2019
We have tested the package using a 500 swaps portfolio and the results, in terms of performance are very encouraging. We analyse them using the amazing profvis tool.
The results are quite interesting as they show we can split the code in three parts:
The first one converts the inputs in tabular form into list of lists of the same format as we described in this post.
The second one is the one that takes the biggest amount of time and it corresponds to the
CashFlowPricingfunction. We have described this function in detail in this post. If we dig into the profvis output, we can notice that the function that takes the most time is the
RQuantlib::advanceone. Unfortunately, this great function doesn’t have a vectorised version hence we need to use functional programming like
mapto calculate multiple future cashflows. Also the
RQuantlib::yearFractionfunction doesn’t have a vectorised version hence I have preferred using another function that does the same but in a vectorised way. This function is
fmdates::year_frac. The difference in performance is sizeable, as you can see from the microbenchmark test below:
today <- Sys.Date() test <- seq(from = 3, to = 50*12, by = 3) cashflows <- today + months(test) microbenchmark::microbenchmark( "Rquantlib" = mapply(function(x) RQuantLib::yearFraction(today, x, 2), cashflows), "fmdates" = fmdates::year_frac(today, cashflows, "act/365") )
## Unit: microseconds ## expr min lq mean median uq max neval ## Rquantlib 2846.262 2939.234 3662.5030 3065.84 3313.7375 20195.020 100 ## fmdates 199.366 219.874 301.3815 274.02 344.1055 1045.776 100
The difference in pricing is quite striking. I have submitted an issue request on the Github page of RQuantLib, let’s hope they’ll want to implement a vectorised version of their cashflow functions.
- The actual pricing of the 500 swaps takes less than 1 second. This is very encouraging because when it will come to curve optimisation (formerly known as bootstrapping) or risk sensitivities and scenarios, parts 1 and 2 will be run once at the beginning of the process while part 3 will be run multiple times.
Finally, I want to warn you: the package is still in a very early version, this is just v0.1.0, hence it is able to price just:
EUR swaps that pay semi-annually on the floating leg
Using a one-curve framework
We are working to improve the toolbox in the next releases
This is the plan of improvements (after I’ll submit the current version on CRAN for everybody to more easily access it - I am thinking, for example, to corporate networks that have the ports to Github blocked):
Extend to all the other currencies, custom frequencies and conventions
Introduce OIS discounting
Calculate additional risk measures like DV01, bucketed DV01 and Convexity
It’s going to be an interesting and long ride, you are welcome on-board!