cluster

聚类算法全集以及内附数据集

#!/bin/bash

cnt=${1:-20}

#-----------------------------------------------------------------------

function cluster ()
{
  rm -f x y z
  if [[ $3 = "quick" ]]; then x="-g2"; else x=""; fi
  args="-e10000 -H-4 -T1e-6 $x $1.pat x";
  for (( i = 0; i < $cnt; i++ )); do
    mcli -c$2 -A$3 $4 $5 -S$(( (7*i)+1 )) $args 2> y
    e=`sed -n "s/.*\[\([0-9]*\) epoch(s)\].*/\1/p" y`
    echo $e >> z
    if (( e > 1000 )); then break; fi;
  done
  n=`gawk 'BEGIN { e = s = n = 0; }
           { e += $1; s += $1*$1; n++; }
           END { if (n <= 0) n = 1;
                 s = sqrt((s -e*e/n)/n); e /= n; 
                 printf(":  %.1f (~%.1f)\n", e, s); }' z`
  echo $1 $2 $3 $4 $5 $n
  echo $1 $2 $3 $4 $5 $n >> cluster.out
  rm -f x y z
}

#-----------------------------------------------------------------------

function single ()
{
  cluster $1 $2 none      $3
  cluster $1 $2 expand    $3 -g1.05
  cluster $1 $2 expand    $3 -g1.1
  cluster $1 $2 expand    $3 -g1.2
  cluster $1 $2 expand    $3 -g1.3
  cluster $1 $2 expand    $3 -g1.4
  cluster $1 $2 expand    $3 -g1.5
  cluster $1 $2 expand    $3 -g1.6
  cluster $1 $2 expand    $3 -g1.7
  cluster $1 $2 expand    $3 -g1.8
  cluster $1 $2 expand    $3 -g1.9
  cluster $1 $2 expand    $3 -g2.0
  cluster $1 $2 momentum  $3 -m0.05
  cluster $1 $2 momentum  $3 -m0.1
  cluster $1 $2 momentum  $3 -m0.15
  cluster $1 $2 momentum  $3 -m0.2
  cluster $1 $2 momentum  $3 -m0.25
  cluster $1 $2 momentum  $3 -m0.3
  cluster $1 $2 momentum  $3 -m0.35
  cluster $1 $2 momentum  $3 -m0.4
  cluster $1 $2 momentum  $3 -m0.45
  cluster $1 $2 momentum  $3 -m0.5
  cluster $1 $2 momentum  $3 -m0.55
  cluster $1 $2 momentum  $3 -m0.6
  cluster $1 $2 momentum  $3 -m0.65
  cluster $1 $2 momentum  $3 -m0.7
  cluster $1 $2 momentum  $3 -m0.75
  cluster $1 $2 momentum  $3 -m0.8
  cluster $1 $2 momentum  $3 -m0.85
  cluster $1 $2 momentum  $3 -m0.9
  cluster $1 $2 adaptive  $3
  cluster $1 $2 resilient $3
  cluster $1 $2 quick     $3
}

#-----------------------------------------------------------------------

function all ()
{
  echo "----- $1 -----" >> cluster.out
  single abalone 3 $1
  single abalone 6 $1
  single breast  2 $1
  single iris    3 $1
  single wine    3 $1
  single wine    6 $1
}

#-----------------------------------------------------------------------

function abalone3 ()
{
  cluster abalone 3 none
  cluster abalone 3 expand       -g1.5
  cluster abalone 3 momentum     -m0.3
  cluster abalone 3 adaptive
  cluster abalone 3 resilient
  cluster abalone 3 quick

  cluster abalone 3 none      -v
  cluster abalone 3 expand    -v -g1.6
  cluster abalone 3 momentum  -v -m0.45
  cluster abalone 3 adaptive  -v
  cluster abalone 3 resilient -v
  cluster abalone 3 quick     -v

  cluster abalone 3 none      -V
  cluster abalone 3 expand    -V -g1.8
  cluster abalone 3 momentum  -V -m0.6
  cluster abalone 3 adaptive  -V
  cluster abalone 3 resilient -V
  cluster abalone 3 quick     -V
}

#-----------------------------------------------------------------------

