|
|
An application of using the original
distribution of NAG library.
In a Ch shell or Ch IDE, run the following command:
ch g02brce.c < g02brce.d
You will get the output:
g02brc Example Program Results
Correlation coefficients:
1.00000 0.85714 0.12849
0.71429 1.00000 0.33040
0.10287 0.41148 1.00000
The contents for files g02brce.c and g02brce.d are listed below.
<![CDATA[
/* nag_ken_spe_corr_coeff(g02brc) Example Program
*
* Copyright 1992 Numerical Algorithms Group.
*
* Mark 3, 1993.
*/
#include <nag.h>
#include <stdio.h>
#include <nag_stdlib.h>
#include <nagg02.h>
#define NMAX 20
#define MMAX 20
int main(void)
{
Integer i, j, m, n, tdx, tdc;
double x[NMAX][MMAX], corr[MMAX][MMAX];
Integer svar[MMAX], sobs[NMAX];
Integer *svarptr, *sobsptr;
char s, w;
tdx = MMAX;
tdc = MMAX;
Vprintf("g02brc Example Program Results\n");
/* Skip heading in data file */
Vscanf(" %*[^\n]");
/* Read data */
Vscanf("%ld%ld\n",&m,&n);
if (m>=1 && m<=MMAX && n>=1 && n<=NMAX)
{
for (i=0; i<n; i++)
for (j=0; j<m; j++)
Vscanf("%lf",&x[i][j]);
/* Read flag specifying if svar is to be supplied */
Vscanf(" %c",&s);
if (s == 'S' || s == 's')
{
/* Assign pointer to svar and read in values for svar */
svarptr = svar;
for (i=0; i<m; i++)
Vscanf("%ld",&svar[i]);
}
else
{
/* Assign pointer to NULL and discard rest of line */
svarptr = (Integer *)0;
Vscanf(" %*[^\n]"); /* skip rest of line */
}
/* Read flag specifying if sobs is to be supllied */
Vscanf(" %c",&w);
if (w == 'W' || w == 'w')
{
/* Assign pointer to sobs and read in values for sobs */
sobsptr = sobs;
for (i=0; i<n; i++)
Vscanf("%ld",&sobs[i]);
}
else
{
/* Assign pointer to NULL and discard rest of line */
sobsptr = (Integer *)0;
Vscanf(" %*[^\n]"); /* skip rest of line */
}
/* Calculate the Kendall and Spearman coefficients */
g02brc(n, m, &x[0][0], tdx, svarptr, sobsptr, &corr[0][0], tdc,
NAGERR_DEFAULT);
Vprintf("\nCorrelation coefficients:\n\n");
for (i=0; i<m; i++)
{
for (j=0; j<m; j++)
Vprintf("%8.5f ", corr[i][j]);
Vprintf("\n");
}
}
else
{
Vfprintf(stderr, "One or both of m and n are out of range:"
" m = %-3ld while n = %-3ld\n", m, n);
return EXIT_FAILURE;
}
return EXIT_SUCCESS;
}
]]>
The content of
file "g02brc.d" is as follows.
<![CDATA[
g02brc Example Program Data
3 7
1.0 2.0 4.0
7.0 7.0 3.0
2.0 3.0 4.0
4.0 4.0 5.0
5.0 6.0 7.0
3.0 1.0 3.0
6.0 5.0 5.0
s 1 1 1
w 1 1 1 1 1 1 1
]]>
|