p <- cbind(x = seq(0,1,length.out = 100), 
           y = seq(0,1,length.out=100))

plot(p,type="l")
polygon(x=c(0,1,1),y=c(0,1,0), col = "grey")
text(0.6,0.3,"Postive Region")