Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples
diff RScript.r @ 66:ef3ac4f431bb draft
Uploaded
author | davidvanzessen |
---|---|
date | Fri, 20 Nov 2015 08:27:13 -0500 |
parents | 7a63e90cfc3d |
children | 40c72b9ffc79 |
line wrap: on
line diff
--- a/RScript.r Mon Nov 09 08:19:53 2015 -0500 +++ b/RScript.r Fri Nov 20 08:27:13 2015 -0500 @@ -17,7 +17,7 @@ #require(xtable) cat("<tr><td>Reading input</td></tr>", file=logfile, append=T) dat = read.table(inFile, header=T, sep="\t", dec=".", fill=T, stringsAsFactors=F) -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", "Related_to_leukemia_clone", "Clone_Sequence")] +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")] dat$dsPerM = 0 dat = dat[!is.na(dat$Patient),] dat$Related_to_leukemia_clone = F