Hapinfo to MHL and MF.job

From ZhangLabWiki
Jump to navigation Jump to search
#!/bin/csh
#PBS -N test
#PBS -q pdafm
#PBS -l nodes=1:ppn=1
#PBS -l walltime=14:00:00
#PBS -o test.o
#PBS -e test.e
#PBS -V
#PBS -M shihcheng.guo@gmail.com
#PBS -m abe
#PBS -A k4zhang-group
cd /home/shg047/oasis/Estellar2016/hapinfo
perl ~/bin/hapinfo2mhl.pl ./ > Estellar2016.hapinfo2mhl.march28.txt
Rscript --vanilla ~/bin/matrixPlot.R -i "Estellar2016.hapinfo2mhl.march28.txt" -o "Estellar2016.hapinfo2mhl.march28"
perl ~/bin/hapinfo2mf.pl ./ > Estellar2016.hapinfo2mf.march28.txt
Rscript --vanilla ~/bin/matrixPlot.R -i "Estellar2016.hapinfo2mf.march28.txt" -o "Estellar2016.hapinfo2mf.march28"