Есть JS библиотеки для работы с статистическими рассчетами (мне нужно распределение Пуассона)?

Question

[coderisimo]

Собственно, вопрос в заголовке. Подойдет и просто отдельная реализация данного распределения.

Спасибо

Answer 1

[Roman_Kh]

Есть, например, jStat, в котором как раз есть распределение Пуассона.

Comments

[coderisimo]

Спасибо! Если понадобятся специализированные расчеты обязательно воспользуюсь. А пока я уже выдрал нужный кусок и адаптировал под свои нужды ))
Может кому-то пригодится.

// calculate Poisson Distribution

//первая функция это ОНО 
   function PoissonDistrib(key_var , lambda_var){
   
   var u = PosV(parseFloat(lambda_var));
   var k = PosV(parseInt(key_var));
   
   var result = Fixed(PoissonTerm( u, k ),8,4);
   
   return result*100; 
   
   }



function Fixed( s, wid, dec ) {
   // many combinations of possibilities

   // maybe prepare for upcoming truncate
   var z = 1
   if (dec > 0) {
      z /= Math.pow( 10, dec );
      if (s < -z)  s -= 0.5 * z;
      else
         if (s > z)  s += 0.5 * z;
         else
            s = 0;
      }

   // assure a string
   s = "" + s;

   // chop neg, if any
   var neg = 0;
   if (s.charAt(0) == "-") {
      neg = 2;
      s = s.substring( 1, s.length );
      }

   // chop exponent, if any
   var exp = "";
   var e = s.lastIndexOf( "E" );
   if (e < 0)  e = s.lastIndexOf( "e" );
   if (e > -1) {
      exp = s.substring( e, s.length );
      s = s.substring( 0, e );
      }

   // if dec > 0 assure "."; dp == index of "."
   var dp = s.indexOf( ".", 0 );
   if (dp == -1) {
      dp = s.length;
      if (dec > 0) {
         s += ".";
         dp = s.length - 1;
         }
      }

   // assure leading digit
   if (dp == 0) {
      s = '0' + s;
      dp = 1;
      }

   // not enough dec pl?  add 0's
   while ((dec > 0) && ((s.length - dp - 1) < dec))
      s += "0";

   // too many dec pl?  take a substring
   var places = s.length - dp - 1;
   if (places > dec)
      if (dec == 0)
         s = s.substring( 0, dp );
      else
         s = s.substring( 0, dp + dec + 1 );

   // recover exponent, if any
   s += exp;

   // recover neg, if any
   if (neg > 0)
      s = "-" + s;

   // if not enough width, add spaces IN FRONT
   //    too many places?  tough!
   while (s.length < wid)
      s = " " + s;

   return s
   }

   
   
   
   function DoPoi( aform ) {
   var u = PosV(parseFloat(aform.u.value));
   aform.u.value = Fixed(u,10,4);
   var k = PosV(parseInt(aform.k.value));
   aform.k.value = Fixed(k,8,0);
   aform.tnk.value = Fixed(PoissonTerm( u, k ),8,4);
   var t = PoissonP( u, k );
   aform.puk.value = Fixed(t,8,4);
   aform.quk.value = Fixed(1-t,8,4);
   }

   
function PoissonTerm( u, k ) {
   // by logs
   return  Math.exp( (k * Math.log(u)) - u - LnFact(k) );
   }

   
   function LnFact( x ) {
   // ln(x!) by Stirling's formula
   //   see Knuth I: 111
   if (x <= 1)  x = 1;

   if (x < 12)
      return Math.log( Fact(Math.round(x)) );
   else {
      var invx = 1 / x;
      var invx2 = invx * invx;
      var invx3 = invx2 * invx;
      var invx5 = invx3 * invx2;
      var invx7 = invx5 * invx2;

      var sum = ((x + 0.5) * Math.log(x)) - x;
      sum += Math.log(2*Math.PI) / 2;
      sum += (invx / 12) - (invx3 / 360);
      sum += (invx5 / 1260) - (invx7 / 1680);

      return sum;
      }
   }

function Fact( x ) {
   // x factorial
   var  t=1;
   while (x > 1)
      t *= x--;
   return t;
   }

   
function PosV( x ) {
   if (x < 0)  x = -x;
   return x;
   }