# HG changeset patch # User davidvanzessen # Date 1410781252 14400 # Node ID 8313c6cc65c50e4b0566e620c9f26c09ed00b901 # Parent 9641f3dfc590de6f5a17b32da98935508b598097 Uploaded diff -r 9641f3dfc590 -r 8313c6cc65c5 RScript.r --- a/RScript.r Mon Sep 15 05:38:46 2014 -0400 +++ b/RScript.r Mon Sep 15 07:40:52 2014 -0400 @@ -36,7 +36,7 @@ patients = split(dat, dat$Patient, drop=T) rm(dat) intervalReads = rev(c(0,10,25,50,100,1000,10000)) -intervalFreq = rev(c(0,0.01,0.1,0.5,1,5)) +intervalFreq = rev(c(0,0.01,0.05,0.1,0.5,1,5)) V_Segments = c(".*", "IGHV", "IGHD", "IGKV", "IGKV", "IgKINTR", "TRGV", "TRDV", "TRDD" , "TRBV") J_Segments = c(".*", ".*", ".*", "IGKJ", "KDE", ".*", ".*", ".*", ".*", ".*") Titles = c("Total", "IGH-Vh-Jh", "IGH-Dh-Jh", "Vk-Jk", "Vk-Kde" , "Intron-Kde", "TCRG", "TCRD-Vd-Dd", "TCRD-Dd-Dd", "TCRB-Vb-Jb")