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("ROC1a.pdf")

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

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

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

axis(side=1, at=xt,col.lab="purple","x(t) = .41", cex.axis=1.5,
             line=1,tck=0,lwd=2,col.axis="purple") 

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

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

dev.off()