52 lines
1.1 KiB
C
52 lines
1.1 KiB
C
|
|
/* Driver for routine kendl1 */
|
|
|
|
#include <stdio.h>
|
|
#define NRANSI
|
|
#include "nr.h"
|
|
#include "nrutil.h"
|
|
|
|
#define NDAT 200
|
|
|
|
int main(void)
|
|
{
|
|
int i,j;
|
|
long idum;
|
|
float prob,tau,z,*data1,*data2;
|
|
static char *txt[5]={"RAN0","RAN1","RAN2","RAN3","RAN4"};
|
|
|
|
data1=vector(1,NDAT);
|
|
data2=vector(1,NDAT);
|
|
/* Look for correlations in RAN0, RAN1, RAN2, RAN3 and RAN4 */
|
|
printf("\nPair correlations of RAN0 ... RAN4\n\n");
|
|
printf("%9s %17s %16s %18s\n",
|
|
"Program","Kendall tau","Std. Dev.","Probability");
|
|
for (i=1;i<=5;i++) {
|
|
idum=(-1357);
|
|
for (j=1;j<=NDAT;j++) {
|
|
if (i == 1) {
|
|
data1[j]=ran0(&idum);
|
|
data2[j]=ran0(&idum);
|
|
} else if (i == 2) {
|
|
data1[j]=ran1(&idum);
|
|
data2[j]=ran1(&idum);
|
|
} else if (i == 3) {
|
|
data1[j]=ran2(&idum);
|
|
data2[j]=ran2(&idum);
|
|
} else if (i == 4) {
|
|
data1[j]=ran3(&idum);
|
|
data2[j]=ran3(&idum);
|
|
} else if (i == 5) {
|
|
data1[j]=ran4(&idum);
|
|
data2[j]=ran4(&idum);
|
|
}
|
|
}
|
|
kendl1(data1,data2,NDAT,&tau,&z,&prob);
|
|
printf("%8s %17.6f %17.6f %17.6f\n",txt[i-1],tau,z,prob);
|
|
}
|
|
free_vector(data2,1,NDAT);
|
|
free_vector(data1,1,NDAT);
|
|
return 0;
|
|
}
|
|
#undef NRANSI
|