public void LogProbabilityMassFunctionTest()
{
MultinomialDistribution dist = new MultinomialDistribution(5, 0.25, 0.25, 0.25, 0.25);
int[] observation = { 1, 1, 1, 2 };
double actual = dist.LogProbabilityMassFunction(observation);
double expected = System.Math.Log(0.058593750);
Assert.AreEqual(expected, actual, 1e-6);
}