Mercurial > repos > davidvanzessen > argalaxy_tools

comparison report_clonality/RScript.r @ 36:fd55088e48d2 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Fri, 22 Jan 2016 10:22:24 -0500
parents f2010de70741
children 2e0a7c35082e
comparison
equal deleted inserted replaced
35:d7580295e13b 36:fd55088e48d2
490 490
491 for(smpl in unique(PRODF$Sample)){ 491 for(smpl in unique(PRODF$Sample)){
492 PRODF.sample = PRODF[PRODF$Sample == smpl,] 492 PRODF.sample = PRODF[PRODF$Sample == smpl,]
493 493
494 fltr = PRODF.sample$Top.V.Gene == "" 494 fltr = PRODF.sample$Top.V.Gene == ""
495 if(sum(fltr) > 0){ 495 if(any(fltr, na.rm=T)){
496 PRODF.sample[fltr, "Top.V.Gene"] = "NA" 496 PRODF.sample[fltr, "Top.V.Gene"] = "NA"
497 } 497 }
498 498
499 fltr = PRODF.sample$Top.D.Gene == "" 499 fltr = PRODF.sample$Top.D.Gene == ""
500 if(sum(fltr) > 0){ 500 if(any(fltr, na.rm=T)){
501 PRODF.sample[fltr, "Top.D.Gene"] = "NA" 501 PRODF.sample[fltr, "Top.D.Gene"] = "NA"
502 } 502 }
503 503
504 fltr = PRODF.sample$Top.J.Gene == "" 504 fltr = PRODF.sample$Top.J.Gene == ""
505 if(sum(fltr) > 0){ 505 if(any(fltr, na.rm=T)){
506 PRODF.sample[fltr, "Top.J.Gene"] = "NA" 506 PRODF.sample[fltr, "Top.J.Gene"] = "NA"
507 } 507 }
508 508
509 v.d = table(PRODF.sample$Top.V.Gene, PRODF.sample$Top.D.Gene) 509 v.d = table(PRODF.sample$Top.V.Gene, PRODF.sample$Top.D.Gene)
510 v.j = table(PRODF.sample$Top.V.Gene, PRODF.sample$Top.J.Gene) 510 v.j = table(PRODF.sample$Top.V.Gene, PRODF.sample$Top.J.Gene)