MARbest3worst3.r

### MARbest3worst3


data=read.csv("C:/empiricalResultsData/rbind43DataFrame526.csv") # use "rbind43DataFrame526.csv" as data set
dataframe2=data.frame(data,rowMeans(data[,3:14]))
dRanking=dataframe2[order(dataframe2[,15]),] 

dateC=c("2007-10-23", "2007-11-02",  "2007-11-09", "2007-12-06", "2008-01-11",  "2008-02-04", "2008-02-22", "2008-03-18", "2008-04-04",  "2008-04-07",  "2008-05-30",  "2008-07-01")
plot(1:12,seq(min(data[1:43,3:14]),max(data[1:43,3:14]),length.out=12),xlab="",ylab="", col="White",axes=FALSE, xaxt="n", yaxt="n")
box() 
axis(1,at=1:12, labels=FALSE)
axis(2,at=c(0.272800, 1.649975, 3.027150, 4.404325, 5.781500))  
mtext("Relative Difference Measure (RDM)", side=2, line=2.7, at=3.027150, font=2,cex=2)  
mtext("Pricing Date", side=1, line=3.7, at=6, font=2,cex=2)  
text(cex=1, x=1:12-.45, y=0.272800-.6, dateC, xpd=TRUE, srt=30)  
#####

lines(1:12,dRanking[1,3:14],pch=1,lty=1,col=1, lwd=9)
lines(1:12,dRanking[2,3:14],pch=1,lty=2,col=2, lwd=9)
lines(1:12,dRanking[3,3:14],pch=1,lty=3,col=3, lwd=9)
lines(1:12,dRanking[41,3:14],pch=1,lty=4,col=4, lwd=9)
lines(1:12,dRanking[42,3:14],pch=1,lty=5,col=5, lwd=9)
lines(1:12,dRanking[43,3:14],pch=1,lty=6,col=6, lwd=9)

legend('topright', c("C_g-j (best)","C_g-g","C_f-j","C_t38","C_t-t","C_c (worst)"), lty=c(1:6),col=c(1:6), lwd=9,cex=1.5)
###