Hi,
Is there a way to use a color ramp when using stplot and simultaneously
assign a specific color for a single value?  For example, I am trying to
plot disease cases and would like 0s to be plotted as grey and the rest of
the cases to be colored by a color ramp.

#this plots cases including 0s as a color ramp:
stplot(stfdf[,,"lyme_cases"], t, col.regions = brewer.pal(9, "YlOrRd"),
cuts=9 )[28]
cuts=c(1,seq(20,160,by=20), 186)

I have tried defining cuts so that any polygon with less than 1 case is
colored grey, but this does not work and instead plots any polygon with <20
cases as grey.

stplot(stfdf[,,"lyme_cases"], t, col.regions = c("grey", brewer.pal(9,
"YlOrRd")), cuts=c(0.1,seq(20,160,by=20), 186) )[28]

Thank you very much for your help!

	[[alternative HTML version deleted]]

