---
title: "usaldusvahemik"
output: html_document
---


```{r}
ounad10=read.table("http://www.tlu.ee/~jaagup/andmed/muu/ounad/ounad10.txt", header=TRUE, sep=",")
ounad10
t.test(ounad10$diameeter)
t.test(ounad10$diameeter, conf.level=0.9)
 #90% tõenäosusega on sarnaste õunte diameetrite 
 #aritmeetiline keskmine 2.54-3.97
t.test(ounad10$diameeter, conf.level=0.99)
t.test(ounad10$diameeter, alternative="greater")
 #95% tõenäosusega on õunte keskmine diameeter suurem kui 2.54
#Leidke, millisest väärtusest on õunte keskmine diameeter
#väiksem (less) 90% tõenäosusega

#Harjutus: Tehke samad tehed läbi sama kataloogi failiga ounad100.txt
#Võrrelge, kuidas vahemikud erinevad
```

```{r}
ounad100=read.table("http://www.tlu.ee/~jaagup/andmed/muu/ounad/ounad100.txt", header=TRUE, sep=",")
t.test(ounad100$diameeter)
attributes(t.test(ounad100$diameeter))
vahemik=t.test(ounad100$diameeter)[["conf.int"]]
vahemik=round(vahemik, 1)
paste("95% tõenäosusega on õunte keskmine diameeter vahemikus ",
      vahemik[1], " kuni ", vahemik[2], " cm ")

plot(density(ounad100$diameeter), main="Õunte läbimõõt",
       xlab="Keskmine ja keskmise 90% usaldusvahemik")
abline(v=mean(ounad100$diameeter))
vahemik=t.test(ounad100$diameeter, conf.level=0.9)[["conf.int"]]
abline(v=vahemik, lty=2)

#Harjutus: koostage sarnane joonis kümne õunaga
```
Diameetrite võrdlemine õunasortide järgi
```{r}
boxplot(ounad100$diameeter~ounad100$ounasort)
kuldrenetid=ounad100[ounad100$ounasort=="Kuldrenett", "diameeter"]
sibulounad=ounad100[ounad100$ounasort=="Liivi sibul", "diameeter"]
kuldrenetid
sibulounad
plot(density(sibulounad), xlim=c(0, 10))
lines(density(kuldrenetid), col="brown")
#abline(v=mean(sibulounad))
segments(mean(sibulounad), 0, mean(sibulounad), 0.3)
abline(v=mean(kuldrenetid), col="brown")
t.test(kuldrenetid, sibulounad)
# 95% tõenäosusega on kuldrenetid suuremad -0.1518338 kuni  2.0770093 cm
# Kuna võimalus kõikuda on mõlemale poole, siis selle tõenäosusega
# ma ei või midagi väita
t.test(kuldrenetid, sibulounad, conf.level=0.9)
#90% tõenäosusega on kuldrenetid suuremad kui sibulõunad

#Harjutus: võtke 1000 õunaga fail ning näidake, 
#millise tõenäosusega võite väita, et üks õunasort on suurem kui teine
ounad1000=read.table("http://www.tlu.ee/~jaagup/andmed/muu/ounad/ounad1000.txt", header=TRUE, sep=",")
kuldrenetid=ounad1000[ounad1000$ounasort=="Kuldrenett", "diameeter"]
sibulounad=ounad1000[ounad1000$ounasort=="Liivi sibul", "diameeter"]
t.test(kuldrenetid, sibulounad)
#95% tõenäosusega on kuldrenetid suuremad  0.8614071 kuni 1.2355042 cm
t.test(kuldrenetid, sibulounad, conf.level=0.99)
#99% tõenäosusega on kuldrenetid suuremad vähemalt 0.8 cm
```