joonised

  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

  #hist(ilm$RH1H)
  hist(ilm$RH1H, 
       xlab="Õhuniiskus protsentides", ylab="Esinemistundide arv",
       main="Õhuniiskus Tallinn-Harku jaamas 2014")

  ?hist

## starting httpd help server ...

##  done

  #las - määrab telgede siltide suundi 0..3
 hist(ilm$RH1H, labels=TRUE, las=1, breaks=5,
       ylim=c(0, 7000), xlim=c(40, 100) ) 

  hist(ilm$RH1H[ilm$RH1H>40])

  hist(ilm[ilm$Kuu==1, "RH1H"])

Kaks histogrammi kõrvuti ühel joonisel

  par(mfrow=c(1, 2))
  hist(ilm[ilm$Kuu==1, "RH1H"])
  hist(ilm[ilm$Kuu==7, "RH1H"])

  par(mfrow=c(1, 1))

Kordused

 for(kuunr in 1:12){
   print(kuunr)
 }

## [1] 1
## [1] 2
## [1] 3
## [1] 4
## [1] 5
## [1] 6
## [1] 7
## [1] 8
## [1] 9
## [1] 10
## [1] 11
## [1] 12

Tundide arv jaanuaris, kus õhuniiskus ületab 90%

  length(ilm[ilm$Kuu==1 & ilm$RH1H>90, "RH1H"])

## [1] 402

for (kuunr in 1:12){
  print(paste(kuunr, ' ', 
      length(ilm[ilm$Kuu==kuunr & ilm$RH1H>90, "RH1H"])))
}

## [1] "1   402"
## [1] "2   453"
## [1] "3   198"
## [1] "4   120"
## [1] "5   252"
## [1] "6   298"
## [1] "7   165"
## [1] "8   289"
## [1] "9   319"
## [1] "10   351"
## [1] "11   559"
## [1] "12   501"

Suurima niiskete tundide arvu kuus leidmine

kogused=c()
for(kuunr in 1:12){
  kogused[kuunr]=length(ilm[ilm$Kuu==kuunr & ilm$RH1H>95, "RH1H"])
}
kogused

##  [1] 239 300 146  56 151 219  97 187 240 281 324 245

max(kogused)

## [1] 324

Joonised kuude kaupa vastavalt maksimumkogusele

 kuunimed=c("jaanuar", "veebruar", "märts",
            "aprill", "mai", "juuni", 
            "juuli", "august", "september", 
            "oktoober", "november", "detsember")
 par(mfrow=c(4, 3))
 for (kuunr in 1:12){
   hist(ilm[ilm$Kuu==kuunr, "RH1H"], ylim=c(0, max(kogused)), main=kuunimed[kuunr])
 }

 par(mfrow=c(1, 2))

 par(mfrow=c(4, 3))
 for (kuunr in 1:12){
   hist(ilm[ilm$Kuu==kuunr, "RH1H"])
 }

 par(mfrow=c(1, 2))