set.seed(100) #algv?¤?¤rtus juhuarvude loomiseks
pikkused=round(rnorm(20, 170, 10))
pikkused
mean(pikkused)
sd(pikkused)
curve(dnorm, -5, 5)
abline(v=-1)
pnorm(0)
pnorm(-1) #?¼he standardh?¤lbe jagu v?¤hem kui keskv?¤?¤rtus
1-pnorm(-1) #84% andmetest on suuremad kui keskv?¤?¤rtus 
            #miinus ?¼ks standardh?¤lve
pnorm(2)
pnorm(-2)
pnorm(2)-pnorm(-2) #andmete hulk +- 2 standardh?¤lbe kaugusel
pnorm(165, mean=170, sd=10) #31% pikkustest jaotusej?¤rgsetest pikkustest
     #on v?¤iksemad kui 165
qnorm(0.4, mean=170, sd=10) #40 pikkustest on v?¤iksemad kui 167,5

(tabamused=rbinom(150, 10, 0.3))
hist(tabamused, breaks=10, ylim=c(0, 8))
help(pbinom)
pbinom(2, 10, 0.3) #millise tõenäosusega saan sisse kaks või vähem viset kümnest
1-pbinom(2, 10, 0.3) #tõenäosus, et saab sisse rohkem kui kaks viset
qbinom(0.8, 10, 0.3) #80 juhtudest saan sisse 4 või vähem viset
qbinom(0.38, 10, 0.3)

pikkused
pikkused=round(rnorm(50, 170, 10))
mean(pikkused)
sd(pikkused)
t.test(pikkused)
t.test(pikkused, conf.level=0.99)
t.test(pikkused, mu=165)
help(t.test)
t.test(pikkused, mu=169)
t.test(pikkused, mu=169, alternative="greater")

andmed=read.table("http://www.tlu.ee/~jaagup/andmed/ilm/harkutund.txt", sep=",", header=TRUE)
head(andmed)
andmed$TA1H
mean(andmed$TA1H)
mean(andmed$TA1H, na.rm=TRUE)
fivenum(andmed$TA1H, na.rm=TRUE)
attach(andmed) #teeb tabeli tulbad otse kättesaadavaks
Kuu
detach(andmed) #kaotab otsetee tulpadeni

marts=subset(andmed, Kuu==3)$TA1H #märtsikuu keskmised temperatuurid
november=subset(andmed, Kuu==11)$TA1H

fivenum(marts)
fivenum(november)
boxplot(cbind(marts, november))
t.test(marts, november)
t.test(marts, november, conf.level=0.2)
skewness(marts)
skewness(november)
qqnorm(marts)
qqline(marts)
hist(marts)
hist(november, breaks=5)

hist(november, breaks=5)
shapiro.test(marts) #väike p-väärtus näitab, et tegemist pole normaaljaotusega
help(wilcox.test)
wilcox.test(marts, november)
jaanuar=subset(andmed, Kuu==1)$TA1H
which(is.na(jaanuar))
length(jaanuar)
wilcox.test(marts, jaanuar)
which(is.na(andmed$TA1H)) #järjekorranumbrid, millised temperatuurid on puudu
andmed[which(is.na(andmed$TA1H)), ]$TA1H=mean(andmed$TA1H, na.rm=TRUE)
  #Nende järjekorranumbritega kirjete temperatuuritulpadele omistatakse 
  #aasta keskmise temperatuuri väärtus

 
nrow(andmed)
plot(andmed$TAN1H, andmed$TAX1H)
#tervikandmed=andmed[which(!is.na(andmed)),]
tervikandmed=na.omit(andmed)
cor(tervikandmed$TAN1H, tervikandmed$TAX1H)
cor(tervikandmed$TAN1H, tervikandmed$RH1H)
plot(tervikandmed$TAN1H, tervikandmed$RH1H)
abline(lm(tervikandmed$RH1H~tervikandmed$TAN1H))