# používá soubor brno_t.csv z předchozího ukázkového příkladu 
# upraví původní hodnoty tzv. standardizací (mají půměr nula a směr. odch. jedna)
z_JUL <- scale(teploty_brno$JUL)
mean(z_JUL)
sd(z_JUL)
barvy<-ifelse(z_JUL>=0, "red","blue")
plot(teploty_brno$YEAR,z_JUL, type="h", lwd=2, las=1,col=barvy,
	xlab="Roky", ylab="Teplotní anomálie")

# prida horizontalni caru
abline(h=0,lty=2,col="tomato")
# přidá řadu shlazených nízkofrekvenčním filtrem
lines(ksmooth(teploty_brno$YEAR, z_JUL, 'normal', bandwidth = 12), col="darkgreen", lwd=3)

# prida legendu
legend(x="bottomright", legend=c("teplé", "chladné"), bty="n", lwd=2, col=c("red", "blue"))

 # vyznaci abs min
min1.t<-ifelse(z_JUL<=sort(z_JUL)[1],z_JUL,NA)
points(teploty_brno$YEAR,min1.t, col="blue",pch=16)
 
# vyznaci abs max
max1.t<-ifelse(z_JUL>sort(z_JUL)[length(z_JUL)-1],z_JUL,NA)
points(teploty_brno$YEAR,max1.t, col="red",pch=16)