diag

it is the Data Mining Algorithm source code.

#!/bin/bash

#-----------------------------------------------------------------------
# default parameters
#-----------------------------------------------------------------------

file=""

width=52
height=38

trial="0.75white:"
mean="black:"
none="black:"
limit="red:"
red="red:"
green="0.8green:"
blue="blue:"

#-----------------------------------------------------------------------
# option evaluation
#-----------------------------------------------------------------------

while getopts "o:w:h:c" option; do
  case "$option" in
  o) file="$OPTARG";;
  w) width="$OPTARG";;
  h) height="$OPTARG";;
  c) trial="0.75white:"
     mean="black:"
     none="0.6white:"
     limit="0.6white:dashed withdots scaled 0.5"
     red="black:dashed evenly"
     green="black:dashed withdots scaled 0.5"
     blue="black:";;
  esac
done

#-----------------------------------------------------------------------
# functions
#-----------------------------------------------------------------------

function start ()
{
  out="$1"                      # note the output file name
  echo    "input dataplot.mp"                     > "$out"
  echo    ""                                     >> "$out"
  echo    "wd := "$width"mm;"                    >> "$out";
  echo    "ht := "$height"mm;"                   >> "$out"
  echo    ""                                     >> "$out"
  echo -n "%-----------------------------------" >> "$out"
  echo    "------------------------------------" >> "$out"
  echo    ""                                     >> "$out"
  i=0                           # initialize the figure identifier
}  # start ()

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

function term ()
{
  echo "end." >> "$out"
}  # term ()

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

list=${@:$OPTIND}
if [[ $list == ""     ]]; then list="all";     fi
if [[ $file == "diag" ]]; then file="diag.mp"; fi
if [[ $file != ""     ]]; then start "$file";  fi

for data in $list; do

#-----------------------------------------------------------------------
# iris
#-----------------------------------------------------------------------

if [[ $data == "all" || $data == "iris" ]]; then
if [[ $file == "" ]]; then start iris.mp; fi
echo -n "processing iris ... "

mpdiag $i \
  "log(difference):botlft | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   none00:diff:$trial:trials:0 | \
   none01:diff:$trial::-1 | \
   none02:diff:$trial::-1 | \
   none03:diff:$trial::-1 | \
   none04:diff:$trial::-1 | \
   none05:diff:$trial::-1 | \
   none06:diff:$trial::-1 | \
   none07:diff:$trial::-1 | \
   none08:diff:$trial::-1 | \
   none09:diff:$trial::-1 | \
   none10:diff:$trial::-1 | \
   none11:diff:$trial::-1 | \
   none12:diff:$trial::-1 | \
   none13:diff:$trial::-1 | \
   none14:diff:$trial::-1 | \
   none15:diff:$trial::-1 | \
   none16:diff:$trial::-1 | \
   none17:diff:$trial::-1 | \
   none18:diff:$trial::-1 | \
   none19:diff:$trial::-1 | \
   none:diff:$mean:average:1" \
   iris.trl "$out"
let i++

