@@ 1374-1388 (lines=15) @@ | ||
1371 | $PreviousValue = 1; |
|
1372 | $Done = false; |
|
1373 | $k = $m + 1; |
|
1374 | while ((!$Done) && ($k <= $trials)) { |
|
1375 | $CurrentValue = $PreviousValue * ($trials - $k + 1) * $probability / ($k * (1 - $probability)); |
|
1376 | $TotalUnscaledProbability += $CurrentValue; |
|
1377 | if ($k == $Guess) { |
|
1378 | $UnscaledPGuess += $CurrentValue; |
|
1379 | } |
|
1380 | if ($k <= $Guess) { |
|
1381 | $UnscaledCumPGuess += $CurrentValue; |
|
1382 | } |
|
1383 | if ($CurrentValue <= $EssentiallyZero) { |
|
1384 | $Done = true; |
|
1385 | } |
|
1386 | $PreviousValue = $CurrentValue; |
|
1387 | ++$k; |
|
1388 | } |
|
1389 | ||
1390 | $PreviousValue = 1; |
|
1391 | $Done = false; |
|
@@ 1393-1407 (lines=15) @@ | ||
1390 | $PreviousValue = 1; |
|
1391 | $Done = false; |
|
1392 | $k = $m - 1; |
|
1393 | while ((!$Done) && ($k >= 0)) { |
|
1394 | $CurrentValue = $PreviousValue * $k + 1 * (1 - $probability) / (($trials - $k) * $probability); |
|
1395 | $TotalUnscaledProbability += $CurrentValue; |
|
1396 | if ($k == $Guess) { |
|
1397 | $UnscaledPGuess += $CurrentValue; |
|
1398 | } |
|
1399 | if ($k <= $Guess) { |
|
1400 | $UnscaledCumPGuess += $CurrentValue; |
|
1401 | } |
|
1402 | if ($CurrentValue <= $EssentiallyZero) { |
|
1403 | $Done = true; |
|
1404 | } |
|
1405 | $PreviousValue = $CurrentValue; |
|
1406 | --$k; |
|
1407 | } |
|
1408 | ||
1409 | $PGuess = $UnscaledPGuess / $TotalUnscaledProbability; |
|
1410 | $CumPGuess = $UnscaledCumPGuess / $TotalUnscaledProbability; |