Skip to content
Prev 349775 / 398513 Next

stripchart

Hi Catalin,
I'm not quite sure which values you want to use, but does this start to do
what you want?

Jq90<-quantile(spei$J,0.9)
Jq10<-quantile(spei$J,0.1)
require(plotrix)
plot(spei$year,spei$J,xlab="Year",ylab="Temerature anomaly",
 main="Temperature anomalies by year (1901 to 2013)",
 pch=ifelse(spei$J>=Jq90|spei$J<=Jq10,19,1),
 col=color.scale(spei$J,extremes=c("blue","red")))
abline(h=c(Jq10,Jq90),lty=2)
mtext(c("q10","q90"),side=4,line=0.5,at=c(Jq10,Jq90),las=1)

The stripchart that I got was just three wodges of points.

Jim


On Tue, Apr 7, 2015 at 6:01 AM, catalin roibu <catalinroibu at gmail.com>
wrote: