## Time-stamp: <2013-08-28 12:10 Koistinen Petri P M>

##-------------------------------------------------------------
##
## Ehdollinen todennäköisyys simuloiduissa nopanheitossa ---
## esimerkki ehdollisen todennäköisyyden frekvenssitulkinnasta
##
##-------------------------------------------------------------

## Simuloidaan n = 100000 nopanheittoa
n <- 100000
x <- sample(6, n, replace = TRUE)

## Poimitaan ne heitot, joissa silmäluku on parillinen
even <- (x %% 2) == 0
y <- x[even]

## Tällaisia tuloksia saatiin, kun silmäluku oli parillinen:
table(y)

## Arvio tn:lle, että silmäluku on yksi, kun silmäluku on parillinen:

sum(y == 1) / length(y)

## Arvio tn:lle, että silmaluku on kaksi, kun silmäluku on parillinen

sum(y == 2) / length(y)


## Local Variables:
## coding: utf-8
## End: