#!/usr/bin/env -S davinci -f  

###########################  was:  #!/usr/local/bin/davinci -f


#
# cluster analysis
# usage: cluster filename
#
verbose = 0
scale=3
echo( "Cluster analysis: " + basename($1))
data = load($1)
cdata = cluster(data)
d = ccount(cdata, threshold=3)
echo( "             3/9: " + (d[1]*100) +  "%  " + d[2] + "/" + d[3])
d = ccount(cdata, threshold=4)
echo( "             4/9: " + (d[1]*100) +  "%  " + d[2] + "/" + d[3])
d = ccount(cdata, threshold=5)
echo( "             5/9: " + (d[1]*100) +  "%  " + d[2] + "/" + d[3])
d = ccount(cdata, threshold=6)
echo( "             6/9: " + (d[1]*100) +  "%  " + d[2] + "/" + d[3])
d = ccount(cdata, threshold=7)
echo( "             7/9: " + (d[1]*100) +  "%  " + d[2] + "/" + d[3])
d = ccount(cdata, threshold=8)
echo( "             8/9: " + (d[1]*100) +  "%  " + d[2] + "/" + d[3])
d = ccount(cdata, threshold=9)
echo( "             9/9: " + (d[1]*100) +  "%  " + d[2] + "/" + d[3])
