Mercurial > repos > insilico-bob > ngchm
comparison mda_heatmap_gen/heatmap.sh @ 0:73aa1f80c237 draft
Upload
| author | insilico-bob |
|---|---|
| date | Thu, 07 Apr 2016 10:54:23 -0400 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| -1:000000000000 | 0:73aa1f80c237 |
|---|---|
| 1 echo $1 $2 $3 $4 $5 $6 $7 $8 $9 ${10} ${11} ${12} ${13} ${14} ${15} | |
| 2 | |
| 3 #create temp directory for row and col order and dendro files. | |
| 4 tdir=${11}/$(date +%s) | |
| 5 echo $tdir | |
| 6 mkdir $tdir | |
| 7 | |
| 8 #run R to cluster matrix | |
| 9 output="$(R --slave --vanilla --file=${11}/CHM.R --args $3 $4 $5 $6 $7 $8 $9 $tdir/ROfile.txt $tdir/COfile.txt $tdir/RDfile.txt $tdir/CDfile.txt)" | |
| 10 rc=$?; if [[ $rc != 0 ]]; then echo $output; exit $rc; fi | |
| 11 | |
| 12 | |
| 13 #there are a variable number of triplicate parameters for classification bars | |
| 14 count=0 | |
| 15 classifcations='' | |
| 16 | |
| 17 for i in "$@"; do | |
| 18 if [ $count -gt 10 ] | |
| 19 then | |
| 20 classifications=$classifications' '$i | |
| 21 fi | |
| 22 | |
| 23 count=$((count+1)) | |
| 24 done | |
| 25 | |
| 26 #call java program to generate NGCHM viewer files. | |
| 27 java -jar ${11}/GalaxyMapGen.jar $1 $2 DataLayer1 $3 linear Row Column $4 $5 $6 $tdir/ROfile.txt $tdir/RDfile.txt $7 $8 $9 $tdir/COfile.txt $tdir/CDfile.txt ${10} $classifications | |
| 28 | |
| 29 #clean up tempdir | |
| 30 rm -rf $tdir |
