<<<<<<< HEAD

    gammaA = new Gamma(a, 1);
    gammaB = new Gamma(b, 1);
  }
  /**
   * Returns the probability of value under this distribution.
   */
  public double getProb(List args, Object value) {
    if (!(value instanceof Number)) {
      throw new IllegalArgumentException(
          "Beta CPD defines a distribution over objects"
              + " of class Number, not " + value.getClass() + ".");
    } else {
      double x = ((Number) value).doubleValue();
      return ((Math.pow(x, (a - 1)) * Math.pow((1 - x), (b - 1))) / beta(a, b));
    }
  }

  /**
   * Returns the log of the probability of value under this distribution.
   */
  public double getLogProb(List args, Object value) {
    if (!(value instanceof Number)) {
      throw new IllegalArgumentException(
          "Beta CPD defines a distribution over objects"
              + " of class Number, not " + value.getClass() + ".");
    } else {
      double x = ((Number) value).doubleValue();
      return (((a - 1) * Math.log(x)) + ((b - 1) * Math.log(1 - x)) - Math
          .log(beta(a, b)));
    }
  }

  /**
   * Returns a double sampled according to this distribution. Takes time
   * equivalent to the distrib.Gamma sampling function. (Reference: A Guide To
   * Simulation, 2nd Ed. Bratley, Paul, Bennett L. Fox and Linus E. Schrage.)
   */
  public Object sampleVal(List args) {
    LinkedList l = new LinkedList();
    double y = ((Double) gammaA.sampleVal(l)).doubleValue();
    double z = ((Double) gammaB.sampleVal(l)).doubleValue();
    return new Double(y / (y + z));
  }

  /**
   * Returns the Beta function of reals a and b B(a,b) = Gamma(a)Gamma(b) /
   * Gamma(a+b) Reference: Numerical Recipes in C
   * http://www.library.cornell.edu/nr/cbookcpdf.html
   */
  public static double beta(double a, double b) {
    return ((Gamma.gamma(a) * Gamma.gamma(b)) / Gamma.gamma(a + b));
  }

  public String toString() {
    return getClass().getName();
  }

  private double a;
  private double b;
  private Gamma gammaA;
  private Gamma gammaB;

  /*
   * (non-Javadoc)
   * 
   * @see blog.distrib.CondProbDistrib#setParams(java.util.List)
   */
  @Override
  public void setParams(List

  /**
    gammaA = new Gamma(a, 1);
    gammaB = new Gamma(b, 1);
  }

   * Returns the probability of value under this distribution.
   */
  public double getProb(List args, Object value) {
    if (!(value instanceof Number)) {
      throw new IllegalArgumentException(
          "Beta CPD defines a distribution over objects"
              + " of class Number, not " + value.getClass() + ".");
    } else {
      double x = ((Number) value).doubleValue();
      return ((Math.pow(x, (a - 1)) * Math.pow((1 - x), (b - 1))) / beta(a, b));
    }
  }

  /**
   * Returns the log of the probability of value under this distribution.
   */
  public double getLogProb(List args, Object value) {
    if (!(value instanceof Number)) {
      throw new IllegalArgumentException(
          "Beta CPD defines a distribution over objects"
              + " of class Number, not " + value.getClass() + ".");
    } else {
      double x = ((Number) value).doubleValue();
      double t1 = 0;
      double t2 = 0;
      if (a != 1) {
        t1 = (a - 1) * Math.log(x);
      }
      if (b != 1) {
        t2 = (b - 1) * Math.log(1 - x);
      }
      return t1 + t2 - Math.log(beta(a, b));
    }
  }

  /**
   * Returns a double sampled according to this distribution. Takes time
   * equivalent to the distrib.Gamma sampling function. (Reference: A Guide To
   * Simulation, 2nd Ed. Bratley, Paul, Bennett L. Fox and Linus E. Schrage.)
   */
  public Object sampleVal(List args, Type childType) {
    LinkedList l = new LinkedList();
    double y = ((Double) gammaA.sampleVal(l, childType)).doubleValue();
    double z = ((Double) gammaB.sampleVal(l, childType)).doubleValue();
    return new Double(y / (y + z));
  }

  /**
   * Returns the Beta function of reals a and b B(a,b) = Gamma(a)Gamma(b) /
   * Gamma(a+b) Reference: Numerical Recipes in C
   * http://www.library.cornell.edu/nr/cbookcpdf.html
   */
  public static double beta(double a, double b) {
    return ((Gamma.gamma(a) * Gamma.gamma(b)) / Gamma.gamma(a + b));
  }

  public String toString() {
    return getClass().getName();
  }

  private double a;
  private double b;
  private Gamma gammaA;
  private Gamma gammaB;

  /*
   * (non-Javadoc)
   * 
   * @see blog.distrib.CondProbDistrib#setParams(java.util.List)
   */
  @Override
  public void setParams(List