hist(rnorm(3000, m=.15, sd=.025), 20, col=rgb(0,1,0,.25), xlim=c(0,1),
     xlab="Y = 0 (3000:Green)  ||   Y = 1 (1000:Red)",
     main="Distribution of Predicted Probability (Y = 1)")
hist(rnorm(1000, m=.35, sd=.025), 20, col=rgb(1,0,0,.25), add=TRUE)
abline(v=0.5,col='lightgray')


With the Distribution of Predicted Probability shown above :