// g02bg Example Program Text // C# version, NAG Copyright 2008 using System; using NagLibrary; using System.Globalization; using System.IO; namespace NagDotNetExamples { public class G02BGE { static string datafile = "ExampleData/g02bge.d"; static void Main(String[] args) { if (args.Length == 1) { datafile = args[0]; } StartExample(); } public static void StartExample() { try { DataReader sr = new DataReader(datafile); int i, j; const int m = 4; const int n = 5; const int nvars = 3; double[,] r = new double[nvars, nvars]; double[,] ssp = new double[nvars, nvars]; double[] std = new double[nvars]; double[,] x = new double[n, m]; double[] xbar = new double[nvars]; int[] kvar = new int[nvars]; int ifail; Console.WriteLine("g02bg Example Program Results"); // Skip heading in data file sr.Reset(); kvar[0] = 4; kvar[1] = 1; kvar[2] = 2; sr.Reset(); for (i = 1; i <= n; i++) { for (j = 1; j <= m; j++) { x[i - 1, j - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture); } } Console.WriteLine(" "); Console.WriteLine(" {0}{1,5}", "Number of variables (columns) =", m); Console.WriteLine(" {0}{1,5}", "Number of cases (rows) =", n); Console.WriteLine(" "); Console.Write(" {0}", "Data matrix is:-"); Console.Write("\r\n "); for (j = 1; j <= m; j++) { Console.Write("{0,10}", j); } Console.WriteLine(" "); for (i = 1; i <= n; i++) { Console.Write(" " + i + " "); for (j = 1; j <= m; j++) { Console.Write(" {0,8:f4}", x[i - 1, j - 1]); } Console.WriteLine(" "); } Console.WriteLine(" "); Console.WriteLine(" "); // G02.g02bg(n, m, x, nvars, kvar, xbar, std, ssp, r, out ifail); // if (ifail < 0) { Console.WriteLine(" "); Console.WriteLine("** g02bg failed with ifail = {0,5}", ifail); } else { if (ifail != 0) { Console.Write(" {0}{1,5}", "Routine fails, ifail =", ifail); } else { Console.Write(" {0}", "Variable Mean St. dev.\n"); for (i = 1; i <= nvars; i++) { Console.Write(" {0} {1,10:f4}{2,10:f4}", kvar[i - 1], xbar[i - 1], std[i - 1]); Console.WriteLine(" "); } Console.WriteLine(" "); Console.WriteLine(" "); Console.Write(" {0}", "Sums of squares and cross-products of deviations\n"); for (i = 1; i <= nvars; i++) { Console.Write(" {0,14}", kvar[i - 1]); } Console.WriteLine(" "); for (i = 1; i <= nvars; i++) { Console.Write(" {0}", kvar[i - 1]); for (j = 1; j <= nvars; j++) { Console.Write(" \t{0,10:f4}", ssp[i - 1, j - 1]); } Console.WriteLine(" "); } Console.WriteLine(" "); Console.Write(" {0}", "Correlation coefficients\n"); Console.Write(" "); for (i = 1; i <= nvars; i++) { Console.Write(" {0,11}", kvar[i - 1]); } Console.WriteLine(" "); for (i = 1; i <= nvars; i++) { Console.Write(" {0}", kvar[i - 1]); for (j = 1; j <= nvars; j++) { Console.Write(" \t{0,10:f4}", r[i - 1, j - 1]); } Console.WriteLine(" "); } Console.WriteLine(" "); } } // } catch (Exception e) { Console.WriteLine(e.Message); Console.Write( "Exception Raised"); } } } }