h=hist(rnorm(200),plot=FALSE)
n=length(h$breaks)
level=(h$breaks[n]+h$breaks[n-3])/2
clr <- ifelse(h$breaks < level, "white", "red")[-n]
plot(h, col=clr)
segments(level,0,level,50,lty=2,lwd=2)