mpdiag $i \
  "log(difference):botlft | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vnone00:diff:$trial:trials:0 | \
   vnone01:diff:$trial::-1 | \
   vnone02:diff:$trial::-1 | \
   vnone03:diff:$trial::-1 | \
   vnone04:diff:$trial::-1 | \
   vnone05:diff:$trial::-1 | \
   vnone06:diff:$trial::-1 | \
   vnone07:diff:$trial::-1 | \
   vnone08:diff:$trial::-1 | \
   vnone09:diff:$trial::-1 | \
   vnone10:diff:$trial::-1 | \
   vnone11:diff:$trial::-1 | \
   vnone12:diff:$trial::-1 | \
   vnone13:diff:$trial::-1 | \
   vnone14:diff:$trial::-1 | \
   vnone15:diff:$trial::-1 | \
   vnone16:diff:$trial::-1 | \
   vnone17:diff:$trial::-1 | \
   vnone18:diff:$trial::-1 | \
   vnone19:diff:$trial::-1 | \
   vnone:diff:$mean:average:1" \
   iris.trl "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   Vnone00:diff:$trial:trials:0 | \
   Vnone01:diff:$trial::-1 | \
   Vnone02:diff:$trial::-1 | \
   Vnone03:diff:$trial::-1 | \
   Vnone04:diff:$trial::-1 | \
   Vnone05:diff:$trial::-1 | \
   Vnone06:diff:$trial::-1 | \
   Vnone07:diff:$trial::-1 | \
   Vnone08:diff:$trial::-1 | \
   Vnone09:diff:$trial::-1 | \
   Vnone10:diff:$trial::-1 | \
   Vnone11:diff:$trial::-1 | \
   Vnone12:diff:$trial::-1 | \
   Vnone13:diff:$trial::-1 | \
   Vnone14:diff:$trial::-1 | \
   Vnone15:diff:$trial::-1 | \
   Vnone16:diff:$trial::-1 | \
   Vnone17:diff:$trial::-1 | \
   Vnone18:diff:$trial::-1 | \
   Vnone19:diff:$trial::-1 | \
   Vnone:diff:$mean:average:1" \
   iris.trl "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   none:diff:$none:none:0 | \
   e1.1:diff:$blue:1.1/1.3/1.5:1 | \
   e1.3:diff:$blue::-1 | \
   e1.5:diff:$blue::-1 | \
   e1.6:diff:$limit:1.6:2" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   none:diff:$none:none:0 | \
   m0.1:diff:$blue:0.1/0.15/0.4:1 | \
   m0.15:diff:$blue::-1 | \
   m0.4:diff:$blue::-1 | \
   m0.5:diff:$limit:0.5:2" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   none:diff:$none:none:0 | \
   resilient:diff:$green:resilient:2 | \
   quick:diff:$red:quick:3 | \
   adaptive:diff:$blue:adaptive:1" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vKnone:diff:$none:none:0 | \
   vKe1.1:diff:$blue:1.1/1.3/1.6:1 | \
   vKe1.3:diff:$blue::-1 | \
   vKe1.6:diff:$blue::-1 | \
   vKe1.7:diff:$limit:1.7:2" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vKnone:diff:$none:none:0 | \
   vKm0.1:diff:$blue:0.1/0.15/0.4:1 | \
   vKm0.15:diff:$blue::-1 | \
   vKm0.4:diff:$blue::-1 | \
   vKm0.5:diff:$limit:0.5:2" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vKnone:diff:$none:none:0 | \
   vKresilient:diff:$green:resilient:2 | \
   vKquick:diff:$red:quick:3 | \
   vKadaptive:diff:$blue:adaptive:1" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vnone:diff:$none:none:0 | \
   ve1.1:diff:$blue:1.1/1.2/1.6:1 | \
   ve1.2:diff:$blue::-1 | \
   ve1.6:diff:$blue::-1 | \
   ve1.7:diff:$limit:1.7:2" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vnone:diff:$none:none:0 | \
   vm0.05:diff:$blue:0.05/0.1/0.3:1 | \
   vm0.1:diff:$blue::-1 | \
   vm0.3:diff:$blue::-1 | \
   vm0.4:diff:$limit:0.4:2" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vnone:diff:$none:none:0 | \
   vresilient:diff:$green:resilient:2 | \
   vquick:diff:$red:quick:3 | \
   vadaptive:diff:$blue:adaptive:1" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   VKnone:diff:$none:none:0 | \
   VKe1.2:diff:$blue:1.2/1.4/1.9:1 | \
   VKe1.4:diff:$blue::-1 | \
   VKe1.9:diff:$blue::-1 | \
   VKe2.0:diff:$limit:2.0:2" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   VKnone:diff:$none:none:0 | \
   VKm0.2:diff:$blue:0.2/0.3/0.5:1 | \
   VKm0.3:diff:$blue::-1 | \
   VKm0.5:diff:$blue::-1 | \
   VKm0.6:diff:$limit:0.6:2" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   VKnone:diff:$none:none:0 | \
   VKresilient:diff:$green:resilient:2 | \
   VKquick:diff:$red:quick:3 | \
   VKadaptive:diff:$blue:adaptive:1" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   Vnone:diff:$none:none:0 | \
   Ve1.2:diff:$blue:1.2/1.4/1.9:1 | \
   Ve1.4:diff:$blue::-1 | \
   Ve1.9:diff:$blue::-1 | \
   Ve2.0:diff:$limit:2.0:2" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   Vnone:diff:$none:none:0 | \
   Vm0.2:diff:$blue:0.2/0.3/0.8:1 | \
   Vm0.3:diff:$blue::-1 | \
   Vm0.8:diff:$blue::-1 | \
   Vm0.9:diff:$limit:0.9:2" \
   iris.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   Vnone:diff:$none:none:0 | \
   Vresilient:diff:$green:resilient:2 | \
   Vquick:diff:$red:quick:3 | \
   Vadaptive:diff:$blue:adaptive:1" \
   iris.res "$out"
