---
title: "joonised"
output: html_document
---

```{r}
  ilm=read.table("http://www.tlu.ee/~jaagup/andmed/ilm/harkutund.txt" , header=TRUE, sep=",")
  head(ilm)
  #hist(ilm$RH1H)
  hist(ilm$RH1H, 
       xlab="Õhuniiskus protsentides", ylab="Esinemistundide arv",
       main="Õhuniiskus Tallinn-Harku jaamas 2014")
  ?hist
  #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

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

Kordused


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

Tundide arv jaanuaris, kus õhuniiskus ületab 90%
```{r}
  length(ilm[ilm$Kuu==1 & ilm$RH1H>90, "RH1H"])
for (kuunr in 1:12){
  print(paste(kuunr, ' ', 
      length(ilm[ilm$Kuu==kuunr & ilm$RH1H>90, "RH1H"])))
}
```
  Suurima niiskete tundide arvu kuus leidmine
```{r}
kogused=c()
for(kuunr in 1:12){
  kogused[kuunr]=length(ilm[ilm$Kuu==kuunr & ilm$RH1H>95, "RH1H"])
}
kogused
max(kogused)
```
  Joonised kuude kaupa vastavalt maksimumkogusele
```{r fig.height=12}
 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],
        xlab="%", ylab="N")
 }
 par(mfrow=c(1, 2))
 #lattice
 #ggplot2
```


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


```{r}
   tapply(ilm$RH1H, ilm$Kuu, min, na.rm=TRUE)
   barplot(tapply(ilm$RH1H, ilm$Kuu, min, na.rm=TRUE))
   kuunimed=c("jaanuar", "veebruar", "märts",
            "aprill", "mai", "juuni", 
            "juuli", "august", "september", 
            "oktoober", "november", "detsember")
   kuivimad=tapply(ilm$RH1H, ilm$Kuu, min, na.rm=TRUE)
   barplot(kuivimad, names.arg = kuunimed, las=2,
        main="Kuivim mõõdetud tund kuus")
   
   #Koostage tulpdiagramm kuude keskmiste niiskustega (mean)
   
```

 Karp ja vurrud
```{r}
  boxplot(ilm$RH1H)
  boxplot(ilm[ilm$Kuu==1, "RH1H"])
  boxplot(ilm$RH1H~ilm$Kuu)

```