TU1KeyRateDuration.java

  1. package org.drip.sample.treasuryfuturesrisk;

  2. import java.util.*;

  3. import org.drip.analytics.date.JulianDate;
  4. import org.drip.feed.loader.*;
  5. import org.drip.historical.sensitivity.TenorDurationNodeMetrics;
  6. import org.drip.market.exchange.TreasuryFuturesContractContainer;
  7. import org.drip.numerical.common.FormatUtil;
  8. import org.drip.service.env.EnvManager;
  9. import org.drip.service.product.TreasuryFuturesAPI;

  10. /*
  11.  * -*- mode: java; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
  12.  */

  13. /*!
  14.  * Copyright (C) 2018 Lakshmi Krishnamurthy
  15.  * Copyright (C) 2017 Lakshmi Krishnamurthy
  16.  * Copyright (C) 2016 Lakshmi Krishnamurthy
  17.  *
  18.  *  This file is part of DRIP, a free-software/open-source library for buy/side financial/trading model
  19.  *      libraries targeting analysts and developers
  20.  *      https://lakshmidrip.github.io/DRIP/
  21.  *  
  22.  *  DRIP is composed of four main libraries:
  23.  *  
  24.  *  - DRIP Fixed Income - https://lakshmidrip.github.io/DRIP-Fixed-Income/
  25.  *  - DRIP Asset Allocation - https://lakshmidrip.github.io/DRIP-Asset-Allocation/
  26.  *  - DRIP Numerical Optimizer - https://lakshmidrip.github.io/DRIP-Numerical-Optimizer/
  27.  *  - DRIP Statistical Learning - https://lakshmidrip.github.io/DRIP-Statistical-Learning/
  28.  *
  29.  *  - DRIP Fixed Income: Library for Instrument/Trading Conventions, Treasury Futures/Options,
  30.  *      Funding/Forward/Overnight Curves, Multi-Curve Construction/Valuation, Collateral Valuation and XVA
  31.  *      Metric Generation, Calibration and Hedge Attributions, Statistical Curve Construction, Bond RV
  32.  *      Metrics, Stochastic Evolution and Option Pricing, Interest Rate Dynamics and Option Pricing, LMM
  33.  *      Extensions/Calibrations/Greeks, Algorithmic Differentiation, and Asset Backed Models and Analytics.
  34.  *
  35.  *  - DRIP Asset Allocation: Library for model libraries for MPT framework, Black Litterman Strategy
  36.  *      Incorporator, Holdings Constraint, and Transaction Costs.
  37.  *
  38.  *  - DRIP Numerical Optimizer: Library for Numerical Optimization and Spline Functionality.
  39.  *
  40.  *  - DRIP Statistical Learning: Library for Statistical Evaluation and Machine Learning.
  41.  *
  42.  *  Licensed under the Apache License, Version 2.0 (the "License");
  43.  *      you may not use this file except in compliance with the License.
  44.  *  
  45.  *  You may obtain a copy of the License at
  46.  *      http://www.apache.org/licenses/LICENSE-2.0
  47.  *  
  48.  *  Unless required by applicable law or agreed to in writing, software
  49.  *      distributed under the License is distributed on an "AS IS" BASIS,
  50.  *      WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  51.  *  
  52.  *  See the License for the specific language governing permissions and
  53.  *      limitations under the License.
  54.  */

  55. /**
  56.  * TU1KeyRateDuration demonstrates the Computation of the Key Rate Duration for the TU1 Treasury Futures.
  57.  *
  58.  * @author Lakshmi Krishnamurthy
  59.  */

  60. public class TU1KeyRateDuration {

  61.     public static final void main (
  62.         final String[] astrArgs)
  63.         throws Exception
  64.     {
  65.         EnvManager.InitEnv ("");

  66.         String strFuturesCode = "TU1";

  67.         String strTreasuryCode = TreasuryFuturesContractContainer.TreasuryFuturesContract (strFuturesCode).code();

  68.         String strTreasuryMarkLocation = "C:\\DRIP\\CreditAnalytics\\Daemons\\Transforms\\TreasuryYieldMarks\\" + strTreasuryCode +
  69.             "BenchmarksReconstituted.csv";
  70.         String strPrintLocation = "C:\\DRIP\\CreditAnalytics\\Daemons\\Transforms\\TreasuryFuturesCloses\\" +
  71.             strFuturesCode + "ClosesReconstitutor.csv";

  72.         CSVGrid csvGridTreasuryMark = CSVParser.StringGrid (
  73.             strTreasuryMarkLocation,
  74.             true
  75.         );

  76.         JulianDate[] adtTreasuryMark = csvGridTreasuryMark.dateArrayAtColumn (0);

  77.         double[] adblYield02Y = csvGridTreasuryMark.doubleArrayAtColumn (1);

  78.         double[] adblYield03Y = csvGridTreasuryMark.doubleArrayAtColumn (2);

  79.         double[] adblYield04Y = csvGridTreasuryMark.doubleArrayAtColumn (3);

  80.         double[] adblYield05Y = csvGridTreasuryMark.doubleArrayAtColumn (4);

  81.         double[] adblYield07Y = csvGridTreasuryMark.doubleArrayAtColumn (5);

  82.         double[] adblYield10Y = csvGridTreasuryMark.doubleArrayAtColumn (6);

  83.         double[] adblYield20Y = csvGridTreasuryMark.doubleArrayAtColumn (7);

  84.         double[] adblYield30Y = csvGridTreasuryMark.doubleArrayAtColumn (8);

  85.         Map<JulianDate, Double> mapTreasuryMark02Y = new TreeMap<JulianDate, Double>();

  86.         Map<JulianDate, Double> mapTreasuryMark03Y = new TreeMap<JulianDate, Double>();

  87.         Map<JulianDate, Double> mapTreasuryMark04Y = new TreeMap<JulianDate, Double>();

  88.         Map<JulianDate, Double> mapTreasuryMark05Y = new TreeMap<JulianDate, Double>();

  89.         Map<JulianDate, Double> mapTreasuryMark07Y = new TreeMap<JulianDate, Double>();

  90.         Map<JulianDate, Double> mapTreasuryMark10Y = new TreeMap<JulianDate, Double>();

  91.         Map<JulianDate, Double> mapTreasuryMark20Y = new TreeMap<JulianDate, Double>();

  92.         Map<JulianDate, Double> mapTreasuryMark30Y = new TreeMap<JulianDate, Double>();

  93.         for (int i = 0; i < adtTreasuryMark.length; ++i) {
  94.             mapTreasuryMark02Y.put (adtTreasuryMark[i], adblYield02Y[i]);

  95.             mapTreasuryMark03Y.put (adtTreasuryMark[i], adblYield03Y[i]);

  96.             mapTreasuryMark04Y.put (adtTreasuryMark[i], adblYield04Y[i]);

  97.             mapTreasuryMark05Y.put (adtTreasuryMark[i], adblYield05Y[i]);

  98.             mapTreasuryMark07Y.put (adtTreasuryMark[i], adblYield07Y[i]);

  99.             mapTreasuryMark10Y.put (adtTreasuryMark[i], adblYield10Y[i]);

  100.             mapTreasuryMark20Y.put (adtTreasuryMark[i], adblYield20Y[i]);

  101.             mapTreasuryMark30Y.put (adtTreasuryMark[i], adblYield30Y[i]);
  102.         }

  103.         CSVGrid csvGrid = CSVParser.StringGrid (
  104.             strPrintLocation,
  105.             true
  106.         );

  107.         JulianDate[] adtSpot = csvGrid.dateArrayAtColumn (0);

  108.         double[] adblCleanPrice = csvGrid.doubleArrayAtColumn (2);

  109.         double[] adblCoupon = csvGrid.doubleArrayAtColumn (3);

  110.         JulianDate[] adtEffective = csvGrid.dateArrayAtColumn (4);

  111.         JulianDate[] adtMaturity = csvGrid.dateArrayAtColumn (5);

  112.         JulianDate[] adtExpiry = csvGrid.dateArrayAtColumn (6);

  113.         int iNumCompute = adtSpot.length;
  114.         JulianDate[] adtEffectiveCompute = new JulianDate[iNumCompute];
  115.         JulianDate[] adtMaturityCompute = new JulianDate[iNumCompute];
  116.         double[] adblCouponCompute = new double[iNumCompute];
  117.         JulianDate[] adtExpiryCompute = new JulianDate[iNumCompute];
  118.         JulianDate[] adtSpotCompute = new JulianDate[iNumCompute];
  119.         double[] adblCleanPriceCompute = new double[iNumCompute];
  120.         double[][] aadblComputeYield = new double[iNumCompute][8];
  121.         String[] astrBenchmarkTenor = new String[] {
  122.             "2Y",
  123.             "3Y",
  124.             "4Y",
  125.             "5Y",
  126.             "7Y",
  127.             "10Y",
  128.             "20Y",
  129.             "30Y"
  130.         };

  131.         for (int i = 0; i < iNumCompute; ++i) {
  132.             adtEffectiveCompute[i] = adtEffective[i];
  133.             adtMaturityCompute[i] = adtMaturity[i];
  134.             adblCouponCompute[i] = adblCoupon[i];
  135.             adtExpiryCompute[i] = adtExpiry[i];
  136.             adtSpotCompute[i] = adtSpot[i];
  137.             adblCleanPriceCompute[i] = adblCleanPrice[i];

  138.             aadblComputeYield[i][0] = mapTreasuryMark02Y.get (adtSpotCompute[i]);

  139.             aadblComputeYield[i][1] = mapTreasuryMark03Y.get (adtSpotCompute[i]);

  140.             aadblComputeYield[i][2] = mapTreasuryMark04Y.get (adtSpotCompute[i]);

  141.             aadblComputeYield[i][3] = mapTreasuryMark05Y.get (adtSpotCompute[i]);

  142.             aadblComputeYield[i][4] = mapTreasuryMark07Y.get (adtSpotCompute[i]);

  143.             aadblComputeYield[i][5] = mapTreasuryMark10Y.get (adtSpotCompute[i]);

  144.             aadblComputeYield[i][6] = mapTreasuryMark20Y.get (adtSpotCompute[i]);

  145.             aadblComputeYield[i][7] = mapTreasuryMark30Y.get (adtSpotCompute[i]);
  146.         }

  147.         List<TenorDurationNodeMetrics> lsTDNM = TreasuryFuturesAPI.HorizonKeyRateDuration (
  148.             strTreasuryCode,
  149.             adtEffectiveCompute,
  150.             adtMaturityCompute,
  151.             adblCouponCompute,
  152.             adtExpiryCompute,
  153.             adtSpotCompute,
  154.             adblCleanPriceCompute,
  155.             astrBenchmarkTenor,
  156.             aadblComputeYield
  157.         );

  158.         System.out.println ("SpotDate,ExpiryDate,CTDName,SpotCTDCleanPrice,ExpiryCTDCleanPrice,SpotGSpread,ExpiryGSpread,SpotYield,ExpiryYield,Parallel,2Y,3Y,4Y,5Y,7Y,10Y,20Y,30Y");

  159.         for (TenorDurationNodeMetrics tdnm : lsTDNM) {
  160.             String strTDNM =
  161.                 tdnm.dateSnap() + "," +
  162.                 tdnm.date ("ExpiryDate") + "," +
  163.                 tdnm.c1 ("CTDName") + "," +
  164.                 FormatUtil.FormatDouble (tdnm.r1 ("SpotCTDCleanPrice"), 1, 5, 100.) + "," +
  165.                 FormatUtil.FormatDouble (tdnm.r1 ("ExpiryCTDCleanPrice"), 1, 5, 100.) + "," +
  166.                 FormatUtil.FormatDouble (tdnm.r1 ("SpotGSpread"), 1, 1, 10000.) + "," +
  167.                 FormatUtil.FormatDouble (tdnm.r1 ("ExpiryGSpread"), 1, 1, 10000.) + "," +
  168.                 FormatUtil.FormatDouble (tdnm.r1 ("SpotYield"), 1, 4, 100.) + "," +
  169.                 FormatUtil.FormatDouble (tdnm.r1 ("ExpiryYield"), 1, 4, 100.) + "," +
  170.                 FormatUtil.FormatDouble (tdnm.r1 ("ParallelKRD"), 1, 4, 1.);

  171.             for (Map.Entry<String, Double> meTDNM : tdnm.krdMap().entrySet())
  172.                 strTDNM += "," + FormatUtil.FormatDouble (meTDNM.getValue(), 1, 4, 1.);

  173.             System.out.println (strTDNM);
  174.         }
  175.     }
  176. }