Mercurial > repos > davidvanzessen > argalaxy_tools
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")])