function abalone6 ()
{
  cluster abalone 6 none
  cluster abalone 6 expand       -g1.8
  cluster abalone 6 momentum     -m0.5
  cluster abalone 6 adaptive
  cluster abalone 6 resilient
  cluster abalone 6 quick

  cluster abalone 6 none      -v
  cluster abalone 6 expand    -v -g1.8
  cluster abalone 6 momentum  -v -m0.7
  cluster abalone 6 adaptive  -v
  cluster abalone 6 resilient -v
  cluster abalone 6 quick     -v

  cluster abalone 6 none      -V
  cluster abalone 6 expand    -V -g1.9
  cluster abalone 6 momentum  -V -m0.65
  cluster abalone 6 adaptive  -V
  cluster abalone 6 resilient -V
  cluster abalone 6 quick     -V
}

#-----------------------------------------------------------------------

function breast2 ()
{
  cluster breast 2 none
  cluster breast 2 expand       -g1.2
  cluster breast 2 momentum     -m0.05
  cluster breast 2 adaptive
  cluster breast 2 resilient
  cluster breast 2 quick

  cluster breast 2 none      -v
  cluster breast 2 expand    -v -g1.6
  cluster breast 2 momentum  -v -m0.25
  cluster breast 2 adaptive  -v
  cluster breast 2 resilient -v
  cluster breast 2 quick     -v

  cluster breast 2 none      -V
  cluster breast 2 expand    -V -g1.5
  cluster breast 2 momentum  -V -m0.25
  cluster breast 2 adaptive  -V
  cluster breast 2 resilient -V
  cluster breast 2 quick     -V
}

#-----------------------------------------------------------------------

function iris3 ()
{
  cluster iris 3 none
  cluster iris 3 expand       -g1.4
  cluster iris 3 momentum     -m0.15
  cluster iris 3 adaptive
  cluster iris 3 resilient
  cluster iris 3 quick

  cluster iris 3 none      -v
  cluster iris 3 expand    -v -g1.6
  cluster iris 3 momentum  -v -m0.3
  cluster iris 3 adaptive  -v
  cluster iris 3 resilient -v
  cluster iris 3 quick     -v

  cluster iris 3 none      -V
  cluster iris 3 expand    -V -g1.6
  cluster iris 3 momentum  -V -m0.25
  cluster iris 3 adaptive  -V
  cluster iris 3 resilient -V
  cluster iris 3 quick     -V
}

#-----------------------------------------------------------------------

function wine3 ()
{
  cluster wine 3 none
  cluster wine 3 expand       -g1.7
  cluster wine 3 momentum     -m0.3
  cluster wine 3 adaptive
  cluster wine 3 resilient
  cluster wine 3 quick

  cluster wine 3 none      -v
  cluster wine 3 expand    -v -g1.7
  cluster wine 3 momentum  -v -m0.3
  cluster wine 3 adaptive  -v
  cluster wine 3 resilient -v
  cluster wine 3 quick     -v

  cluster wine 3 none      -V
  cluster wine 3 expand    -V -g1.9
  cluster wine 3 momentum  -V -m0.7
  cluster wine 3 adaptive  -V
  cluster wine 3 resilient -V
  cluster wine 3 quick     -V
}

#-----------------------------------------------------------------------

function wine6 ()
{
  cluster wine 6 none
  cluster wine 6 expand       -g1.9
  cluster wine 6 momentum     -m0.75
  cluster wine 6 adaptive
  cluster wine 6 resilient
  cluster wine 6 quick

  cluster wine 6 none      -v
  cluster wine 6 expand    -v -g2.0
  cluster wine 6 momentum  -v -m0.5
  cluster wine 6 adaptive  -v
  cluster wine 6 resilient -v
  cluster wine 6 quick     -v

  cluster wine 6 none      -V
  cluster wine 6 expand    -V -g1.9
  cluster wine 6 momentum  -V -m0.7
  cluster wine 6 adaptive  -V
  cluster wine 6 resilient -V
  cluster wine 6 quick     -V
}

#-----------------------------------------------------------------------

rm -f cluster.out

all           # fuzzy c-means algorithm
all -v        # axes-parallel Gustafson-Kessel algorithm
all -V        # normal        Gustafson-Kessel algorithm

#abalone3
#abalone6
#breast2
#iris3
#wine3
#wine6