Browse Source

Correct Entropy property

Normalization terms were missing. Corrected as per - https://en.wikipedia.org/wiki/Dirichlet_distribution#Entropy
pull/406/head
Rohit Pandey 10 years ago
committed by GitHub
parent
commit
185cf61fb8
  1. 4
      src/Numerics/Distributions/Dirichlet.cs

4
src/Numerics/Distributions/Dirichlet.cs

@ -237,8 +237,8 @@ namespace MathNet.Numerics.Distributions
{
get
{
var num = _alpha.Sum(t => (t - 1)*SpecialFunctions.DiGamma(t));
return SpecialFunctions.GammaLn(AlphaSum) + ((AlphaSum - Dimension)*SpecialFunctions.DiGamma(AlphaSum)) - num;
var num = _alpha.Sum(t => (t - 1)*SpecialFunctions.DiGamma(t) - SpecialFunctions.GammaLn(t));
return -SpecialFunctions.GammaLn(AlphaSum) + ((AlphaSum - Dimension)*SpecialFunctions.DiGamma(AlphaSum)) - num;
}
}

Loading…
Cancel
Save