beta reg

Dependents

this module uses 5 subunits Log_GammaGammaMachine_EpsilonPoly_EvalStirling

Interface

Overview

These function compute the regularized incomplete beta integral of the arguments, evaluated from zero to x. This is occasionally also called the beta function ratio, given that it returns the probability that a random variable drawn from a beta distribution with parameters a and b will be less that or equal to x.

This ratio (or regularization) function is calculated by dividing the incomplete beta integral by the complete beta integral, i.e.

If the incomplete beta integral is defined by: and the complete beta integral is defined by: where is the Maths/Special/Gamma/Gamma function.

Therefore the regularized beta integral is: where .

All the functions are only valid in the range 0 <= x <= 1. In this implementation a and b are restricted to positive values.

The integral from x to 1 may be obtained by the symmetry relation, i.e.

1 - beta_reg( x, a, b )  =  beta_reg( 1-x, a, b ).

Example:

The following example evaluates the upper regularized incomplete beta integral for 10 points equally spaced in the interval from 0 to 1.

#include <stdio.h>
#include <codecogs/maths/special/gamma/beta_reg.h>
void main()
{
  for (double x=0; x<1; x+=0.1)
  {
    double y = Maths::Special::Gamma::beta_reg(x, 2, 1, true);
    printf("beta_reg(2, %.1lf, 1, true) = %lf\n", x,y);
  }
}

Output:

beta_reg(2, 0.0, 1, true) = 1.000000
beta_reg(2, 0.1, 1, true) = 0.990000
beta_reg(2, 0.2, 1, true) = 0.960000
beta_reg(2, 0.3, 1, true) = 0.910000
beta_reg(2, 0.4, 1, true) = 0.840000
beta_reg(2, 0.5, 1, true) = 0.750000
beta_reg(2, 0.6, 1, true) = 0.640000
beta_reg(2, 0.7, 1, true) = 0.510000
beta_reg(2, 0.8, 1, true) = 0.360000
beta_reg(2, 0.9, 1, true) = 0.190000
beta_reg(2, 1.0, 1, true) = 0.000000

Fraction Exp

Computes using continued fraction expansion that main term need to compute the incomplete Beta integral at the extremes.

Methodolgy

Parameters

x	upper limit of integration. (0<x<1)
a	first shape argument, must be positive.
b	second shape argument, must be positive.
inverse	computes the fraction expansion

Source Code

Source code is available when you agree to a GP Licence or buy a Commercial Licence.

Not a member, then Register with CodeCogs. Already a Member, then Login.

BetaLower Reg Pow

Power series for the regularized incomplete beta integral,, Use when is little and x is not too close to 1.

Parameters

x	upper limit of integration (0<x<1)
a	first shape argument, must be positive
b	second shape argument, must be positive

Authors

Stephen L.Moshier. Copyright 1984, 1987, 1989, 1992, 2000

Source Code

Source code is available when you agree to a GP Licence or buy a Commercial Licence.

Not a member, then Register with CodeCogs. Already a Member, then Login.

Beta Reg

Computes the regularize beta integral using a variety of method depending on the input values:

• When b*x<=1 and x<0.95 then the function is calculated using the power series expansion defined in betaLower_reg_pow.
• If x is greater than the mean of the beta function, i.e.
  

  (12)
  Then the values of a and b are reversed and the calculation proceeds as before only with x=1-x.
• A second attempt is then made to use the power series expansion within the range b*x<=1 and x<0.95.
• Failing this a fractional expansion is used within the tails of this function.

Parameters

x	the value at which to evaluate the function, must be in range 0..1
a	the 1st number of degrees of freedom, must be strictly positive
b	the 2nd number of degrees of freedom, must be strictly positive
upper	Default value = false

Authors

Stephen L.Moshier. Copyright 1984, 1987, 1989, 1992, 2000

Source Code

Source code is available when you agree to a GP Licence or buy a Commercial Licence.

Not a member, then Register with CodeCogs. Already a Member, then Login.