rm(list=ls())

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

t=.4
xt = log((1-t)/t)

mu = p1/(p1+p0)

pdf("ROC2a.pdf")

plot(x,mu,type="l",lwd=3,xlab="",ylab="",xaxt="n",yaxt="n",col="purple") 
mtext(expression(mu(x)),side=2,line=2,cex=1.5)
mtext("x",side=1,line=2,cex=1.5)

## axis(side=1, at=2.5,"x",cex.axis=1.5,line=2)   #cex=1.5,
## tt=round(t,2)
## axis(side=2,at=t,col="purple","t=0.4",las=1,cex=1.5,line=-3)
dev.off()

pdf("ROC2b.pdf")

fpr = rep(0,100)
tpr = rep(0,100)
tt = rep(0,100)
a = rep(0,100)
b = rep(0,100)
for(i in 1:100){
     tt[i] = mu[i]
     fpr[i] = 1-pnorm(x[i],-r,1)
     tpr[i] = 1-pnorm(x[i],r,1)
     }

plot(fpr,tpr,type="l",xlim=c(0,1),ylim=c(0,1),xlab="",ylab="",
       lwd=4,col="purple", xaxt="n",yaxt="n")
polygon(c(fpr,rev(fpr)),c(tpr,rep(0,100)),col="thistle")
## lines(fpr,tpr,lwd=4,col="purple")
## abline(a=0,b=1,col="purple",lwd=2)
## abline(h=1)
## abline(v=.97)
segments(0,0,.97,1,col="purple",lwd=2)
mtext(text="TPR",cex=1.5,side=2,line=2.5)
mtext(text="FPR",cex=1.5,side=1,line=3)
dev.off()