Õunte andmed

ounad=read.table("http://www.tlu.ee/~jaagup/andmed/muu/ounad/ounad100.txt", header=TRUE, sep=",")
head(ounad)

##      ounasort diameeter
## 1 Liivi sibul      3.68
## 2 Liivi sibul      1.79
## 3  Kuldrenett      3.99
## 4 Liivi sibul      3.66
## 5 Liivi sibul      1.95
## 6 Liivi sibul      3.26

table(ounad$ounasort)

## 
##  Kuldrenett Liivi sibul 
##          24          76

pie(table(ounad$ounasort))

barplot(table(ounad$ounasort))

table(ounad$ounasort, ounad$diameeter<5)

##              
##               FALSE TRUE
##   Kuldrenett      9   15
##   Liivi sibul     2   74

barplot(table(ounad$diameeter>5, ounad$ounasort))

barplot(table(ounad$diameeter>5, ounad$ounasort), legend.text=c("alla 5", "üle 5"))

barplot(table(ounad$diameeter>5, ounad$ounasort), beside=TRUE)

boxplot(ounad$diameeter)

boxplot(ounad$diameeter~ounad$ounasort)

hist(ounad$diameeter)

hist(ounad[ounad$ounasort=="Liivi sibul", "diameeter"], col="red")
hist(ounad[ounad$ounasort=="Kuldrenett", "diameeter"], add=TRUE, col="yellow")

table(round(ounad$diameeter))

## 
##  0  1  2  3  4  5  6  7  8 10 
##  1  5 24 37 21  2  5  3  1  1

table(round(ounad$diameeter), ounad$ounasort)

##     
##      Kuldrenett Liivi sibul
##   0           1           0
##   1           3           2
##   2           5          19
##   3           2          35
##   4           4          17
##   5           0           2
##   6           4           1
##   7           3           0
##   8           1           0
##   10          1           0

t(table(round(ounad$diameeter), ounad$ounasort))

##              
##                0  1  2  3  4  5  6  7  8 10
##   Kuldrenett   1  3  5  2  4  0  4  3  1  1
##   Liivi sibul  0  2 19 35 17  2  1  0  0  0

barplot(t(table(round(ounad$diameeter), ounad$ounasort)))

prop.table(t(table(round(ounad$diameeter), ounad$ounasort)))

##              
##                  0    1    2    3    4    5    6    7    8   10
##   Kuldrenett  0.01 0.03 0.05 0.02 0.04 0.00 0.04 0.03 0.01 0.01
##   Liivi sibul 0.00 0.02 0.19 0.35 0.17 0.02 0.01 0.00 0.00 0.00

prop.table(t(table(round(ounad$diameeter), ounad$ounasort)), 2)

##              
##                        0          1          2          3          4
##   Kuldrenett  1.00000000 0.60000000 0.20833333 0.05405405 0.19047619
##   Liivi sibul 0.00000000 0.40000000 0.79166667 0.94594595 0.80952381
##              
##                        5          6          7          8         10
##   Kuldrenett  0.00000000 0.80000000 1.00000000 1.00000000 1.00000000
##   Liivi sibul 1.00000000 0.20000000 0.00000000 0.00000000 0.00000000

prop.table(t(table(round(ounad$diameeter), ounad$ounasort)), 1)

##              
##                        0          1          2          3          4
##   Kuldrenett  0.04166667 0.12500000 0.20833333 0.08333333 0.16666667
##   Liivi sibul 0.00000000 0.02631579 0.25000000 0.46052632 0.22368421
##              
##                        5          6          7          8         10
##   Kuldrenett  0.00000000 0.16666667 0.12500000 0.04166667 0.04166667
##   Liivi sibul 0.02631579 0.01315789 0.00000000 0.00000000 0.00000000

prop.table(table(ounad$ounasort, ounad$diameeter<5))

##              
##               FALSE TRUE
##   Kuldrenett   0.09 0.15
##   Liivi sibul  0.02 0.74

prop.table(table(ounad$ounasort, ounad$diameeter<5), 1)

##              
##                    FALSE       TRUE
##   Kuldrenett  0.37500000 0.62500000
##   Liivi sibul 0.02631579 0.97368421

prop.table(table(ounad$ounasort, ounad$diameeter<5), 2)

##              
##                   FALSE      TRUE
##   Kuldrenett  0.8181818 0.1685393
##   Liivi sibul 0.1818182 0.8314607

barplot(prop.table(table(ounad$diameeter>5, ounad$ounasort)), beside=TRUE)

barplot(prop.table(table(ounad$diameeter>5, ounad$ounasort), 1), beside=TRUE)

barplot(100*prop.table(table(ounad$diameeter>5, ounad$ounasort), 2), beside=TRUE, ylim=c(0, 100), ylab="protsent")