R/durmaturity.R

In GUIDE: GUI for DErivatives in R

durmaturity <-
function(){
  
  
  #discrate<-0.10
  #maxdiscrate<-0.2
  maxmaturity <- 30 # lets keep this fixed
  faceval<-1000 # lets keep this fixed
  
  
  
  #faceval<-as.numeric(panel$facevalue)
  
  # maturity <- panel$maturity
  
  bonddur<- function(panel,maturity,faceval){
    
    
    couprate<-as.numeric(panel$couprate)/100 
    discrate = as.numeric(panel$discrate)/100
    freq<-1 # change thru radio button
    times<-seq(from=1,by=1,length.out=maturity*freq)
    
    pvfactors=1/(1+discrate/freq)^(freq*times)
    
    
    
    #     effrate = discrate/(100*freq)
    #     effperiods = freq*maturity
    #     pv_coupons<-(coupon/effrate)*(1-(1+effrate)^(-effperiods)) # PV of coupons
    #     pv_face<-faceval*(1+effrate)^(-effperiods) # PV of face value
    #     price<-pv_coupons+pv_face # bond price is the sum of both
    #     price <- round(price,2)
    
    coupon<-couprate*faceval/100
    cashflows <- rep(coupon,maturity*freq)
    cashflows[length(cashflows)] = cashflows[length(cashflows)]+faceval
    
    price<-sum(cashflows*pvfactors)
    dur = sum(cashflows*pvfactors*times)/price
    dur<-round(dur,2)
    return(dur)
  }
  
  
  
  
  my.draw <- function(panel) 
  {
    with(panel, {
      bonddurs<-rep(0,maxmaturity)
      for (maturity in 1:maxmaturity){
        bonddurs[maturity]=bonddur(panel,maturity,faceval)
      }
      
      plot(1:maxmaturity, bonddurs, type='l',lwd=2, ylim=c(1,maxmaturity),xlab="Maturity", ylab="Duration",
           frame = TRUE)
      title(paste("Duration and Maturity"))
      #cat(pv)
    })
    panel
    
  }
  
  
  
  
  
  my.redraw <- function(panel)
  {
    rp.tkrreplot(panel, tkrp)
    panel                                                                       
  }
  
  
  
  
  my.panel <- rp.control("Duration and Maturity")
  rp.doublebutton(my.panel,variable=couprate,step=1,title="Coupon (%  p.a.)",initval=10,range=c(1,20),showvalue=TRUE,action=my.redraw)
  rp.doublebutton(my.panel,variable=discrate,step=1,title="Discount Rate (%  p.a.)",initval=10,range=c(1,20),showvalue=TRUE,action=my.redraw)
  #rp.doublebutton(my.panel,variable=maturity,step=0.25,title="Maturity (Yrs)",initval=10,range=c(1,25),showvalue=TRUE,action=my.redraw)
  rp.tkrplot(panel=my.panel, name=tkrp, plotfun=my.draw, hscale = 2, vscale=1.8)
  
  #rp.do(my.panel, my.draw)
}