@@ 1424-1438 (lines=15) @@ | ||
1421 | $PreviousValue = 1; |
|
1422 | $Done = false; |
|
1423 | $k = $m + 1; |
|
1424 | while ((!$Done) && ($k <= $trials)) { |
|
1425 | $CurrentValue = $PreviousValue * ($trials - $k + 1) * $probability / ($k * (1 - $probability)); |
|
1426 | $TotalUnscaledProbability += $CurrentValue; |
|
1427 | if ($k == $Guess) { |
|
1428 | $UnscaledPGuess += $CurrentValue; |
|
1429 | } |
|
1430 | if ($k <= $Guess) { |
|
1431 | $UnscaledCumPGuess += $CurrentValue; |
|
1432 | } |
|
1433 | if ($CurrentValue <= $EssentiallyZero) { |
|
1434 | $Done = true; |
|
1435 | } |
|
1436 | $PreviousValue = $CurrentValue; |
|
1437 | ++$k; |
|
1438 | } |
|
1439 | ||
1440 | $PreviousValue = 1; |
|
1441 | $Done = false; |
|
@@ 1443-1457 (lines=15) @@ | ||
1440 | $PreviousValue = 1; |
|
1441 | $Done = false; |
|
1442 | $k = $m - 1; |
|
1443 | while ((!$Done) && ($k >= 0)) { |
|
1444 | $CurrentValue = $PreviousValue * $k + 1 * (1 - $probability) / (($trials - $k) * $probability); |
|
1445 | $TotalUnscaledProbability += $CurrentValue; |
|
1446 | if ($k == $Guess) { |
|
1447 | $UnscaledPGuess += $CurrentValue; |
|
1448 | } |
|
1449 | if ($k <= $Guess) { |
|
1450 | $UnscaledCumPGuess += $CurrentValue; |
|
1451 | } |
|
1452 | if ($CurrentValue <= $EssentiallyZero) { |
|
1453 | $Done = true; |
|
1454 | } |
|
1455 | $PreviousValue = $CurrentValue; |
|
1456 | --$k; |
|
1457 | } |
|
1458 | ||
1459 | $PGuess = $UnscaledPGuess / $TotalUnscaledProbability; |
|
1460 | $CumPGuess = $UnscaledCumPGuess / $TotalUnscaledProbability; |