# HG changeset patch
# User davidvanzessen
# Date 1489074494 18000
# Node ID b50965edac240ea09ede1d916ab9af83f09c6cfe
# Parent  1f83e14f173bcf767697da558df4b8cb7552eaf8
Uploaded

diff -r 1f83e14f173b -r b50965edac24 report_clonality/RScript.r
--- a/report_clonality/RScript.r	Wed Mar 01 10:54:30 2017 -0500
+++ b/report_clonality/RScript.r	Thu Mar 09 10:48:14 2017 -0500
@@ -827,11 +827,11 @@
   newData = data.frame(data.table(PRODF.no.D)[,list(unique=.N, 
                                                VH.DEL=num_median(.SD$X3V.REGION.trimmed.nt.nb, na.rm=T),
                                                P1=num_median(.SD$P3V.nt.nb, na.rm=T),
-                                               N1=median(.SD$N.REGION.nt.nb, na.rm=T),
+                                               N1=num_median(.SD$N.REGION.nt.nb, na.rm=T),
                                                P2=num_median(.SD$P5J.nt.nb, na.rm=T),
                                                DEL.JH=num_median(.SD$X5J.REGION.trimmed.nt.nb, na.rm=T),
 											   Total.Del=num_median(rowSums(.SD[,c("X3V.REGION.trimmed.nt.nb", "X5J.REGION.trimmed.nt.nb"), with=F], na.rm=T)),
-											   Total.N=median(.SD$N.REGION.nt.nb, na.rm=T),
+											   Total.N=num_median(.SD$N.REGION.nt.nb, na.rm=T),
 											   Total.P=num_median(rowSums(.SD[,c("P3V.nt.nb", "P5J.nt.nb"), with=F], na.rm=T)),
 											   Median.CDR3.l=as.double(as.numeric(median(.SD$CDR3.Length, na.rm=T)))),
                                          by=c("Sample")])
@@ -859,11 +859,11 @@
     newData = data.frame(data.table(UNPROD.no.D)[,list(unique=.N, 
                                                 VH.DEL=num_median(.SD$X3V.REGION.trimmed.nt.nb, na.rm=T),
                                                 P1=num_median(.SD$P3V.nt.nb, na.rm=T),
-                                                N1=median(.SD$N.REGION.nt.nb, na.rm=T),
+                                                N1=num_median(.SD$N.REGION.nt.nb, na.rm=T),
                                                 P2=num_median(.SD$P5J.nt.nb, na.rm=T),
                                                 DEL.JH=num_median(.SD$X5J.REGION.trimmed.nt.nb, na.rm=T),
                                                 Total.Del=num_median(rowSums(.SD[,c("X3V.REGION.trimmed.nt.nb", "X5J.REGION.trimmed.nt.nb"), with=F], na.rm=T)),
-                                                Total.N=median(.SD$N.REGION.nt.nb, na.rm=T),
+                                                Total.N=num_median(.SD$N.REGION.nt.nb, na.rm=T),
                                                 Total.P=num_median(rowSums(.SD[,c("P3V.nt.nb", "P5J.nt.nb"), with=F], na.rm=T)),
                                                 Median.CDR3.l=as.double(as.numeric(median(.SD$CDR3.Length, na.rm=T)))),
 															by=c("Sample")])