rm(list=ls())

r=1
x = seq(-3,3,length=100)
p0 = dnorm(x,-r,1)
p1 = dnorm(x,r,1)

mu = p1/(p1+p0)

t=0.4
xt = log((1-t)/t)

pdf("ROC1b.pdf")

plot(x,p1,type="l",lwd=4,xlab="",ylab="",yaxt="n",col="purple",xaxt="n") 
abline(h=0,col="grey",lwd=1.5)

segments(xt,-.5,xt,dnorm(xt,r,1),col="purple",lwd=3)

mtext(expression(p[1]),side=2,line=1.5,cex=1.8)
mtext("x",side=1,line=.6,cex=1.5,adj=.9)    ## quiqui

segments(r,-.5, r,dnorm(r,r,1),lwd=2,col="purple") 
axis(side=1, at=1,"1", cex.axis=1.7,line=1,tck=0,col="purple")  

Px=seq(xt,3,length=50)
Py=dnorm(Px,r,1)
polygon(c(Px,rev(Px)),c(Py,rep(0,50)),col="thistle")
segments(r,-.5,r,dnorm(r,r,1),lwd=2,col="purple") 
segments(xt,-.5,xt,dnorm(xt,r,1),col="purple",lwd=2)

segments(xt,-.5,xt,dnorm(xt,r,1),lty=2,col="purple",lwd=3)  
axis(side=1, at=xt,"0.41",col.lab="purple", cex.axis=1.5,
           line=.3,tck=0,lwd=2)  

segments(2,-.5, 2,dnorm(2,r,1),lwd=1,lty=2,col="purple") 
axis(side=1, at=2,"2",tck=0, cex.axis=1.5,line=1)  

segments(-1,-.5, -1,dnorm(-1,r,1),lwd=1,lty=2,col="purple") 
axis(side=1, at=-1,"-1",tck=0, cex.axis=1.5,line=1)


dev.off()