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=