Genome Informatics Laboratory at KRIBB

This is an old revision of the document!

---

#!/bin/sh
# by Haeyoung Jeong

FILENAME=$(basename "$1") #sample.hist or sample.pe.histo
#BASE=${FILENAME/.*}
BASE=${FILENAME/.hist*}
echo "File base: $BASE"

echo "set term png" > ${BASE}.gp
echo "set output \"${BASE}.png\"" >> ${BASE}.gp
#echo "set datafile separator \",\"" >> ${BASE}.gp
echo "set logscale x" >> ${BASE}.gp
echo "set logscale y" >> ${BASE}.gp
echo "set grid" >> ${BASE}.gp
echo "set xlabel \"kmer frequency\"" >> ${BASE}.gp
echo "set ylabel \"number of distinct kmer\"" >> ${BASE}.gp
echo "plot \"$1\" using 1:2 with lines" >> ${BASE}.gp

gnuplot ${BASE}.gp