BEGIN {
	 OFS=FS =","
	 All="__all";
   }
  
     {sub(/#.*/,""); gsub(/[ \t]*/,"") } # no comments or white space
/^$/ {next}

/@attribute/ { Nump[++Attr]= $0 ~ /numeric|real|integer/ }
/@/          { next; }
             { store() }
END          { 
			   for(I in Nump) 
				   	print (Nump[Attr] ? relativeSd(I) : relativeEnd(I))
	         }

function keep(    i,c) { 
   Total++;
   c=$NF;
   Classes[c]++;
   for(i=1;i<NF;i++) {
     if ($i=="?") continue;
     if (Nump[i]) {
		 	Sum[  i,$i] += c;
			SumSq[i,$i] += c*c;
			N[    i,$i]++
     } else {
		 Freq[c,i,$i]++;  Freq[All,i,$i]++
         if (++Seen[i,$i]==1) Attributes[i,++Attributes[i,0]] = $i
     }
}
function relativeSd(i) {
	for(c in Classes)
	     tmp[c] += Classes[c]/Total * sd(c,i);
	return tmp[c];
}
function sd(c,i)           { return sd0(SumSq[c,i],SumX[c,i],N[c,i]) }
function sd0(sumSq,sumX,n) { return sqrt((sumSq-((sumX*sumX)/n))/(n-1)) }

function relativeEnt(i,   c,tmp) {
	for(c in Classes)
		tmp[c] = Classes[c]/Total * ent(c,i)
	return tmp[c];
}

function log2(x) { return log(x)/log(2) }

function ent(c,i,  max,j,what,n,p,e) { 
		max=Attributes[i,0];
		for(j=i;j<=max;j++) {
			what = Attributes[i,j];
			n[j] = Freq[c,i,what];
			ns++;
	    }
 		for(j=1;j<=max;j++)
			p = n[j]/ns
			e += -1*p * log2(p) }
		}
		return e
	}

function relativeEnd(i) {
}