public void ChiSquareDistribution()
{
int B = 8;
Random rng = new Random(0);
//BinomialDistribution d = new BinomialDistribution(0.25, 6);
//DiscreteUniformDistribution d = new DiscreteUniformDistribution(0, 4);
//BernoulliDistribution d = new BernoulliDistribution(0.25);
DiscreteDistribution d = new PoissonDistribution(2.0);
Sample s = new Sample();
ChiSquaredDistribution nullDistribution = null;
for (int i = 0; i < 512; i++) {
Histogram h = new Histogram(B);
for (int j = 0; j < 1024; j++) {
int k = d.GetRandomValue(rng);
h.Add(k);
//if (k < h.Count) h[k].Increment();
//h[d.GetRandomValue(rng)].Increment();
}
TestResult cr = h.ChiSquaredTest(d);
nullDistribution = (ChiSquaredDistribution) cr.Distribution;
//Console.WriteLine(((ChiSquaredDistribution) cr.Distribution).DegreesOfFreedom);
s.Add(cr.Statistic);
}
Console.WriteLine(nullDistribution.DegreesOfFreedom);
TestResult kr = s.KuiperTest(nullDistribution);
Console.WriteLine(kr.LeftProbability);
}
