NR<R1  { next }
NR==R1 {for(I=C1;I<=C2;I++) {
	    Min[I] = 10^30
	    Max[I] = -1*Min[I]
            
    }}
{for(I=C1;I<=C2;I++) {
	Mean[I]+= $I;
	N[I]++;
        Min[I]= $I < Min[I] ? $I : Min[I];
	Max[I]= $I > Max[I] ? $I : Max[I];
	if ( 1 ==  ++Unique[I,$I])
	    Count[I]++
	
    }
}
END { 
    for(I in Mean) Mean[I]= Mean[I]/N[I];
    dump("min",Min);
    dump("mean",Mean);
    dump("max",Max);
    dump("count",Count)
}
function dump(label,a,    sep,i) {
    sep = sep ? sep : OFS;
    printf "!" label
    for(i=C1;i<=C2;i++)
	printf("%s%s%s%s", OFS,i,OFS,a[i])
    print ""
}