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 + ""); }