// g02bs Example Program Text // C# version, NAG Copyright 2008 using System; using NagLibrary; using System.Globalization; using System.IO; namespace NagDotNetExamples { public class G02BSE { static string datafile = "ExampleData/g02bse.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, itype, j, ncases; const int m = 3; const int n = 9; double[,] cnt = new double[m, m]; double[,] rr = new double[m, m]; double[,] x = new double[n, m]; double[] xmiss = new double[m]; int[] miss = new int[m]; int ifail; Console.WriteLine("g02bs Example Program Results"); // Skip heading in data file sr.Reset(); 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,7} ", j); } Console.WriteLine(" "); for (i = 1; i <= n; i++) { Console.Write(" {0}", i); for (j = 1; j <= m; j++) { Console.Write(" {0,8:f4}", x[i - 1, j - 1]); } Console.WriteLine(" "); } Console.WriteLine(" "); Console.WriteLine(" "); // // Set up missing values before calling method // miss[0] = 1; miss[1] = 1; miss[2] = 1; xmiss[0] = 0.990e0; xmiss[1] = 9.00e0; xmiss[2] = 0.000e0; itype = 0; // G02.g02bs(n, m, x, miss, xmiss, itype, rr, out ncases, cnt, out ifail); // if (ifail < 0) { Console.WriteLine(" "); Console.WriteLine("** g02bs failed with ifail = {0,5}", ifail); } else { if (ifail != 0) { Console.Write(" {0}{1,5}", "Routine fails, ifail =", ifail); } else { Console.WriteLine(" {0}", "Matrix of rank correlation coefficients:"); Console.WriteLine(" {0}", "Upper triangle -- Spearman's"); Console.WriteLine(" {0}", "Lower triangle -- Kendall's tau"); Console.WriteLine(" "); for (i = 1; i <= m; i++) { Console.Write(" {0,10}", i); } Console.WriteLine(" "); for (i = 1; i <= m; i++) { Console.Write(" {0}", i); for (j = 1; j <= m; j++) { Console.Write(" {0,10:f4}", rr[i - 1, j - 1]); } Console.WriteLine(" "); } Console.WriteLine(" "); Console.WriteLine(" "); Console.Write(" {0}{1,5}", "Minimum number of cases used for any pair of variables:", ncases); Console.WriteLine(" "); Console.Write(" {0}", "Numbers used for each pair are:\n\n"); for (i = 1; i <= m; i++) { Console.Write(" {0,10}", i); } Console.WriteLine(" "); for (i = 1; i <= m; i++) { Console.Write(" {0}", i); for (j = 1; j <= m; j++) { Console.Write(" {0,10:f4}", cnt[i - 1, j - 1]); } Console.WriteLine(" "); } Console.WriteLine(" "); } } // } catch (Exception e) { Console.WriteLine(e.Message); Console.Write( "Exception Raised"); } } } }