Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

comparison RScript.r @ 22:b662fdc7eff4 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 23 Feb 2015 10:37:28 -0500
parents d938aef60589
children 6904186d13b9
comparison
equal deleted inserted replaced
21:d980e2493657 22:b662fdc7eff4
49 49
50 dat = dat[dat$normalized_read_count >= min_cells,] 50 dat = dat[dat$normalized_read_count >= min_cells,]
51 51
52 dat$paste = paste(dat$Sample, dat$Clone_Sequence) 52 dat$paste = paste(dat$Sample, dat$Clone_Sequence)
53 53
54 patients = split(dat, dat$Patient, drop=T)[1:10] 54 patients = split(dat, dat$Patient, drop=T)
55 intervalReads = rev(c(0,10,25,50,100,250,500,750,1000,10000)) 55 intervalReads = rev(c(0,10,25,50,100,250,500,750,1000,10000))
56 intervalFreq = rev(c(0,0.01,0.05,0.1,0.5,1,5)) 56 intervalFreq = rev(c(0,0.01,0.05,0.1,0.5,1,5))
57 V_Segments = c(".*", "IGHV", "IGHD", "IGKV", "IGKV", "IgKINTR", "TRGV", "TRDV", "TRDD" , "TRBV") 57 V_Segments = c(".*", "IGHV", "IGHD", "IGKV", "IGKV", "IgKINTR", "TRGV", "TRDV", "TRDD" , "TRBV")
58 J_Segments = c(".*", ".*", ".*", "IGKJ", "KDE", ".*", ".*", ".*", ".*", ".*") 58 J_Segments = c(".*", ".*", ".*", "IGKJ", "KDE", ".*", ".*", ".*", ".*", ".*")
59 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") 59 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")