#####  Figures EFFICIENCY (score is efficiency)
# actual values. Add this script to the bioeq....R
setwd("~/DocInvest/MyPapers/sqj2012/datasets/razali")
ci1score <- c(0.02, 0.15)
ci2score <- c(0.09, 0.48)

ciq1score <- c(0.25 , 0.61)
ciq2score <- c(-0.09,0.18)
ciq3score <- c(0.14 , 0.79 )
ciq4score <- c(-0.43 , 0.54 )
ciq5score <- c(-0.14 , 0.18 )
ciq6score <- c(0.02 , 0.58 )

# DecisionCI.mean at 20%=(,)
# 
# > cat("DecisionCI.mean at 25%=(",dec25.low,",",dec25.up,")",sep="","\n")
# DecisionCI.mean at 25%=(,)
# 
# > cat("DecisionCI.mean at 50%=(",dec5.low,",",dec5.up,")",sep="","\n")
# DecisionCI.mean at 50%=(,)

## from other R scripts
lowdelta20 <- -0.1327805
lowdelta25 <- -0.1659756
lowdelta50 <- -0.3319512
updelta20 <-  0.1327805
updelta25 <-  0.1659756
updelta50 <-  0.3319512

limxscore <- 0.9

ticky1 <- 0.2
ticky2 <- 0.7
ticky3 <- 1.2
ticky4 <- 1.7
ticky5 <- 2.2
ticky6 <- 2.7
ticky7 <- 3.2
ticky8 <- 3.7

ticktextx1 <- -0.65
ticktextx2 <- -0.5
ticktextx3 <- lowdelta50


namethisfile <- paste(getwd(),"/Fig10",".pdf",sep="")
pdf(file=namethisfile, height=5.5, width=4.5) #, paper="a4",pagecentre=TRUE) #, pagecentre=TRUE)

par(mar=c(3,0,0,0))
plot(0,0,type="n", xlim=c(-limxscore,limxscore), ylim=c(0,6),axes=FALSE, ann=FALSE)

equivintervals <- c(-limxscore,limxscore,lowdelta50,lowdelta20,0,updelta20,updelta50)
equivintervals <- sort(equivintervals)

labels <- c(
  -limxscore,
  expression(-Delta[paste(50,'%',sep="")]),
  expression(-Delta[paste(20,'%',sep="")]),
  0,
  expression(Delta[paste(20,'%',sep="")]),
#   expression(Delta[paste(25,'%',sep="")]),
  expression(Delta[paste(50,'%',sep="")]),
  limxscore)
axis(1, at=equivintervals, labels = labels, pos=0)
mtext(expression(mu[T]-mu[R]),at=0,side=1,line=2)
#mtext("UML-B",at=updelta50,side=1,line=2)

#abline(v=equivintervals,lty="dotted")
segments(c(lowdelta50,updelta50),c(0,0),c(lowdelta50,updelta50),c(4,4), lty=3,lwd=1)
segments(c(lowdelta20,updelta20),c(0,0),c(lowdelta20,updelta20),c(1,1), lty=3,lwd=1)

segments(ci1score[1],0.2,ci1score[2],lty=1,lwd=5, lend=2)
text(ci1score[1],0.2,"|",offset=0.5,cex=1.1,adj=1,font=2)
text(ci1score[2],0.2,"|",offset=0.5,cex=1.1,adj=0,font=2)
text(ticktextx1,0.2,
     "Exp1 Comprehension",
     offset=0.5,cex=0.75,adj=0)
text(0.08,0.2,
     paste("(",
           round(ci1score[1],2),",",
           round(ci1score[2],2),
           ")"),
     pos=3,offset=0.5,cex=0.7,adj=1)

segments(ci2score[1],ticky2,ci2score[2],lty=1,lwd=5, lend=2)
text(ci2score[1],ticky2,"|",offset=0.5,cex=1.1,adj=1,font=2)
text(ci2score[2],ticky2,"|",offset=0.5,cex=1.1,adj=0,font=2)
text(ticktextx1,ticky2,
     "Exp1 Modifications",
     offset=0.5,cex=0.75,adj=0)
