沖ドキ　設定１で３３ゲーム目から引き戻しモード滞在率をシミュレーションで算出しました。

前回の記事で触れた通り、今回は設定１で引き戻しモードに滞在する確率を調べます。

いつもはコードと結果だけ載せるのですが、それだと勉強にならないので、今回は確率を算出するうえで使っている数学的な道具を下にのせます。

f:id:b092203w:20190622185957p:plain — 算出に当たってベイズの定理を使用します

確率を算出するうえで注意しなければならないのは上記のように各モードに滞在しているときにボーナスに当選する確率にモード滞在率を掛けなければいけないことです。

３３ゲームを抜けると、モードは通常Aか通常Bか引き戻しのいずれかなので、

1geki.jp

上のサイトによると天国抜けした直後に通常A、通常B、引き戻し移行する確率はそれぞれ13.28%,3.91%,7.81%なので、３３ゲーム目以降に各モードに滞在している確率は以下の画像のように算出します。

f:id:b092203w:20190622185840p:plain — 各モードに滞在している確率

以上の結果を使って３３ゲーム目以降から引き戻しモードに滞在している確率をシミュレーションで算出します。コードは以下の通りです。

rRB <- function(n){
a <- 1/5.05 #リプレイ
b <- 1/32.13　#角チェリー
c <- 1/128.00　#スイカ
d <- 1/8192.00　#確定役
e <- 1/10922.67　#確定チェリー
f <- 1/32768.00　#中段チェリー

TA <- 0.1328

TB <- 0.0391

TH <- 0.0781

p <- TH/(TA + TB + TH)

q <- (TA + TB)/(TA + TB + TH)
for(k in 2:n){
g <- 0.0458 #スイカ告知(1)
h <- 0.0229 #チェリー告知(1)
i <- 0.0018 #RB告知(1)
j <- 0.0024 #BB告知(1)
x <- (1 - b - c - d - e - f)*(1 - 2*i-2*j)
y <- (1 - b - c - d - e - f)*i
m <- b * (1 - h)
l <- c * (1 - g)
z <- (x + m + l)^(k - 1)*y + (x + m + l)^(k - 2)*(y + c * g)
g <- 0.0183 #スイカ次告知(1)
h <- 0.0092 #チェリー次告知(1)
i <- 0.0007 #RB告知(1)
j <- 0.0009 #BB告知(1)
x <- (1 - b - c - d - e - f)*(1 - 2*i-2*j)
y <- (1 - b - c - d - e - f)*i
m <- b * (1 - h)
l <- c * (1 - 2 * g)
w <- (x + m + l)^(k - 1)*y + (x + m + l)^(k - 2)*(y + c * g)
cat(z*p/(z*p+w*q),"\n")
}
}

rBB <- function(n){
a <- 1/5.05 #リプレイ
b <- 1/32.13　#角チェリー
c <- 1/128.00　#スイカ
d <- 1/8192.00　#確定役
e <- 1/10922.67　#確定チェリー
f <- 1/32768.00　#中段チェリー

TA <- 0.1328

TB <- 0.0391

TH <- 0.0781

p <- TH/(TA + TB + TH)

q <- (TA + TB)/(TA + TB + TH)
for(k in 2:n){
g <- 0.0458 #スイカ告知(1)
h <- 0.0229 #チェリー告知(1)
i <- 0.0018 #RB告知(1)
j <- 0.0024 #BB告知(1)
x <- (1 - b - c - d - e - f)*(1 - 2*i-2*j)
y <- (1 - b - c - d - e - f)*j
v <- b * h
s <- c * g
m <- b * (1 - h)
l <- c * (1 - 2 * g)
z <- (x + m + l)^(k - 1)*y + (x + m + l)^(k - 2)*(y + v + s)
g <- 0.0183 #スイカ告知(1)
h <- 0.0092 #チェリー告知(1)
i <- 0.0007 #RB告知(1)
j <- 0.0009 #BB告知(1)
x <- (1 - b - c - d - e - f)*(1 - 2*i-2*j)
y <- (1 - b - c - d - e - f)*i
v <- b * h
s <- c * g
m <- b * (1 - h)
l<- c * (1 - g)
w <- (x + m + l)^(k - 1)*y + (x + m + l)^(k - 2)*(y + v + s)
cat(z*p/(z*p+w*q),"\n")
}
}