Mercurial > repos > davidvanzessen > complete_immunerepertoire_igg

diff RScript.r @ 8:043fd6613fd9 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 12 Jan 2015 04:32:40 -0500
parents a9053212a462
children 7dbc9ebcefce
line wrap: on
line diff
--- a/RScript.r	Mon Jan 05 09:30:08 2015 -0500
+++ b/RScript.r	Mon Jan 12 04:32:40 2015 -0500
@@ -552,7 +552,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)))
 {
-  newData = data.frame(data.table(inputdata)[,list(unique=.N, 
+  newData = data.frame(data.table(PRODF)[,list(unique=.N, 
                                               VH.DEL=mean(X3V.REGION.trimmed.nt.nb, na.rm=T),
                                               P1=mean(P3V.nt.nb, na.rm=T),
                                               N1=mean(N1.REGION.nt.nb, na.rm=T),