public void ValidateBetaSpecialCaseDensityLn(double x) { var d = new Dirichlet(new[] { 0.1, 0.3 }); var beta = new Beta(0.1, 0.3); AssertHelpers.AlmostEqualRelative(d.DensityLn(new[] { x }), beta.DensityLn(x), 10); }
public void ValidateDensityLn(double[] x) { var d = new Dirichlet(new[] { 0.1, 0.3, 0.5, 0.8 }); AssertHelpers.AlmostEqualRelative(d.DensityLn(x), Math.Log(d.Density(x)), 12); }