BetaDiscreteDistribution.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: The Bio++ Development Group
//
// SPDX-License-Identifier: CECILL-2.1

#ifndef BPP_NUMERIC_PROB_BETADISCRETEDISTRIBUTION_H
#define BPP_NUMERIC_PROB_BETADISCRETEDISTRIBUTION_H


#include "../Random/RandomTools.h"
#include "AbstractDiscreteDistribution.h"

namespace bpp
{
class BetaDiscreteDistribution :
  public AbstractDiscreteDistribution
{
private:
  double alpha_, beta_;

  /* for computation economy */

  double diffln_;

public:
  BetaDiscreteDistribution(size_t n, double alpha = 1, double beta = 1, short discretization = DISCRETIZATION_EQUAL_PROB_WHEN_POSSIBLE);

  BetaDiscreteDistribution(const BetaDiscreteDistribution&) = default;

  BetaDiscreteDistribution& operator=(const BetaDiscreteDistribution&) = default;

  BetaDiscreteDistribution* clone() const override { return new BetaDiscreteDistribution(*this); }

public:
  std::string getName() const override { return "Beta";}

  void fireParameterChanged(const ParameterList& parameters) override;

  double randC() const override
  {
    double x = RandomTools::randBeta(getParameterValue("alpha"),
          getParameterValue("beta"));
    while (!intMinMax_->isCorrect(x))
      x = RandomTools::randBeta(getParameterValue("alpha"),
            getParameterValue("beta"));
    return x;
  }

  double qProb(double x) const override;

  double pProb(double x) const override;

  double Expectation(double a) const override;
};
} // end of namespace bpp.
#endif // BPP_NUMERIC_PROB_BETADISCRETEDISTRIBUTION_H