jaanuar=c(-1.2,0.4,1.1,0.4,1.1,1.6,3.6,5.9,2.9,1.9,6.3,4.2,4,1.9,-2.2,-0.2,2.3,1,-0.4,-0.8,-0.2,-0.7,-0.2,-2.2,-3.2,-4,-6.8,-10.7,-12.1,-11.3,-1)
jada=seq(-16, 7, 4)
plokid=cut(jaanuar, jada)
cbind(table(plokid)) 
table(plokid)/length(jaanuar) #suhteline sagedus
table(plokid)
cumsum(table(plokid)) #kumulatiivne jaotus, kui palju tulemusi on väiksemad kui vastav temperatuur
100*cumsum(table(plokid)/length(jaanuar)) #kumulatiivne jaotus, mitu protsenti tulemusi on väiksemad kui vastav temperatuur
#svg("c:\\jaagup\\2016\\2\\temperatuurid.svg")
pdf("c:\\jaagup\\2016\\2\\temperatuurid.pdf")
hist(jaanuar, breaks=c(-20, -10, -5, 0, 5, 10), main="Jaanuari temperatuuride jaotus")
dev.off()