Mercurial > repos > davidvanzessen > baseline_bayesian_estimation
diff wrapper.sh @ 5:d8de51314d3f draft
Uploaded
author | davidvanzessen |
---|---|
date | Thu, 14 Aug 2014 07:17:26 -0400 |
parents | 38c038c17d0c |
children | 8a4183eed741 |
line wrap: on
line diff
--- a/wrapper.sh Wed Aug 13 09:09:56 2014 -0400 +++ b/wrapper.sh Thu Aug 14 07:17:26 2014 -0400 @@ -1,3 +1,4 @@ + #!/bin/bash dir="$(cd "$(dirname "$0")" && pwd)" @@ -48,7 +49,7 @@ filtered="$PWD/filtered_${id}.txt" cat $PWD/$id/*/1_* > $summaryfile cat $PWD/$id/*/2_* > $gappedfile - Rscript $dir/filter.r $summaryfile $gappedfile "$selection" $filtered + Rscript $dir/filter.r $summaryfile $gappedfile "$selection" $filtered 2>&1 final="$PWD/final_${id}.txt" cat $filtered | cut -f2,4,7 > $final @@ -64,9 +65,9 @@ echo "$workdir/${outID}.txt" -rows=`tail -n +2 $workdir/${outID}.txt | grep -n 'Group' | grep -Eoh '^[0-9]+' | tr '\n' ' '` +rows=`tail -n +2 $workdir/${outID}.txt | grep -v "All sequences combined" | grep -n 'Group' | grep -Eoh '^[0-9]+' | tr '\n' ' '` rows=($rows) -unset rows[${#rows[@]}-1] +#unset rows[${#rows[@]}-1] cd $dir Rscript --verbose $dir/comparePDFs.r $workdir/${outID}.RData $output ${rows[@]} 2>&1