Mercurial > repos > davidvanzessen > prisca

comparison RScript.r @ 7:7ce82833977c draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Tue, 12 Dec 2017 04:53:20 -0500
parents a9d2ed661541
children eb2aa7cffca3
comparison
equal deleted inserted replaced
6:a9d2ed661541 7:7ce82833977c
32 if(!("Total_Read_Count" %in% names(dat))){ 32 if(!("Total_Read_Count" %in% names(dat))){
33 dat$Total_Read_Count = 0 33 dat$Total_Read_Count = 0
34 } 34 }
35 35
36 dat = dat[,c("Patient", "Receptor", "Sample", "Cell_Count", "Clone_Molecule_Count_From_Spikes", "Log10_Frequency", "Total_Read_Count", "J_Segment_Major_Gene", "V_Segment_Major_Gene", "CDR3_Sense_Sequence", "Clone_Sequence")] 36 dat = dat[,c("Patient", "Receptor", "Sample", "Cell_Count", "Clone_Molecule_Count_From_Spikes", "Log10_Frequency", "Total_Read_Count", "J_Segment_Major_Gene", "V_Segment_Major_Gene", "CDR3_Sense_Sequence", "Clone_Sequence")]
37 dat = dat[!(nchar(dat$Clone_Sequence) < 2),]
37 38
38 dat$dsPerM = 0 39 dat$dsPerM = 0
39 dat = dat[!is.na(dat$Patient),] 40 dat = dat[!is.na(dat$Patient),]
40 dat$Related_to_leukemia_clone = F 41 dat$Related_to_leukemia_clone = F
41 42