changeset 30:b50965edac24 draft

Uploaded
author davidvanzessen
date Thu, 09 Mar 2017 10:48:14 -0500 (2017-03-09)
parents 1f83e14f173b
children 3a76faa53c59
files report_clonality/RScript.r
diffstat 1 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- 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")])