ロケミツのサイコロの統計その３（４回連続して６がでるとイカサマ？）

６が連続４回でるとイカサマ？

これは、前回の記事で使った関数を流用できる

prenzoku <- function(m, q, r) {
    pren <- numeric(m - r + 1)

    pren[1] <- q^r
    pren[2:(r + 1)] <- (1 - q) * q^r

    for (i in (r + 2):(m - r + 1)) {
        pren[i] <- (1 - q) * q^r * (1 - sum(pren[1:(i - r - 1)]))
    }

    return(sum(pren))  #答えの確率
}

prenzoku(241, 1/6, 4)

## [1] 0.1424

241回サイコロを振る番組が１００番組あると、４回連続６が出ることは １４番組くらいで起こりうるということ。決してまれとはいえない。

ということで、今回も白判定。

とは言え、本当はいくつもの疑問点について計算を繰り返すのはよくない。 ２が４連続とか、3が４連続とか、１，２，１，２，と繰り返すとか、いくらでも計算する ネタを作ることができる。検定を繰り返せば、そのうち５％以下でしか起こらない珍しい ネタみつかるであろうが、不適切である。