Example 2.14

A fair coin is repeatedly tossed. Estimate the probability that you observe Heads for the third time on the 10th toss.

Answer

event <- replicate(10000, {
  coinToss <- sample(c("H", "T"), 10, replace = TRUE)
  coinToss
  sum(coinToss == "H") # test
  sum(coinToss == "H") == 3
  coinToss[1:9]
  sum(coinToss[1:9] == "H") == 2
  sum(coinToss == "H") == 3 & sum(coinToss[1:9] == "H") == 2
})
mean(event)
[1] 0.0351
LS0tDQp0aXRsZTogIlIgTm90ZWJvb2sgRm9yIGV4YW1wbGUgMi4xNCINCmF1dGhvcjogIkJyYW5kb24gQm95ZWEiDQpkYXRlOiAiOS8xMC8yMDIxIg0Kb3V0cHV0OiBodG1sX25vdGVib29rDQotLS0NCiMjIyBFeGFtcGxlIDIuMTQNCkEgZmFpciBjb2luIGlzIHJlcGVhdGVkbHkgdG9zc2VkLiBFc3RpbWF0ZSB0aGUgcHJvYmFiaWxpdHkgdGhhdCB5b3Ugb2JzZXJ2ZSBIZWFkcyBmb3IgdGhlIHRoaXJkIHRpbWUgb24gdGhlIDEwdGggdG9zcy4NCg0KIyMjIEFuc3dlcg0KDQpgYGB7cn0NCmV2ZW50IDwtIHJlcGxpY2F0ZSgxMDAwMCwgew0KICBjb2luVG9zcyA8LSBzYW1wbGUoYygiSCIsICJUIiksIDEwLCByZXBsYWNlID0gVFJVRSkNCiAgY29pblRvc3MNCiAgc3VtKGNvaW5Ub3NzID09ICJIIikgDQogIHN1bShjb2luVG9zcyA9PSAiSCIpID09IDMNCiAgY29pblRvc3NbMTo5XQ0KICBzdW0oY29pblRvc3NbMTo5XSA9PSAiSCIpID09IDINCiAgc3VtKGNvaW5Ub3NzID09ICJIIikgPT0gMyAmIHN1bShjb2luVG9zc1sxOjldID09ICJIIikgPT0gMg0KfSkNCm1lYW4oZXZlbnQpDQpgYGANCg0KDQo=