vv2<-seq(0,1,length=1000) plot(vv2,dbeta(vv2,6,8),type="l") abline(h=.7,col="red") lines(c(.2,.2),c(0,.7)) lines(vv2[vv2<.2],dbeta(vv2[vv2<.2],6,8),type="h",col="green",lwd=5) lines(c(.658,.658),c(0,.7)) lines(vv2[vv2>.658],dbeta(vv2[vv2>.658],6,8),type="h",col="purple",lwd=5)