I have a plot where the color of circles (pch=19) at present is either red (when values in conditioning column is >4.5) and black if values are missing (NA).
points(bb$Year,bb$Tot,col=ifelse(bb$rate<4.5|is.na(bb$rate),'black','red'),pch=19,cex=1.7)
I need a third color, or better, an open circle (pch=21) when values are <=4.5.
Probably simple, but a bit tricky to me.