changeset 17:ee1bda8c27c8 draft

Uploaded
author davidvanzessen
date Mon, 09 Nov 2015 08:28:31 -0500
parents cf1def1aaa74
children 431797cd74c8
files igblast/test.txt report_clonality/RScript.r report_clonality/r_wrapper.sh
diffstat 3 files changed, 12 insertions(+), 15 deletions(-) [+]
line wrap: on
line diff
--- a/igblast/test.txt	Fri Nov 06 05:12:17 2015 -0500
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,1 +0,0 @@
-test1	test2	test3
--- a/report_clonality/RScript.r	Fri Nov 06 05:12:17 2015 -0500
+++ b/report_clonality/RScript.r	Mon Nov 09 08:28:31 2015 -0500
@@ -563,6 +563,7 @@
 imgtcolumns = c("X3V.REGION.trimmed.nt.nb","P3V.nt.nb", "N1.REGION.nt.nb", "P5D.nt.nb", "X5D.REGION.trimmed.nt.nb", "X3D.REGION.trimmed.nt.nb", "P3D.nt.nb", "N2.REGION.nt.nb", "P5J.nt.nb", "X5J.REGION.trimmed.nt.nb", "X3V.REGION.trimmed.nt.nb", "X5D.REGION.trimmed.nt.nb", "X3D.REGION.trimmed.nt.nb", "X5J.REGION.trimmed.nt.nb", "N1.REGION.nt.nb", "N2.REGION.nt.nb", "P3V.nt.nb", "P5D.nt.nb", "P3D.nt.nb", "P5J.nt.nb")
 if(all(imgtcolumns %in% colnames(inputdata)))
 {
+  print("found IMGT columns, running junction analysis")
   newData = data.frame(data.table(PRODF)[,list(unique=.N, 
                                                VH.DEL=mean(.SD$X3V.REGION.trimmed.nt.nb, na.rm=T),
                                                P1=mean(.SD$P3V.nt.nb, na.rm=T),
@@ -587,7 +588,6 @@
                                                            mean(.SD$P3D.nt.nb, na.rm=T) +
                                                            mean(.SD$P5J.nt.nb, na.rm=T))),
                                          by=c("Sample")])
-  print(newData)
   newData[,sapply(newData, is.numeric)] = round(newData[,sapply(newData, is.numeric)],1)
   write.table(newData, "junctionAnalysisProd.csv" , sep=",",quote=F,na="-",row.names=F,col.names=F)
   
@@ -602,18 +602,16 @@
                                                 N2=mean(.SD$N2.REGION.nt.nb, na.rm=T),
                                                 P4=mean(.SD$P5J.nt.nb, na.rm=T),
                                                 DEL.JH=mean(.SD$X5J.REGION.trimmed.nt.nb, na.rm=T),
-                                                Total.Del=(	mean(.SD$X3V.REGION.trimmed.nt.nb, na.rm=T) + 
-                                                            mean(.SD$X5D.REGION.trimmed.nt.nb, na.rm=T) + 
-                                                            mean(.SD$X3D.REGION.trimmed.nt.nb, na.rm=T) +
-                                                            mean(.SD$X5J.REGION.trimmed.nt.nb, na.rm=T)),
-                                                
-                                                Total.N=(	mean(.SD$N1.REGION.nt.nb, na.rm=T) +
-                                                          mean(.SD$N2.REGION.nt.nb, na.rm=T)),
-                                                
-																								Total.P=(	mean(.SD$P3V.nt.nb, na.rm=T) +
-																													mean(.SD$P5D.nt.nb, na.rm=T) +
-																													mean(.SD$P3D.nt.nb, na.rm=T) +
-																													mean(.SD$P5J.nt.nb, na.rm=T))),
+                                                Total.Del=(mean(.SD$X3V.REGION.trimmed.nt.nb, na.rm=T) + 
+                                                           mean(.SD$X5D.REGION.trimmed.nt.nb, na.rm=T) + 
+                                                           mean(.SD$X3D.REGION.trimmed.nt.nb, na.rm=T) +
+                                                           mean(.SD$X5J.REGION.trimmed.nt.nb, na.rm=T)),
+                                                Total.N=(  mean(.SD$N1.REGION.nt.nb, na.rm=T) +
+                                                           mean(.SD$N2.REGION.nt.nb, na.rm=T)),
+                                                Total.P=(  mean(.SD$P3V.nt.nb, na.rm=T) +
+							   mean(.SD$P5D.nt.nb, na.rm=T) +
+							   mean(.SD$P3D.nt.nb, na.rm=T) +
+						           mean(.SD$P5J.nt.nb, na.rm=T))),
                                           by=c("Sample")])
   newData[,sapply(newData, is.numeric)] = round(newData[,sapply(newData, is.numeric)],1)
   write.table(newData, "junctionAnalysisUnProd.csv" , sep=",",quote=F,na="-",row.names=F,col.names=F)
--- a/report_clonality/r_wrapper.sh	Fri Nov 06 05:12:17 2015 -0500
+++ b/report_clonality/r_wrapper.sh	Mon Nov 09 08:28:31 2015 -0500
@@ -137,7 +137,7 @@
 	echo "</div></div>" >> $outputFile
 fi
 
-hasJunctionData="$(if head -n 1 $inputFile | grep -q '3V-REGION trimmed-nt nb'; then echo 'Yes'; else echo 'No'; fi)"
+hasJunctionData="$(if head -n 1 $inputFile | grep -qE '3V.REGION.trimmed.nt.nb'; then echo 'Yes'; else echo 'No'; fi)"
 
 if [[ "$hasJunctionData" == "Yes" ]] ; then
 	echo "<div class='tabbertab' title='Junction Analysis'>" >> $outputFile