let i++

if [[ $file == "" ]]; then term; fi

echo "done."
fi

#-----------------------------------------------------------------------
# wsel 3
#-----------------------------------------------------------------------

if [[ $data == "all" || $data == "wsel3" ]]; then
if [[ $file == "" ]]; then start wsel3.mp; fi
echo -n "processing wsel3 ... "

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   none00:diff:$trial:trials:0 | \
   none01:diff:$trial::-1 | \
   none02:diff:$trial::-1 | \
   none03:diff:$trial::-1 | \
   none04:diff:$trial::-1 | \
   none05:diff:$trial::-1 | \
   none06:diff:$trial::-1 | \
   none07:diff:$trial::-1 | \
   none08:diff:$trial::-1 | \
   none09:diff:$trial::-1 | \
   none10:diff:$trial::-1 | \
   none11:diff:$trial::-1 | \
   none12:diff:$trial::-1 | \
   none13:diff:$trial::-1 | \
   none14:diff:$trial::-1 | \
   none15:diff:$trial::-1 | \
   none16:diff:$trial::-1 | \
   none17:diff:$trial::-1 | \
   none18:diff:$trial::-1 | \
   none19:diff:$trial::-1 | \
   none:diff:$mean:average:1" \
   wsel3.trl "$out"
let i++

mpdiag $i \
  "log(difference):botlft | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vnone00:diff:$trial:trials:0 | \
   vnone01:diff:$trial::-1 | \
   vnone02:diff:$trial::-1 | \
   vnone03:diff:$trial::-1 | \
   vnone04:diff:$trial::-1 | \
   vnone05:diff:$trial::-1 | \
   vnone06:diff:$trial::-1 | \
   vnone07:diff:$trial::-1 | \
   vnone08:diff:$trial::-1 | \
   vnone09:diff:$trial::-1 | \
   vnone10:diff:$trial::-1 | \
   vnone11:diff:$trial::-1 | \
   vnone12:diff:$trial::-1 | \
   vnone13:diff:$trial::-1 | \
   vnone14:diff:$trial::-1 | \
   vnone15:diff:$trial::-1 | \
   vnone16:diff:$trial::-1 | \
   vnone17:diff:$trial::-1 | \
   vnone18:diff:$trial::-1 | \
   vnone19:diff:$trial::-1 | \
   vnone:diff:$mean:average:1" \
   wsel3.trl "$out"
let i++

