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")