Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

comparison RScript.r @ 6:8313c6cc65c5 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 15 Sep 2014 07:40:52 -0400
parents 9641f3dfc590
children 68c6c7624ffc
comparison
equal deleted inserted replaced
5:9641f3dfc590 6:8313c6cc65c5
34 cat("<tr><td>Removing duplicates</td></tr>", file=logfile, append=T) 34 cat("<tr><td>Removing duplicates</td></tr>", file=logfile, append=T)
35 dat = dat[!duplicated(dat$paste),] 35 dat = dat[!duplicated(dat$paste),]
36 patients = split(dat, dat$Patient, drop=T) 36 patients = split(dat, dat$Patient, drop=T)
37 rm(dat) 37 rm(dat)
38 intervalReads = rev(c(0,10,25,50,100,1000,10000)) 38 intervalReads = rev(c(0,10,25,50,100,1000,10000))
39 intervalFreq = rev(c(0,0.01,0.1,0.5,1,5)) 39 intervalFreq = rev(c(0,0.01,0.05,0.1,0.5,1,5))
40 V_Segments = c(".*", "IGHV", "IGHD", "IGKV", "IGKV", "IgKINTR", "TRGV", "TRDV", "TRDD" , "TRBV") 40 V_Segments = c(".*", "IGHV", "IGHD", "IGKV", "IGKV", "IgKINTR", "TRGV", "TRDV", "TRDD" , "TRBV")
41 J_Segments = c(".*", ".*", ".*", "IGKJ", "KDE", ".*", ".*", ".*", ".*", ".*") 41 J_Segments = c(".*", ".*", ".*", "IGKJ", "KDE", ".*", ".*", ".*", ".*", ".*")
42 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") 42 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")
43 Titles = factor(Titles, levels=Titles) 43 Titles = factor(Titles, levels=Titles)
44 TitlesOrder = data.frame("Title"=Titles, "TitlesOrder"=1:length(Titles)) 44 TitlesOrder = data.frame("Title"=Titles, "TitlesOrder"=1:length(Titles))