text(0.25,ticky2,
     paste("(",
           round(ci2score[1],2),",",
           round(ci2score[2],2),
           ")"),
     pos=3,offset=0.5,cex=0.7,adj=1)

segments(ciq1score[1],ticky3,ciq1score[2],lty=1,lwd=5, lend=2)
text(ciq1score[1],ticky3,"|",offset=0.5,cex=1.1,adj=1,font=2)
text(ciq1score[2],ticky3,"|",offset=0.5,cex=1.1,adj=0,font=2)
text(ticktextx2,ticky3,
     "Exp2 Q1",
     offset=0.5,cex=0.75,adj=1)
text(0.4,ticky3,
     paste("(",
           round(ciq1score[1],2),",",
           round(ciq1score[2],2),
           ")"),
     pos=3,offset=0.5,cex=0.7,adj=1)

segments(ciq2score[1],ticky4,ciq2score[2],lty=1,lwd=5, lend=2)
text(ciq2score[1],ticky4,"|",offset=0.5,cex=1.1,adj=1,font=2)
text(ciq2score[2],ticky4,"|",offset=0.5,cex=1.1,adj=0,font=2)
text(ticktextx2,ticky4,
     "Exp2 Q2",
     offset=0.5,cex=0.75,adj=1)
text(0.05,ticky4,
     paste("(",
           round(ciq2score[1],2),",",
           round(ciq2score[2],2),
           ")"),
     pos=3,offset=0.5,cex=0.7,adj=1)


segments(ciq3score[1],ticky5,ciq3score[2],lty=1,lwd=5, lend=2)
text(ciq3score[1],ticky5,"|",offset=0.5,cex=1.1,adj=1,font=2)
text(ciq3score[2],ticky5,"|",offset=0.5,cex=1.1,adj=0,font=2)
text(ticktextx2,ticky5,
     "Exp2 Q3",
     offset=0.5,cex=0.75,adj=1)
text(0.5,ticky5,
     paste("(",
           round(ciq3score[1],2),",",
           round(ciq3score[2],2),
           ")"),
     pos=3,offset=0.5,cex=0.7,adj=1)

segments(ciq4score[1],ticky6,ciq4score[2],lty=1,lwd=5, lend=2)
text(ciq4score[1],ticky6,"|",offset=0.5,cex=1.1,adj=1,font=2)
text(ciq4score[2],ticky6,"|",offset=0.5,cex=1.1,adj=0,font=2)
text(ticktextx2,ticky6,
     "Exp2 Q4",
     offset=0.5,cex=0.75,adj=1)
text(0,ticky6,
     paste("(",
           round(ciq4score[1],2),",",
           round(ciq4score[2],2),
           ")"),
     pos=3,offset=0.5,cex=0.7,adj=1)

segments(ciq5score[1],ticky7,ciq5score[2],lty=1,lwd=5, lend=2)
text(ciq5score[1],ticky7,"|",offset=0.5,cex=1.1,adj=1,font=2)
text(ciq5score[2],ticky7,"|",offset=0.5,cex=1.1,adj=0,font=2)
text(ticktextx2,ticky7,
     "Exp2 Q5",
     offset=0.5,cex=0.75,adj=1)
text(0,ticky7,
     paste("(",
           round(ciq5score[1],2),",",
           round(ciq5score[2],2),
           ")"),
     pos=3,offset=0.5,cex=0.7,adj=1)

segments(ciq6score[1],ticky8,ciq6score[2],lty=1,lwd=5, lend=2)
text(ciq6score[1],ticky8,"|",offset=0.5,cex=1.1,adj=1,font=2)
text(ciq6score[2],ticky8,"|",offset=0.5,cex=1.1,adj=0,font=2)
text(ticktextx2,ticky8,
     "Exp2 Q6",
     offset=0.5,cex=0.75,adj=1)
text(0.25,ticky8,
     paste("(",
           round(ciq6score[1],2),",",
           round(ciq6score[2],2),
           ")"),
     pos=3,offset=0.5,cex=0.7,adj=1)


dev.off()  # for pdf et al