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