private static void runs_sample_test()
//****************************************************************************80
//
// Purpose:
//
// RUNS_SAMPLE_TEST tests RUNS_SAMPLE.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 03 April 2016
//
// Author:
//
// John Burkardt
//
{
const int SAMPLE_NUM = 1000;
int i;
int seed = 123456789;
int[] x = new int[SAMPLE_NUM];
Console.WriteLine("");
Console.WriteLine("RUNS_SAMPLE_TEST");
Console.WriteLine(" RUNS_MEAN computes the Runs mean;");
Console.WriteLine(" RUNS_SAMPLE samples the Runs distribution;");
Console.WriteLine(" RUNS_VARIANCE computes the Runs variance");
const int m = 10;
const int n = 5;
Console.WriteLine("");
Console.WriteLine(" PDF parameter M = " + m + "");
Console.WriteLine(" PDF parameter N = " + n + "");
double mean = Runs.runs_mean(m, n);
double variance = Runs.runs_variance(m, n);
Console.WriteLine(" PDF mean = " + mean + "");
Console.WriteLine(" PDF variance = " + variance + "");
for (i = 0; i < SAMPLE_NUM; i++)
{
x[i] = Runs.runs_sample(m, n, ref seed);
}
mean = typeMethods.i4vec_mean(SAMPLE_NUM, x);
variance = typeMethods.i4vec_variance(SAMPLE_NUM, x);
int xmax = typeMethods.i4vec_max(SAMPLE_NUM, x);
int xmin = typeMethods.i4vec_min(SAMPLE_NUM, x);
Console.WriteLine("");
Console.WriteLine(" Sample size = " + SAMPLE_NUM + "");
Console.WriteLine(" Sample mean = " + mean + "");
Console.WriteLine(" Sample variance = " + variance + "");
Console.WriteLine(" Sample maximum = " + xmax + "");
Console.WriteLine(" Sample minimum = " + xmin + "");
}