package?day02;

public?class?SPT?{
	/**
	?*?
	?*?多元線性回歸分析
	?*
	?*?
	?*?@param?x[m][n]?每一列存放?m?個自變量的觀察值
	?*?
	?*?@param?y[n]?存放隨即變量?y?的?n?個觀察值
	?*?
	?*?@param?m?自變量的個數
	?*?
	?*?@param?n?觀察數據的組數
	?*?
	?*?@param?a?返回回歸系數?a0...am
	?*?
	?*?@param?dt[4]?dt[0]?偏差平方和?qdt[1]?平均標準偏差?s?dt[2]?返回復相關系數?r?dt[3]?返回回歸平方和?u
	?*?
	?*?@param?v[m]?返回?m?個自變量的偏相關系數
	?*/
	public?static?void?sqt2（double[][]?x?double[]?y?int?m?int?n?double[]?a?double[]?dt?double[]?v）?{
		int?i?j?k?mm;
		double?q?e?u?p?yy?s?r?pp;
		double[]?b?=?new?double[（m?+?1）?*?（m?+?1）];
		mm?=?m?+?1;
		b[mm?*?mm?-?1]?=?n;
		for?（j?=?0;?j?<=?m?-?1;?j++）?{
			p?=?0.0;
			for?（i?=?0;?i?<=?n?-?1;?i++）
				p?=?p?+?x[j][i];
			b[m?*?mm?+?j]?=?p;
			b[j?*?mm?+?m]?=?p;
		}
		for?（i?=?0;?i?<=?m?-?1;?i++）
			for?（j?=?i;?j?<=?m?-?1;?j++）?{
				p?=?0.0;
				for?（k?=?0;?k?<=?n?-?1;?k++）
					p?=?p?+?x[i][k]?*?x[j][k];
				b[j?*?mm?+?i]?=?p;
				b[i?*?mm?+?j]?=?p;
			}
		a[m]?=?0.0;
		for?（i?=?0;?i?<=?n?-?1;?i++）
			a[m]?=?a[m]?+?y[i];
		for?（i?=?0;?i?<=?m?-?1;?i++）?{
			a[i]?=?0.0;
			for?（j?=?0;?j?<=?n?-?1;?j++）
				a[i]?=?a[i]?+?x[i][j]?*?y[j];
		}
		chlk（b?mm?1?a）;
		yy?=?0.0;
		for?（i?=?0;?i?<=?n?-?1;?i++）
			yy?=?yy?+?y[i]?/?n;
		q?=?0.0;
		e?=?0.0;
		u?=?0.0;
		for?（i?=?0;?i?<=?n?-?1;?i++）?{
			p?=?a[m];
			for?（j?=?0;?j?<=?m?-?1;?j++）
				p?=?p?+?a[j]?*?x[j][i];
			q?=?q?+?（y[i]?-?p）?*?（y[i]?-?p）;
			e?=?e?+?（y[i]?-?yy）?*?（y[i]?-?yy）;
			u?=?u?+?（yy?-?p）?*?（yy?-?p）;
		}
		s?=?Math.sqrt（q?/?n）;
		r?=?Math.sqrt（1.0?-?q?/?e）;
		for?（j?=?0;?j?<=?m?-?1;?j++）?{
			p?=?0.0;
			for?（i?=?0;?i?<=?n?-?1;?i++）?{
				pp?=?a[m];
				for?（k?=?0;?k?<=?m?-?1;?k++）
					if?（k?!=?j）
						pp?=?pp?+?a[k]?*?x[k][i];
				p?=?p?+?（y[i]?-?pp）?*?（y[i]?-?pp）;
			}
			v[j]?=?Math.sqrt（1.0?-?q?/?p）;
		}
		dt[0]?=?q;
		dt[1]?=?s;
		dt[2]?=?r;
		dt[3]?=?u;
	}

	private?static?int?chlk（double[]?a?int?n?int?m?double[]?d）?{
		int?i?j?