Example usage for org.apache.commons.math3.special Beta regularizedBeta

Introduction

In this page you can find the example usage for org.apache.commons.math3.special Beta regularizedBeta.

Prototype

public static double regularizedBeta(double x, double a, double b) 

Source Link

Document

Returns the <a href="http://mathworld.wolfram.com/RegularizedBetaFunction.html"> regularized beta function</a> I(x, a, b).

Usage

From source file:de.tudarmstadt.ukp.experiments.argumentation.sequence.significance.LiddellsExactTest.java

public static double pValue(int b, int d) {
    double r = Math.max(b, d);
    double s = Math.min(b, d);

    double f = r / (s + 1);

    double mm = 2 * (s + 1);
    double nn = 2 * r;

    double betaPar1 = 1.0 / (1 + mm * f / nn);
    double betaPar2 = nn / 2;
    double betaPar3 = mm / 2;

    double beta = Beta.regularizedBeta(betaPar1, betaPar2, betaPar3);

    double pValue = 2 * (1 - (1 - beta));

    //        if (pValue < 0.1) {
    //            System.out.println(r);
    //            System.out.println(s);
    //        }//from  w  w  w.java 2  s .c  o m

    return pValue;
}

---

From source file:beast.math.distributions.NegativeBinomialDistribution.java

public double cdf(double x) {
    double r = -1 * (mean * mean) / (mean - stdev * stdev);
    double p = mean / (stdev * stdev);
    return Beta.regularizedBeta(p, r, x + 1);
}

---

From source file:experiment.PascalDistribution_bug.java

/** {@inheritDoc} */
public double cumulativeProbability(int x) {
    double ret;/* w w  w .  ja  v  a  2  s.  c  om*/
    if (x < 0) {
        ret = 0.0;
    } else {
        ret = Beta.regularizedBeta(probabilityOfSuccess, numberOfSuccesses, x + 1.0);
    }
    return ret;
}

---

From source file:beast.math.distributions.BetaDistribution.java

public double cdf(double x) {
    if (x <= 0) {
        return 0;
    } else if (x >= 1) {
        return 1;
    } else {//  w w w .j a  v  a 2 s  .  co  m
        return Beta.regularizedBeta(x, alpha, beta);
    }

}

---

From source file:beast.math.distributions.BetaDistribution.java

public double cumulativeProbability(double x) {
    if (x <= 0) {
        return 0;
    } else if (x >= 1) {
        return 1;
    } else {//  w ww.j av a2 s  . c om
        return Beta.regularizedBeta(x, alpha, beta);
    }
}

---

From source file:org.nd4j.linalg.api.rng.distribution.impl.BinomialDistribution.java

/**
 * {@inheritDoc}/*from w  w w .j  av a2  s.  co m*/
 */
public double cumulativeProbability(int x) {

    double ret;
    if (x < 0) {
        ret = 0.0;
    } else if (x >= numberOfTrials) {
        ret = 1.0;
    } else {
        ret = 1.0 - Beta.regularizedBeta(probabilityOfSuccess, x + 1.0, numberOfTrials - x);
    }
    return ret;
}

---

From source file:org.nd4j.linalg.api.rng.distribution.impl.BinomialDistribution.java

@Override
public double cumulativeProbability(double x) {

    double ret;//from w w  w. j  a v a2 s  . co m
    if (x < 0) {
        ret = 0.0D;
    } else if (x >= this.numberOfTrials) {
        ret = 1.0D;
    } else {
        ret = 1.0D - Beta.regularizedBeta(this.probabilityOfSuccess, x + 1.0D, (this.numberOfTrials - x));
    }

    return ret;
}

---

From source file:org.nd4j.linalg.jcublas.rng.distribution.BinomialDistribution.java

@Override
public double cumulativeProbability(double x) {
    double ret;/*  ww w. ja v a 2 s  . c  o m*/
    if (x < 0) {
        ret = 0.0;
    } else if (x >= numberOfTrials) {
        ret = 1.0;
    } else {
        ret = 1.0 - Beta.regularizedBeta(probabilityOfSuccess, x + 1.0, numberOfTrials - x);
    }
    return ret;
}

---

From source file:statalign.utils.BetaDistribution.java

/** {@inheritDoc} */
@Override/*from  w  w w. j a va 2s. com*/
public double cumulativeProbability(double x) {
    if (x <= 0) {
        return 0;
    } else if (x >= 1) {
        return 1;
    } else {
        return Beta.regularizedBeta(x, alpha, beta);
    }
}

---