lapsed=read.table("http://www.tlu.ee/~jaagup/andmed/muu/5klass.txt", header=TRUE, sep=",")
lapsed
## eesnimi pikkus mass sugu
## 1 Juku 170 45 m
## 2 Kati 160 35 n
## 3 Mati 160 72 m
## 4 Madis 165 53 m
## 5 Mati 163 60 m
## 6 Katrin 165 43 n
## 7 Siim 141 38 m
## 8 Martin 159 46 m
## 9 Kadri 164 57 n
## 10 Katariina 148 35 n
## 11 Maria 143 38 n
## 12 Marta 169 55 n
## 13 Madis 156 65 m
## 14 Mihkel 165 69 m
## 15 Tiina 170 38 n
## 16 Miia 145 68 n
## 17 Siim 151 55 m
## 18 Priit 156 63 m
## 19 Kristjan 164 59 m
## 20 Kristi 155 53 n
## 21 Kristiina 158 62 n
## 22 Killu 164 49 n
## 23 Mart 170 69 m
## 24 Kert 143 36 m
## 25 Gert 152 67 m
## 26 Lauri 156 53 m
## 27 Moonika 164 58 n
## 28 Jaanika 165 59 n
## 29 Jaanus 164 63 m
## 30 Jaan 162 65 m
p=sort(lapsed$pikkus)
p
## [1] 141 143 143 145 148 151 152 155 156 156 156 158 159 160 160 162 163
## [18] 164 164 164 164 164 165 165 165 165 169 170 170 170
p[1] #esimene laps
## [1] 141
length(p) #laste arv
## [1] 30
p[length(p)] #Pikim laps
## [1] 170
#Leidke lühemalt poolt kolmanda lapse pikkus
p[3]
## [1] 143
#Leidke pikemalt poolt neljanda lapse pikkus
p[28]
## [1] 170
p[30-3]
## [1] 169
p[length(p)-3]
## [1] 169
p[1:5] #Viis lühemat last
## [1] 141 143 143 145 148
p[(length(p)-4):length(p)] #Viis pikemat last
## [1] 165 169 170 170 170
p[26:30] #sama 30 lapse puhul
## [1] 165 169 170 170 170
#Leidke kümne lühema lapse pikkuste keskmine
1:10
## [1] 1 2 3 4 5 6 7 8 9 10
p[1:10]
## [1] 141 143 143 145 148 151 152 155 156 156
mean(p[1:10])
## [1] 149
#Leidke kümne pikema lapse pikkuste keskmine
(length(p)-9):length(p)
## [1] 21 22 23 24 25 26 27 28 29 30
p[(length(p)-9):length(p)]
## [1] 164 164 165 165 165 165 169 170 170 170
mean(p[(length(p)-9):length(p)])
## [1] 166.7
p
## [1] 141 143 143 145 148 151 152 155 156 156 156 158 159 160 160 162 163
## [18] 164 164 164 164 164 165 165 165 165 169 170 170 170
eelmisepikkused=c(p[2:length(p)], NA)
eelmisepikkused
## [1] 143 143 145 148 151 152 155 156 156 156 158 159 160 160 162 163 164
## [18] 164 164 164 164 165 165 165 165 169 170 170 170 NA
vahed=eelmisepikkused-p
vahed
## [1] 2 0 2 3 3 1 3 1 0 0 2 1 1 0 2 1 1 0 0 0 0 1 0
## [24] 0 0 4 1 0 0 NA
table(vahed)
## vahed
## 0 1 2 3 4
## 13 8 4 3 1
#Leidke pikkused ja pikkuste vahed vaid poiste kohta rivis
poistepikkused=rev(sort(lapsed[lapsed$sugu=="m", "pikkus"]))
poistepikkused
## [1] 170 170 165 165 164 164 163 162 160 159 156 156 156 152 151 143 141
ilmaesimeseta=poistepikkused[2:length(poistepikkused)]
ilmaesimeseta=c(ilmaesimeseta, NA) #et oleks sama pikk loetelu arvutamiseks
ilmaesimeseta
## [1] 170 165 165 164 164 163 162 160 159 156 156 156 152 151 143 141 NA
vaheeelmisega=poistepikkused-ilmaesimeseta
vaheeelmisega
## [1] 0 5 0 1 0 1 1 2 1 3 0 0 4 1 8 2 NA
table(vaheeelmisega) #vahede esinemissagedused
## vaheeelmisega
## 0 1 2 3 4 5 8
## 5 5 2 1 1 1 1
ilm=read.table("http://www.tlu.ee/~jaagup/andmed/ilm/harkutund.txt", header=TRUE, sep=",")
head(ilm)
## Kuu Paev Kell PR1H RH1H TA1H TAN1H TAX1H WD1H WS1H WSX1H
## 1 1 1 0:00 0 86 4.4 4.4 4.5 221 5.6 10.2
## 2 1 1 1:00 0 87 4.3 4.1 4.5 214 4.3 9.0
## 3 1 1 2:00 0 89 3.9 3.8 4.1 205 3.9 9.1
## 4 1 1 3:00 0 89 3.6 3.4 3.8 208 4.3 10.5
## 5 1 1 4:00 0 89 3.3 3.0 3.4 211 4.4 10.1
## 6 1 1 5:00 0 88 2.8 2.6 3.0 207 4.4 11.9
#Tulp TA1H tunni keskmine temperatuur
jaanuaritemp=ilm[ilm$Kuu==1, "TA1H"]
hist(jaanuaritemp) #histogramm, sagedused kraadivahemike kaupa
jargminetund=c(jaanuaritemp[2:length(jaanuaritemp)], NA)
head(jargminetund)
## [1] 4.3 3.9 3.6 3.3 2.8 2.5
vahed=jargminetund-jaanuaritemp #temperatuurierinevus
head(vahed)
## [1] -0.1 -0.4 -0.3 -0.3 -0.5 -0.3
hist(vahed)
