Question

我一直在研究一个PHP页面，该页面计算不同结果的概率，同时从一个由两类人（+和-）组成的较大组中随机选择一个样本组。

例如，考虑到n美国人，它可以计算从美国各地随机选择的一组1000人中有0（或0.15）吸烟者的概率是吸烟者（+）。

在与10000人以下的人群合作时效果非常好，但是当涉及到像1000000这样的人口较多时，它会对所有概率回应0，除非精度（ .之后的数字增加到3000。即使在那种情况下也需要永远。

该代码的工作原理是计算0正数的概率，并对其进行一些计算以得出1正的概率，依此类推。虽然大多数这些概率都是无用的。

我一直在想，如果我能找出一种快速计算非常大因子（如99.999%）的精确（1000000!或更高）值的方法，那就不会有了需要从0开始，计算可以从需要的地方开始，并且精度非常低，以减少所需的时间。

以下是代码：

<html>
<body>
    <form method="get">
        target population:
        <input type="text" name="tpop"><br><br>
        selected population:
        <input type="text" name="selected"><br><br>
        fraction of positives:
        <input type="text" name="fop"><br><br>
        output margin:
        <input type="text" name="margin"><br><br>
        precision:
        <input type="text" name="precision"><br><br>
        <input type="submit">
    </form>
</body>
</html>
<?php
set_time_limit(0);
if (isset($_GET["precision"],$_GET["tpop"],$_GET["selected"],$_GET["fop"],$_GET["margin"])&&$_GET["precision"]!=''&&$_GET["tpop"]!=''&&$_GET["selected"]!=''&&$_GET["fop"]!=''&&$_GET["margin"]!=''&&$_GET["tpop"]>=$_GET["selected"]&&$_GET["fop"]>=0&&$_GET["fop"]<=1){
    $tpop=$_GET["tpop"];
    $selected=$_GET["selected"];
    $fop=$_GET["fop"];
    $margin=$_GET["margin"];
    $ioioio=0;
    $precision=$_GET["precision"];

    $minor=($selected*$fop)-$margin;
    $maxor=$minor+(2*$margin);
    $popopo=bcmul($tpop,$fop);
    echo '<br><br>min is'.$minor.'max is'.$maxor.'<br><br>';

    $mmm=bcsub($tpop,$selected,$precision);
    $rea=bcsub($mmm,1,$precision);
    $fops=bcmul($tpop,$fop);
    $trss=bcsub($tpop,$fops,$precision);
    $trss=bcsub($trss,$selected,$precision);
    $trss=bcadd($trss,1,$precision);
    while($rea>=$trss){
    $mmm=bcmul($mmm,$rea,$precision);
    $rea=bcsub($rea,1,$precision);
}

$nnn=$tpop;
$sfg=bcsub($nnn,1,$precision);
$ugt=bcmul($tpop,$fop,$precision);
$uyt=bcsub($tpop,$ugt);
$uyt=bcadd($uyt,1,$precision);
while($sfg>=$uyt){
    $nnn=bcmul($nnn,$sfg,$precision); 
    $sfg=bcsub($sfg,1,$precision);
}

$zero=bcdiv($mmm,$nnn,$precision);

echo '0==>'.$zero.'<br><br>';
$a=$selected;
$b=($tpop-($tpop*$fop)-$selected+1);
$c=1;
$d=($tpop*$fop);
$i=1;

$origzero=$zero;
$save=$zero;

while($i<=$selected){

    if($d<=0){
        echo $i.'==>impossible<br><br>';
        $a--;
        $b++;
        $c++;
        $d--;
        $i++;
    }else{
        $zero=bcmul($zero,$a,$precision);
        $zero=bcmul($zero,$d,$precision);
        $zero=bcdiv($zero,$b,$precision);
        $zero=bcdiv($zero,$c,$precision);
        if($i>=$minor){
            if($i<=$maxor){
                if($i<=$popopo){
                    $ioioio=bcadd($ioioio,$zero,$precision);
                    echo 'following value is included in p value<br>';
                    echo $i.'==>'.$zero.'<br><br>';
                }
            }
        }
        $save=bcadd($save,$zero,$precision);

        $a--;
        $b++;
        $c++;
        $d--;
        $i++;
    }
}
echo 'precision==> '.$save.'<br><br>';
$savee = bcsub(1,$save,$precision);
echo '1-precision==> '.$savee.'<br><br>';
if($minor<0||$maxor>$selected){
    echo 'p value==>margin larger than surronding probabilties select an smaller margin to calculate p value';
} elseif($minor>0){
    echo 'p value==>'.$ioioio;
} else{
    $ioioiop=bcadd($ioioio,$origzero,$precision);
    echo 'p value(0included)==> '.$ioioiop;}
}
?>

亲爱的@ shukshin.ivan 非常感谢您的回复，这正是我所寻求的：] 这是一个如何为其他人工作的例子：

$x=950000;
$x =2*$x+1;
$P=pi();
$x =(log(2.0*$P)+log($x/2.0)*$x-$x-(1.0-7.0/(30.0*$x*$x))/(6.0*$x))/2.0;
$x=$x/log(10);
$ex=floor($x);              
$x=pow(10,$x-$ex);
$res=$x.'A';
$res=substr($x,0,6).'E'.$ex;
echo $res;

Answer 1

您可以使用Stirlings approximation。它在大数字上相当精确。意思是阶乘可以计算为近似值

一组其他算法can be found here。

PHP中

1 个答案: