# HG changeset patch
# User davidvanzessen
# Date 1385999807 18000
# Node ID dea2d3353a42f10e936769b700d4fde9c914fe02
# Parent 33ebd6f0d25616ac3e57ca81cb3cad97879483cc
Uploaded
diff -r 33ebd6f0d256 -r dea2d3353a42 r_wrapper.sh
--- a/r_wrapper.sh Mon Dec 02 10:56:40 2013 -0500
+++ b/r_wrapper.sh Mon Dec 02 10:56:47 2013 -0500
@@ -17,10 +17,42 @@
samples=`cat $outputDir/samples.txt`
count=1
-echo "
$clonalType
" >> $outputFile
+echo "$clonalType
" >> $outputFile
for sample in $samples; do
- echo " | " >> $outputFile
+ clonalityScore="$(cat $outputDir/ClonalityScore_$sample.csv)"
+ echo "
|
" >> $outputFile
echo "$sample |
" >> $outputFile
+
+ #if its a 'new' merged file with replicate info
+ if head -n 1 $inputFile | grep -q "Replicate"
+ 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
+ done < $outputDir/ReplicateReads_$sample.csv
+
+ #sum of reads and reads squared
+ while IFS=, read readsSum squaredSum
+ do
+ echo "Sum | $readsSum | $squaredSum | " >> $outputFile
+ done < $outputDir/ReplicateSumReads_$sample.csv
+
+ echo "
|
" >> $outputFile
+
+ #overview
+ echo "Coincidence Type | Raw Coincidence Freq | Coincidence Weight | Coincidences, Weighted | " >> $outputFile
+ while IFS=, read type count weight weightedCount
+ do
+ echo "$type | $count | $weight | $weightedCount | " >> $outputFile
+ done < $outputDir/ClonalityOverView_$sample.csv
+ echo "
|
" >> $outputFile
+ fi
+
echo "V-D Heatmap: | V-J Heatmap: | D-J Heatmap: |
" >> $outputFile
mv "$outputDir/HeatmapVD_$sample.png" "$outputDir/VD_$sample.png"
echo " | " >> $outputFile