pikkused=c(165, 166, 166, 168, 169, 169, 173, 174, 175, 190)
mean(pikkused)
## [1] 171.5
median(pikkused)
## [1] 169
barplot(pikkused, ylim=c(150, 200))
plot(pikkused, ylim=c(150, 200))
abline(h=mean(pikkused))
abline(h=median(pikkused), lty=2)
boxplot(pikkused)
hist(pikkused)
fivenum(pikkused)
## [1] 165 166 169 174 190
#Lisage mõned suured pikkused,
#vaadake, kuidas muutub boxploti kuju
pikkused=c(pikkused, 173, 178, 178, 179)
boxplot(pikkused)
hist(pikkused)
ilm=read.table("http://www.tlu.ee/~jaagup/andmed/ilm/harkutund.txt",
header=TRUE, sep=",")
jaantemp=ilm[ilm$Kuu==1, "TA1H"]
hist(jaantemp)
boxplot(jaantemp)
fivenum(jaantemp)
## [1] -18.6 -12.8 -8.4 0.3 6.3
#25% mõõtmistest oli 2014. aasta jaanuaris tempeartuur 0.3 kuni 6.3 kraadi
#50% külmematest mõõtmistest oli -18.6 kuni -8.4 kraadi
#50% soojematest mõõtmistest oli -8.4 kuni 6.3 kraadi
#50% keskmise väärtusega mõõtmistest oli -12.8 kuni 0.3 kraadi
#75% mõõtmistest oli temperatuuriga alla 0.3 kraadi
jaanuar=ilm[ilm$Kuu==1, ]
boxplot(jaanuar$TA1H~jaanuar$Kell, las=2)
#Kirjutage igaüks kaks lauset,
#mida saate joonise põhjal väita
Kell 3 öösel oli 50% mõõtmistest temperatuur -12 kuni 0 kraadi
for(kuunr in 1:12){
print(kuunr)
kuuandmed=ilm[ilm$Kuu==kuunr, ]
boxplot(kuuandmed$TA1H~kuuandmed$Kell, las=2, main=paste("Kuu nr ",kuunr),
xlab="Kellaaeg", ylab="Temperatuur")
}
## [1] 1
## [1] 2
## [1] 3
## [1] 4
## [1] 5
## [1] 6
## [1] 7
## [1] 8
## [1] 9
## [1] 10
## [1] 11
## [1] 12
Kuvage kuude kaupa iga päeva kohta keskmise tuulekiiruse (WS1H) jaotus
for(kuunr in 1:12){
print(kuunr)
kuuandmed=ilm[ilm$Kuu==kuunr, ]
boxplot(kuuandmed$WS1H~kuuandmed$Paev, las=2, main=paste("Kuu nr ",kuunr),
xlab="Kuupäev", ylab="Temperatuur")
}
## [1] 1
## [1] 2
## [1] 3
## [1] 4
## [1] 5
## [1] 6
## [1] 7
## [1] 8
## [1] 9
## [1] 10
## [1] 11
## [1] 12