---
title: "ounad"
output: html_document
---

```{r}
  ounad=read.table("http://www.tlu.ee/~jaagup/andmed/muu/ounad/antoonovka2.txt",
                    sep=",", header=TRUE)
  head(ounad)
  hist(ounad$august)
  density(ounad$august)
  plot(density(ounad$august), xlim=c(0, 10))
  lines(density(ounad$september), col="blue")
  abline(v=mean(ounad$august))
  abline(v=mean(ounad$september), col="blue")
  
```
 Juhuslikult valitud 100 õuna põhjal võib 95% tõenäosusega
 väita, et selle aia Antoonovkate keskmine diameeter augustis jääb 4,6 ja 5 cm vahele.
```{r}
  t.test(ounad$august)
  t.test(ounad$august, conf.level=0.99)
  t.test(ounad$august, conf.level=0.9)

  plot(density(ounad$august))
  abline(v=t.test(ounad$august)$conf.int[[1]])
  abline(v=t.test(ounad$august)$conf.int[[2]])
```

Võrdlemine väärtusega
p-value = 0.04012 näitab, et tõenäosus õunte keskmise diameetri võrdumiseks viiega on 4%
96% tõenäosusega võin väita, et selle aia selle aasta 
Antoonovkate augusti keskmine diameeter on väiksem kui 5cm
```{r}
  t.test(ounad$august, mu=5)
```

Mis tõenäosusega võiksin väita, et selle aia selle aasta 
Antoonovkate augusti keskmine diameeter on väiksem kui 5,3 cm

```{r}
  t.test(ounad$august, mu=5.3)

```
p-value = 1.076e-06 = nullhüpoteesi kehtmise tõenäosus.

0.9999989 tõenäosusega ehk 99.99989% tõenäosusega võin selle juhuvalimi põhjal väita, et augusti Antoonovkate diameeter aias on väiksem kui 5.3


Suvalised õunad augustis ja septembris.
95% tõenäosusega võin väita, et augustis olid õunad 2,1-1,6 cm väiksema diameetriga
```{r}
  t.test(ounad$august, ounad$september)

```

Samad õunad mõõdetud kahel kuul
95% tõenäosusega võin väita, et augustis olid õunad 1,8-2,0 cm väiksema diameetriga
```{r}
  t.test(ounad$august, ounad$september, paired=TRUE)

```
Harjutus:

http://www.tlu.ee/~jaagup/andmed/muu/ounad/ounad100.txt

Näidake, mida saate nende andmete ja t-testi põhjal väita Kuldreneti ja Liivi Sibula kohta

```{r}
  ounad2=read.table("http://www.tlu.ee/~jaagup/andmed/muu/ounad/ounad100.txt", sep=",", header=TRUE);
  kuldrenetid=ounad2[ounad2$ounasort=="Kuldrenett", "diameeter"]
  kuldrenetid
  t.test(kuldrenetid)
  hist(kuldrenetid, ylim=c(0, 15))
```


```{r}
  subset(
      read.table(
         "http://www.tlu.ee/~jaagup/andmed/muu/ounad/ounad100.txt", 
         sep=",", header=TRUE), ounasort=="Kuldrenett"
  )$diameeter
  tihedused=read.table("http://www.tlu.ee/~jaagup/andmed/muu/ounad/tihedused.txt", sep=",", header=TRUE)
  ounad2
  koos=merge(ounad2, tihedused) #ainuke kattuv tulp on ounasort
  koos=merge(ounad2, tihedused, by.x="ounasort", by.y="ounasort") 
  koos$mass=(4/3*pi*(koos$diameeter/2)**3)*koos$tihedus
  head(koos)
```