#include <NAG\OCN_g08.h>
void ks_test()
{
	// suppose the active window is a worksheet
	Worksheet wks = Project.ActiveLayer();
	if(!wks)
		return;
	
	// suppose the test data is stored in column(A)
	Dataset dsA(wks,0); 
	vector vData =dsA;
	
	int nRet;
	
	double par[2];
	Nag_ParaEstimates estima;
	estima = Nag_ParaEstimated;
	
	Nag_TestStatistics dtype;
	dtype = Nag_TestStatisticsDAbs;
	double dD, dZ, dP;
	
	// perform a one sample Kolmogorov-Smirnov test to examine if the data obeys Normal distribution
	nRet = nag_1_sample_ks_test(vData.GetSize(), vData, Nag_Normal, par, estima, dtype, &dD, &dZ, &dP);
	
	// output test results, where "Tail probability" tells the probability that the data is from a Normal distribution
	printf("Test against Normal distribution with parameters estimated from the data\n");
	printf("\n");
	printf("Mean			= %8.4f\n", par[0]);
	printf("Variance		= %8.4f\n", par[1]);
	printf("Test statistics D	= %8.4f\n", dD);
	printf("Z statistic		= %8.4f\n", dZ);
	printf("Tail probability	= %8.4f\n", dP);
}