User Tools

Site Tools


histo2plot.sh

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

histo2plot.sh.1573539205.txt.gz · Last modified: 2021/03/17 13:09 (external edit)