---
title: "usaldusvahemikud"
output: html_document
---
```{r}
ounad10=read.table("http://www.tlu.ee/~jaagup/andmed/muu/ounad/ounad10.txt", header=TRUE, sep=",")
ounad10
hist(ounad10$diameeter)
t.test(ounad10$diameeter)
```
Selgitusi testi vastuse kohta:
Aritmeetiline keskmine: 3.257
Nende andmete põhjal 95% tõenäosusega (95-l juhul sajast arvan õigesti, et )
on sarnaste õunte
diameetrite aritmeetiline keskmine vahemikus
2.375835 kuni 4.138165
One Sample t-test
data: ounad10$diameeter
t = 8.3615, df = 9, p-value = 1.552e-05
alternative hypothesis: true mean is not equal to 0
95 percent confidence interval:
2.375835 4.138165
sample estimates:
mean of x
3.257
```{r}
t.test(ounad10$diameeter, conf.level=0.9) #Piirdun 90% täpsusega
#saan veidi kitsama vahemiku - 2.542958 kuni 3.971042
#Leidke keskmise diameetri hinnangu usaldusvahemik 80% täpsusega
```
```{r}
ounad100=read.table("http://www.tlu.ee/~jaagup/andmed/muu/ounad/ounad100.txt", header=TRUE, sep=",")
t.test(ounad100$diameeter)
t.test(ounad100$diameeter, mu=3) #pakutav keskmine 2cm
t.test(ounad100$diameeter, mu=3, alternative="greater")
```
p-value = 0.05145
5,1% tõenäosusega ei ole õunte keskmine diameeter suurem kui 3 cm
ehk 94.9% tõenäosusega on õunte keskmine diameeter suurem kui 3 cm
95 percent confidence interval:
2.997822 Inf
95% tõenäosusega on õunte keskmine diameeter suurem kui 2,997
Harjutus:
Võtke tuhande õunaga fail
Kuvage keskmise usaldusvahemikud 99, 95, 80 % täpsusega
Leidke, millisest väärtusest on 90% tõenäosusega õunte keskmine suurem
Leidke, millisest väärtusest on 99% tõenäosusega õunte keskmine väiksem (less)
```{r}
kuldrenetid=ounad100[ounad100$ounasort=="Kuldrenett", "diameeter"]
sibulounad=ounad100[ounad100$ounasort=="Liivi sibul", "diameeter"]
kuldrenetid
sibulounad
t.test(kuldrenetid, sibulounad)
#näha kummagi õunasordi keskmine diameeter
#95% tõenäosusega on nende keskmiste vahe
# -0.1518338 kuni 2.0770093 cm
#Kuna keskmine on nullpunkti mõlemal poolel, siis ma ei saa suuruste
#kohta praegu midagi väita
t.test(kuldrenetid, sibulounad, conf.level=0.9)
#90% tõenäosusega on kulrenetid suuremad
# 0.03841949 kuni 1.88675595 cm
plot(density(sibulounad), xlim=c(0, 10))
lines(density(kuldrenetid), col="green")
abline(v=mean(sibulounad))
abline(v=mean(kuldrenetid), col="green")
```