R/ycfwdrates.R

In Techtonique/esgtoolkit: Toolkit for Monte Carlo Simulations

# Instantaneous forward rates ---------------------------------------------


# Instantaneous forward rates
esgfwdrates <- function(in.maturities, in.zerorates,
                        n, horizon, 
                        out.frequency = c("annual", "semi-annual", 
                                          "quarterly", "monthly", 
                                          "weekly", "daily"),  
                        method = c("fmm", "periodic", "natural", 
                                   "monoH.FC", "hyman", "HCSPL", 
                                   "SW"),
                        ...)
{
  if(is.null(in.maturities) || is.null(in.zerorates))
    stop("Zero rates and maturities must be provided")
  
  if (floor(n) != n) 
    stop("'n' must be an integer")
  
  if (floor(horizon) != horizon) 
    stop("'horizon' must be an integer")
  
  max.in.maturities <- max(in.maturities)
  
  if (horizon > (max.in.maturities + .Machine$double.eps))
    stop("not enough maturities for the provided horizon")
  
  frequency <- match.arg(out.frequency)
  delta <- switch(frequency, 
                  "annual" = 1, 
                  "semi-annual" = 0.5, 
                  "quarterly" = 0.25, 
                  "monthly" = 1/12,
                  "weekly" = 1/52,
                  "daily" = 1/252)
  method <- match.arg(method)
  
  # simply coumpounded zero-coupon price
  pricefromeuribor <- function(t, T, L)
  {
    # Brigo P. 7
    tau <- T-t
    return(as.vector(1/(1 + L*tau)))    
  } 
  
  p <- pricefromeuribor(t = 0, T = in.maturities, L = in.zerorates)
  tt <- seq(from = min(in.maturities), to = max.in.maturities, by = delta)
  nn <- length(tt)
  
  # yc <- ycinter(matsin = in.maturities, matsout = tt, p = p, 
  #               typeres="prices", ...)  
  # ZC.prices <- fitted(yc)
  if (method %in% c("fmm", "periodic", "natural", "monoH.FC", "hyman"))
    ZC.prices <- stats::spline(x = in.maturities, y = p, xout = tt, 
                             method =  method, ...)$y
  
  if (base::identical(method, "HCSPL"))
    ZC.prices <- hermitecubicspline(p = p, 
                                    matsin = in.maturities, 
                                    matsout = tt, 
                                    typeres="prices")$values
  
  if (base::identical(method, "SW"))
    ZC.prices <- tZC_SW(p = p, 
                        u = in.maturities, 
                        t = tt, 
                        UFR = 0.0345, 
                        typeres="prices", 
                        ...)$values
    
  ZC.prices <- c(1 + ((p[1] - 1)/(in.maturities[1] - 0))*(seq(0, 1, by = delta) - 0), 
                 ZC.prices[-1])
  
  tt <- seq(0, max(in.maturities), by = delta)
  
  nn <- length(tt)
  
  fwdrates <- ts(data = (1/delta)*(ZC.prices[-nn]/ZC.prices[-1]-1), 
                 start = 0,
                 deltat = delta)     
  
  return(ts(replicate(n, fwdrates), start = 0,
            deltat = delta))
}