张宏良博客?'':(张洪量?'':张虹亮'blog)

来源：https://github.com/jkromski/wilson-score-interval

<?php
	class WilsonScoreInterval{
		var $positive,$negative,$confidence,$n,$p,$z;
		function setVar($positive,$negative,$confidence=0.95){//2个参数分别为总评数和差评数
			$this->positive=intval($positive);
			$this->negative=intval($negative);
			$this->confidence=floatval($confidence);
			$this->n=$this->positive+$this->negative;
			$this->z=$this->pnorm(1-(1-$this->confidence)/2);
			$this->p=$this->n>0?1.0*$this->positive/$this->n:null;
		}
		function getLowerBound(){return $this->n>0?$this->multiplier()*($this->innerLeft()-$this->innerRight()):0;}
		function getUpperBound(){return $this->n>0?$this->multiplier()*($this->innerLeft()+$this->innerRight()):0;}
		function score(){return $this->getLowerBound();}
		function getInterval(){return array(0=>$this->getLowerBound(),1=>$this->getUpperBound());}
		function pnorm($qn){
			$b=array(
				1.570796288,0.03706987906,-0.8364353589e-3,-0.2250947176e-3,0.6841218299e-5,0.5824238515e-5,
				-0.104527497e-5,0.8360937017e-7,-0.3231081277e-8,0.3657763036e-10,0.6936233982e-12
			);
			if($qn<0.0 || $qn>1.0 || $qn==0.5) {return 0.0;}
			$w1=$qn>0.5?1.0-$qn:$qn;
			$w3=-log(4.0*$w1*(1.0-$w1));
			$w1=$b[0];
			for($i=1;$i<=10;$i++){$w1+=$b[$i]*pow($w3,$i);}
			return $qn>0.5?sqrt($w1*$w3):-sqrt($w1*$w3);
		}
		function multiplier(){return 1/(1+(pow($this->z,2)/$this->n));}
		function innerLeft(){return $this->p+(pow($this->z,2)/(2*$this->n));}
		function innerRight(){return $this->z*sqrt(($this->p*(1-$this->p)/$this->n)+(pow($this->z,2)/(4*pow($this->n,2))));}
	}
?>

使用方法：

<?php
	require 'wilsonscore.class.php';
	$WilsonScoreInterval=new WilsonScoreInterval();

	$WilsonScoreInterval->setVar(10,1);//2个参数分别为总评数和差评数
	$wilsonScore=substr($WilsonScoreInterval->score(),0,5);
	die($wilsonScore.'x');
?>

顺道贴上GitHub上的原味Wilson score代码类

<?php
namespace ScoreInterval;
class WilsonScoreInterval{
    protected $positive; // Count of Positive results
    protected $negative; // Count of Negative results
    protected $confidence; // Level of Confidence
    protected $n; // Count of Total results
    protected $p; // Positive results out of total results
    protected $z; // Quantile of the standard normal distribution
    /**
     * @param mixed $positive
     * @param mixed $negative
     * @param mixed $confidence
     */
    public function __construct($positive, $negative, $confidence = 0.95){
        $this->positive = intval($positive);
        $this->negative = intval($negative);
        $this->confidence = floatval($confidence);
        $this->n = $this->positive + $this->negative;
        $this->z = $this->pnorm(1 - (1 - $this->confidence) / 2);
        $this->p = $this->n > 0 ? 1.0 * $this->positive / $this->n : null;
    }
    /**
     * @return float|int
     */
    public function getLowerBound(){ return $this->n > 0 ? $this->multiplier() * ($this->innerLeft() - $this->innerRight()) : 0; }
    /**
     * @return float|int
     */
    public function getUpperBound(){ return $this->n > 0 ? $this->multiplier() * ($this->innerLeft() + $this->innerRight()) : 0; }
    /**
     * @return float|int
     */
    public function score(){ return $this->getLowerBound(); }
    /**
     * @return array
     */
    public function getInterval(){ return array(0 => $this->getLowerBound(), 1 => $this->getUpperBound()); }
    /**
     * @link https://github.com/abscondment/statistics2/blob/master/lib/statistics2/base.rb#L89
     * @param $qn
     * @return float
     */
    protected function pnorm($qn){
        $b = array(
            1.570796288, 0.03706987906, -0.8364353589e-3, -0.2250947176e-3,
            0.6841218299e-5, 0.5824238515e-5, -0.104527497e-5, 0.8360937017e-7,
            -0.3231081277e-8, 0.3657763036e-10, 0.6936233982e-12
        );
        if ($qn < 0.0 || $qn > 1.0 || $qn == 0.5) {
            return 0.0;
        }
        $w1 = $qn > 0.5 ? 1.0 - $qn : $qn;
        $w3 = - log(4.0 * $w1 * (1.0 - $w1));
        $w1 = $b[0];
        for ($i = 1; $i <= 10; $i++) {
            $w1 += $b[$i] * pow($w3, $i);
        }
        return $qn > 0.5 ? sqrt($w1 * $w3) : -sqrt($w1 * $w3);
    }
    /**
     * @return float
     */
    private function multiplier(){ return 1 / (1 + (pow($this->z, 2) / $this->n)); }
    /**
     * @return float
     */
    private function innerLeft(){ return $this->p + (pow($this->z, 2) / (2 * $this->n)); }
    /**
     * @return float
     */
    private function innerRight(){ return $this->z * sqrt(($this->p * (1 - $this->p) / $this->n) + (pow($this->z, 2) / (4 * pow($this->n, 2)))); }
}