scala - 为什么 Breeze 二项分布在应该返回 1.0 时返回 NaN？

Question

下面给出的是 NaN，但应该给出 1.0，为什么？

breeze.stats.distributions.Binomial(10, 1.0).probabilityOf(10)

同样对于Binomial(1, 0.0).probabilityOf(0)

Answer 1

从微风源。

case class Binomial(n: Int, p: Double)(implicit rand: RandBasis=Rand) extends DiscreteDistr[Int] with Moments[Double, Double] {
  require(n > 0, "n must be positive!");
  require(p >= 0.0, "p must be non-negative!");
  def probabilityOf(k: Int) = exp(logProbabilityOf(k));

  override def toString() = "Binomial(" + n + ", " + p + ")";

  override def logProbabilityOf(k: Int) = {
    require(n >= k);
    require(k >= 0);
    lgamma(n+1) - lgamma(k+1) - lgamma(n-k+1) + k * log(p) + (n-k) * log(1-p)
  }

0 的日志未定义。在您的第一个示例中，您将 p 设为 1，因此log(1-p) = log(0) = undefined在您的第二个示例中，最后一个，log(p) = log(0) = undefined

由于这些概率都应该是1，我认为这是微风源代码中的一个错误