process.awk

来自「Scalable k-means software and test datas」· AWK 代码 · 共 27 行

{	if(!(i%5))
		d = 1e10;
	if($1 < d)
		d = $1;
	i++;
	if(!(i%5))
	{	# print d;
		d1 += d;
		d2 += d*d;
		n++;
	}
	t += $2;
	t2 += $2*$2;
	tn++;
}

END {
	m = d1/n;
	d = sqrt((d2/n - m*m)*n/(n - 1));

	mt = t/tn;
	md = sqrt((t2/tn - mt*mt)*tn/(tn - 1));

	print m, d, d/sqrt(n), mt, md, md/sqrt(tn);
}