Mercurial > repos > davidvanzessen > complete_immunerepertoire_igg

comparison 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
comparison
equal deleted inserted replaced
7:a9053212a462 8:043fd6613fd9
550 } 550 }
551 551
552 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") 552 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")
553 if(all(imgtcolumns %in% colnames(inputdata))) 553 if(all(imgtcolumns %in% colnames(inputdata)))
554 { 554 {
555 newData = data.frame(data.table(inputdata)[,list(unique=.N, 555 newData = data.frame(data.table(PRODF)[,list(unique=.N,
556 VH.DEL=mean(X3V.REGION.trimmed.nt.nb, na.rm=T), 556 VH.DEL=mean(X3V.REGION.trimmed.nt.nb, na.rm=T),
557 P1=mean(P3V.nt.nb, na.rm=T), 557 P1=mean(P3V.nt.nb, na.rm=T),
558 N1=mean(N1.REGION.nt.nb, na.rm=T), 558 N1=mean(N1.REGION.nt.nb, na.rm=T),
559 P2=mean(P5D.nt.nb, na.rm=T), 559 P2=mean(P5D.nt.nb, na.rm=T),
560 DEL.DH=mean(X5D.REGION.trimmed.nt.nb, na.rm=T), 560 DEL.DH=mean(X5D.REGION.trimmed.nt.nb, na.rm=T),