沖ドキ　32ゲーム以内の当選は天国即抜けとみるべきか通常モードでの当選とみるべきか

f:id:b092203w:20190610192050p:plain — RBの3ゲームを天国単発とみるべきか通常モードでの当選とみるべきか

上の画像はデータロボセブン様から引用した某ホールの沖ドキの履歴です。

ここで問題になってくるのが３ゲームで当選したRBで、これを天国単発とみるべきか通常モードでの当選とみるべきかで、天国スルー回数が天国単発なら１回か通常当選なら４回で大きく違ってきます。

そこで今回は設定１で32ゲーム以内にボーナス解除した場合に天国と通常それぞれに滞在している確率を求めました。

確率を求めるときに使用したコードは以下の通りです。

#1ゲーム目にビッグボーナスに当選した時に天国に滞在する確率

rBB1 <- function(){
a <- 1/5.05 #リプレイ
b <- 1/32.13　#角チェリー
c <- 1/128.00　#スイカ
d <- 1/8192.00　#確定役
e <- 1/10922.67　#確定チェリー
f <- 1/32768.00　#中段チェリー
g <- 0.1250 #スイカ告知(6)
h <- 0.0625 #チェリー告知(6)
i <- 0.0183 #RB告知(6)
j <- 0.0427 #BB告知(6)
z <- (1 - b - c - d - e - f)*j
g <- 0.0183 #スイカ告知(1)
h <- 0.0092 #チェリー告知(1)
i <- 0.0007 #RB告知(1)
j <- 0.0009 #BB告知(1)
w <- (1 - b - c - d - e - f)*j
cat(z/(z+w),"\n")
}

#1ゲーム目にレギュラーボーナスに当選した時に天国に滞在する確率

rRB1 <- function(){
a <- 1/5.05 #リプレイ
b <- 1/32.13　#角チェリー
c <- 1/128.00　#スイカ
d <- 1/8192.00　#確定役
e <- 1/10922.67　#確定チェリー
f <- 1/32768.00　#中段チェリー
g <- 0.1250 #スイカ告知
h <- 0.0625 #チェリー告知
i <- 0.0183 #RB告知
j <- 0.0427 #BB告知
z <- (1 - b - c - d - e - f)*i
g <- 0.0183 #スイカ告知(1)
h <- 0.0092 #チェリー告知(1)
i <- 0.0007 #RB告知(1)
j <- 0.0009 #BB告知(1)
w <- (1 - b - c - d - e - f)*i
cat(z/(z+w),"\n")
}

#2ゲーム目から31ゲーム目にビッグボーナスに当選した時に天国に滞在する確率

rBB <- function(){
a <- 1/5.05 #リプレイ
b <- 1/32.13　#角チェリー
c <- 1/128.00　#スイカ
d <- 1/8192.00　#確定役
e <- 1/10922.67　#確定チェリー
f <- 1/32768.00　#中段チェリー
for(k in 2:31){
g <- 0.1250 #スイカ告知(6)
h <- 0.0625 #チェリー告知(6)
i <- 0.0183 #RB告知(6)
j <- 0.0427 #BB告知(6)
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)
n <- c * (1 - 2 * g)
z <- (x + m + n)^(k - 1)*y + (x + m + n)^(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)
n <- c * (1 - g)
w <- (x + m + n)^(k - 1)*y + (x + m + n)^(k - 2)*(y + v + s)
cat(z/(z+w),"\n")
}
}

#2ゲーム目から31ゲーム目にレギュラーボーナスに当選した時に天国に滞在する確率

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　#中段チェリー
for(k in 2:n){
g <- 0.1250 #スイカ次告知
h <- 0.0625 #チェリー次告知
i <- 0.0183 #RB告知
j <- 0.0427 #BB告知
x <- (1 - b - c - d - e - f)*(1 - 2*i-2*j)#ハズレorリプレイ
y <- (1 - b - c - d - e - f)*i
m <- b * (1 - h)
n <- c * (1 - 2 * g)
z <- (x + m + n)^(k - 1)*y + (x + m + n)^(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)
n <- c * (1 - 2 * g)
w <- (x + m + n)^(k - 1)*y + (x + m + n)^(k - 2)*(y + c * g)
cat(z/(z+w),"\n")
}
}

#32ゲーム目にビッグボーナスに当選した時に天国に滞在する確率

rBB32 <- function(){
a <- 1/5.05 #リプレイ
b <- 1/32.13　#角チェリー
c <- 1/128.00　#スイカ
d <- 1/8192.00　#確定役
e <- 1/10922.67　#確定チェリー
f <- 1/32768.00　#中段チェリー
g <- 0.1250 #スイカ告知
h <- 0.0625 #チェリー告知
i <- 0.0183 #RB告知
j <- 0.0427 #BB告知
z <- (1 - b - c - d - e - f)*j
for(k in 2:31){
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)
n <- c * (1 - 2 * g)
cat((x + m + n)^(k - 1)*y + (x + m + n)^(k - 2)*(y + v + s),"\n")
z <- z + (x + m + n)^(k - 1)*y + (x + m + n)^(k - 2)*(y + c * g)
}
z <- 1 - z
cat(z,"\n")
g <- 0.0183 #スイカ次告知(1)
h <- 0.0092 #チェリー次告知(1)
i <- 0.0007 #RB告知(1)
j <- 0.0009 #BB告知(1)
w <- (x + m + n)^31*y + (x + m + n)^30 * (y + v + s)
cat(z/(z+w),"\n")
}

#32ゲーム目にレギュラーボーナスに当選した時に天国に滞在する確率

rRB32 <- function(){
a <- 1/5.05 #リプレイ
b <- 1/32.13　#角チェリー
c <- 1/128.00　#スイカ
d <- 1/8192.00　#確定役
e <- 1/10922.67　#確定チェリー
f <- 1/32768.00　#中段チェリー
g <- 0.1250 #スイカ告知
h <- 0.0625 #チェリー告知
i <- 0.0183 #RB告知
j <- 0.0427 #BB告知
z <- (1 - b - c - d - e - f)*i
for(k in 2:31){
x <- (1 - b - c - d - e - f)*(1 - 2*i-2*j)
y <- (1 - b - c - d - e - f)*i
m <- b * (1 - h)
n <- c * (1 - 2 * g)
cat((x + m + n)^(k - 1)*y + (x + m + n)^(k - 2)*(y + c * g),"\n")
z <- z + (x + m + n)^(k - 1)*y + (x + m + n)^(k - 2)*(y + c * g)
}
z <- 1 - z
cat(z,"\n")
g <- 0.0183 #スイカ次告知(1)
h <- 0.0092 #チェリー次告知(1)
i <- 0.0007 #RB告知(1)
j <- 0.0009 #BB告知(1)
w <- (x + m + n)^31*y + (x + m + n)^30 * (y + c * g)
cat(z/(z+w),"\n")
}