DataMuseum.dk

top - download

    Length: 547 (0x223)
    Types: TextFile
    Names: »anovahist«

└─⟦87ddcff64⟧ Bits:30001253 CPHDIST85 Tape, 1985 Autumn Conference Copenhagen
    └─ ⟦this⟧ »cph85dist/stat/scripts/anovahist« 

TextFile

#!/bin/sh
#compare levels of a variable using histograms
#the data is assumed to be in the last column (x[N])
#Usage: anovahist column datafile intwidth min

var=$1
dfile=$2
intwidth=${3-1}
min=${4-0}

l1=`dm s$var < $dfile | sort | uniq`
for v1 in $l1
do
		(
		echo "$v1"
		dm "s$var='$v1' ? x[N] : KILL" < $dfile |
			desc -i $intwidth -m $min -h |
			grep -v Freq |
			colrm 1 12
		) > /tmp/$$$v1
done
series -1 53 | dm "INLINE=1?'Level:':$min+x1*$intwidth" > /tmp/$$
pr -h "variable: $var in datafile: $dfile" -m - /tmp/$$*
/bin/rm -f /tmp/$$