Verjetnost, da v talonu dobiš enega, dva, tri ali štiri “kralje”.
kartopivec <- function(karte, talon, iscemKarte, stPonovitev) {
# karte = katere karte imamo na razpolago talon = koliko kart potegnemo iz
# kupčka iscemKarte = točno katere karte iščem? stPonovitev = kolikokrat
# ponovimo simulacijo vleka 6 kart
# Vlečemo karte brez ponovitve in izračunamo koliko 'kraljev'' imamo. To
# ponovimo stPonovitev-krat.
simulacija <- replicate(stPonovitev, {
poteg <- sample(karte, size = talon, replace = FALSE)
sum(iscemKarte %in% poteg) # rezultat ene simulacije je število kraljev v talonu
})
# Izračunamo sorazmeren delež števila 'kraljev' v talonu glede na število
# simulacij.
table(simulacija)/stPonovitev
}
Imamo navaden kup 54 kart. Potegnemo 6 kart. Kolikokrat bomo videli 1, 2, 3 ali 4 “kralje”? Poskus je zasnovan na 100000 ponovitvah.
kartopivec(karte = 1:54, talon = 6, iscemKarte = c(1, 2, 3, 4), stPonovitev = 1e+05)
## simulacija
## 0 1 2 3 4
## 0.61647 0.32723 0.05327 0.00298 0.00005
Enako kot zgoraj, le, da smo 4 karte dali na stran. Predpostavka je, da te 4 karte niso “kralji”.
kartopivec(karte = 1:50, talon = 6, iscemKarte = c(1, 2, 3, 4), stPonovitev = 1e+05)
## simulacija
## 0 1 2 3 4
## 0.59085 0.34286 0.06240 0.00382 0.00007