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