# HG changeset patch
# User davidvanzessen
# Date 1446804737 18000
# Node ID cf1def1aaa74d2977e865a618aff0cbf9bc11f3d
# Parent 6e068ef43336ca934f962304776bb747857246f9
Uploaded
diff -r 6e068ef43336 -r cf1def1aaa74 complete.sh
--- a/complete.sh Fri Aug 07 08:42:58 2015 -0400
+++ b/complete.sh Fri Nov 06 05:12:17 2015 -0500
@@ -52,7 +52,8 @@
mergerInput+=($parsedFileName)
count=$((count+1))
done
-wait
+
+echo "
after loooop |
" >> $html
echo "after loop"
diff -r 6e068ef43336 -r cf1def1aaa74 igblast/test.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/igblast/test.txt Fri Nov 06 05:12:17 2015 -0500
@@ -0,0 +1,1 @@
+test1 test2 test3
diff -r 6e068ef43336 -r cf1def1aaa74 report_clonality/RScript.r
--- a/report_clonality/RScript.r Fri Aug 07 08:42:58 2015 -0400
+++ b/report_clonality/RScript.r Fri Nov 06 05:12:17 2015 -0500
@@ -563,10 +563,6 @@
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("MEAN P3V.nt.nb:")
- print(PRODF$P3V.nt.nb)
- print(mean(PRODF$P3V.nt.nb, na.rm=T))
- print(head(PRODF))
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),
@@ -591,6 +587,8 @@
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)
newData = data.frame(data.table(UNPROD)[,list(unique=.N,
@@ -617,6 +615,7 @@
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)
}