# Time-stamp: <2012-09-07 17:54 Petri Koistinen>

# Arvioidaan ehdollista todennäköisyyttä simuloiduissa nopanheitossa

# 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)

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

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