//LinearRegressionND.cpp	n元線性回歸分析

#include <iostream>			//模板類輸入輸出流iostream標準頭文件
#include <valarray>			//模板類valarray的標準頭文件
#include <Matrix.h>			//模板類Matrix頭文件
#include <Statistic.h>		//數據處理與回歸分析頭文件
using namespace std;		//名字空間

void main(void)
{	
	valarray<long double> a(4);
	valarray<long double> v(3);
	valarray<long double> dt(4);
	long double x0(192), h(2);

	long double x[3][5] = 
	{
		{1.1,1.0,1.2,1.1,0.9},
        {2.0,2.0,1.8,1.9,2.1},
		{3.2,3.2,3.0,2.9,2.9}
	};
	
	long double y[5] = {10.1,10.2,10.0,10.1,10.0};
	 
	matrix<long double> xx(&x[0][0], 3, 5);
	valarray<long double> yy(y, 5);

    LinearRegressionND(xx,yy,a,dt,v);

	cout << endl;

	for (int i=0; i<4; i++)
		cout << " a(" << i << ") = " << a[i] << endl;
	
	cout << endl << " q = " << dt[0] << "\t s = " << dt[1];
								cout <<  "\t r = " << dt[2] << endl;
	
	cout << endl;
	
	for (i=0; i<3; i++)
		cout << " v(" << i << ") = " << v[i] << endl;

	cout << endl;

	cout << " u = " << dt[3] << endl << endl;
}