# HG changeset patch
# User davidvanzessen
# Date 1386586013 18000
# Node ID e6209052aaf4b078581fbf15b9ef3b9df3d6b95d
# Parent 3da85f1a7f420b483de298ad9863774fb0fe7973
Uploaded
diff -r 3da85f1a7f42 -r e6209052aaf4 r_wrapper.sh
--- a/r_wrapper.sh Mon Dec 09 05:42:56 2013 -0500
+++ b/r_wrapper.sh Mon Dec 09 05:46:53 2013 -0500
@@ -17,23 +17,24 @@
samples=`cat $outputDir/samples.txt`
count=1
-echo "
$clonalType
" >> $outputFile
+echo "$clonalType
" >> $outputFile
+hasReplicateColumn="$(if head -n 1 $inputFile | grep -q 'Replicate'; then echo 'Yes'; else echo 'No'; fi)"
for sample in $samples; do
clonalityScore="$(cat $outputDir/ClonalityScore_$sample.csv)"
echo " |
" >> $outputFile
echo "$sample |
" >> $outputFile
-
+
+ echo "$hasReplicateColumn"
#if its a 'new' merged file with replicate info
- if head -n 1 $inputFile | grep -q "Replicate"
- then
- echo "Clonality Score: $clonalityScore |
" >> $outputFile
+ if [[ "$hasReplicateColumn" == "Yes" ]] ; then
+ echo "Clonality Score: $clonalityScore |
" >> $outputFile
#replicate,reads,squared
echo "Replicate ID | Number of Reads | Reads Squared | " >> $outputFile
while IFS=, read replicate reads squared
do
- echo "$replicate | $reads | $squared | " >> $outputFile
+ echo "$replicate | $reads | $squared | " >> $outputFile
done < $outputDir/ReplicateReads_$sample.csv
#sum of reads and reads squared
|