mpdiag $i \
  "log(difference):botlft | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   Vnone00:diff:$trial:trials:0 | \
   Vnone01:diff:$trial::-1 | \
   Vnone02:diff:$trial::-1 | \
   Vnone03:diff:$trial::-1 | \
   Vnone04:diff:$trial::-1 | \
   Vnone05:diff:$trial::-1 | \
   Vnone06:diff:$trial::-1 | \
   Vnone07:diff:$trial::-1 | \
   Vnone08:diff:$trial::-1 | \
   Vnone09:diff:$trial::-1 | \
   Vnone10:diff:$trial::-1 | \
   Vnone11:diff:$trial::-1 | \
   Vnone12:diff:$trial::-1 | \
   Vnone13:diff:$trial::-1 | \
   Vnone14:diff:$trial::-1 | \
   Vnone15:diff:$trial::-1 | \
   Vnone16:diff:$trial::-1 | \
   Vnone17:diff:$trial::-1 | \
   Vnone18:diff:$trial::-1 | \
   Vnone19:diff:$trial::-1 | \
   Vnone:diff:$mean:average:1" \
   wsel3.trl "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   none:diff:$none:none:0 | \
   e1.2:diff:$blue:1.2/1.4/1.6:1 | \
   e1.4:diff:$blue::-1 | \
   e1.6:diff:$blue::-1 | \
   e1.7:diff:$limit:1.7:2" \
   wsel3.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   none:diff:$none:none:0 | \
   m0.1:diff:$blue:0.1/0.15/0.4:1 | \
   m0.15:diff:$blue::-1 | \
   m0.4:diff:$blue::-1 | \
   m0.5:diff:$limit:0.5:2" \
   wsel3.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   none:diff:$none:none:0 | \
   resilient:diff:$green:resilient:2 | \
   quick:diff:$red:quick:3 | \
   adaptive:diff:$blue:adaptive:1" \
   wsel3.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vKnone:diff:$none:none:0 | \
   vKe1.1:diff:$blue:1.1/1.3/1.6:1 | \
   vKe1.3:diff:$blue::-1 | \
   vKe1.6:diff:$blue::-1 | \
   vKe1.7:diff:$limit:1.7:2" \
   wsel3.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vKnone:diff:$none:none:0 | \
   vKm0.1:diff:$blue:0.1/0.2/0.4:1 | \
   vKm0.2:diff:$blue::-1 | \
   vKm0.4:diff:$blue::-1 | \
   vKm0.5:diff:$limit:0.5:2" \
   wsel3.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vKnone:diff:$none:none:0 | \
   vKresilient:diff:$green:resilient:2 | \
   vKquick:diff:$red:quick:3 | \
   vKadaptive:diff:$blue:adaptive:1" \
   wsel3.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vnone:diff:$none:none:0 | \
   ve1.1:diff:$blue:1.1/1.3/1.6:1 | \
   ve1.3:diff:$blue::-1 | \
   ve1.6:diff:$blue::-1 | \
   ve1.7:diff:$limit:1.7:2" \
   wsel3.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vnone:diff:$none:none:0 | \
   vm0.05:diff:$blue:0.05/0.1/0.2:1 | \
   vm0.1:diff:$blue::-1 | \
   vm0.2:diff:$blue::-1 | \
   vm0.3:diff:$limit:0.3:2" \
   wsel3.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   vnone:diff:$none:none:0 | \
   vresilient:diff:$green:resilient:2 | \
   vquick:diff:$red:quick:3 | \
   vadaptive:diff:$blue:adaptive:1" \
   wsel3.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   VKnone:diff:$none:none:0 | \
   VKe1.2:diff:$blue:1.2/1.4/1.8:1 | \
   VKe1.4:diff:$blue::-1 | \
   VKe1.8:diff:$blue::-1 | \
   VKe1.9:diff:$limit:1.9:2" \
   wsel3.res "$out"
let i++

mpdiag $i \
  "log(difference) | \
   1:0:0:*:*:-9.999:-9.999:*:* | \
   VKnone:diff:$none:none:0 | \
   VKm0.1:diff:$blue:0.1/0.2/0.6:1 | \