Peter Donnelly'nin konuşma esinlenerek TED o belli desen sikke fırlatır bir dizi görünmesinin ne kadar süreceğini anlatılır ki, ben 'htt' iki paterni 'hth' ve Verilen R. aşağıdaki komut dosyasını oluşturan, onu Bu kalıplardan birine basmadan önce ortalama ne kadar sürdüğünü (yani kaç jeton attığını) hesaplar.
coin <- c('h','t')
hit <- function(seq) {
miss <- TRUE
fail <- 3
trp <- sample(coin,3,replace=T)
while (miss) {
if (all(seq == trp)) {
miss <- FALSE
}
else {
trp <- c(trp[2],trp[3],sample(coin,1,T))
fail <- fail + 1
}
}
return(fail)
}
n <- 5000
trials <- data.frame("hth"=rep(NA,n),"htt"=rep(NA,n))
hth <- c('h','t','h')
htt <- c('h','t','t')
set.seed(4321)
for (i in 1:n) {
trials[i,] <- c(hit(hth),hit(htt))
}
summary(trials)
Özet istatistikler aşağıdaki gibidir,
hth htt
Min. : 3.00 Min. : 3.000
1st Qu.: 4.00 1st Qu.: 5.000
Median : 8.00 Median : 7.000
Mean :10.08 Mean : 8.014
3rd Qu.:13.00 3rd Qu.:10.000
Max. :70.00 Max. :42.000
Konuşmada, ortalama madeni para atma sayısının iki kalıp için farklı olacağı; simülasyonumdan görülebileceği gibi. Konuşmayı birkaç kez izlememe rağmen, neden böyle olacağını tam olarak anlamadım. Anladığım kadarıyla 'hth' kendini örtüşüyor ve sezgisel olarak 'htt' den daha kısa sürede 'hth' vurduğunu düşünürdüm, ama durum bu değil. Biri bana bunu açıklarsa çok sevinirim.