User Tools

Site Tools


Sidebar

This is the sidebar. Without it, the main text is too wide!

myillu_04_plotkmer.sh

Script <<myIllu_plotKmer.sh>>

#!/bin/sh

FILENAME=$(basename "$1") #sample.hist or sample.pe.fq.hist
#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
myillu_04_plotkmer.sh.txt · Last modified: 2016/05/07 13:05 by hyjeong