# HG changeset patch # User davidvanzessen # Date 1410773926 14400 # Node ID 9641f3dfc590de6f5a17b32da98935508b598097 # Parent f11df36f43bb1537d5d618570c40e3f9ec927168 Uploaded diff -r f11df36f43bb -r 9641f3dfc590 RScript.r --- a/RScript.r Mon Sep 15 05:37:16 2014 -0400 +++ b/RScript.r Mon Sep 15 05:38:46 2014 -0400 @@ -35,7 +35,6 @@ dat = dat[!duplicated(dat$paste),] patients = split(dat, dat$Patient, drop=T) rm(dat) -patients = patients[1:5] intervalReads = rev(c(0,10,25,50,100,1000,10000)) intervalFreq = rev(c(0,0.01,0.1,0.5,1,5)) V_Segments = c(".*", "IGHV", "IGHD", "IGKV", "IGKV", "IgKINTR", "TRGV", "TRDV", "TRDD" , "TRBV")