Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

changeset 61:77a090cd0e02 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Fri, 30 Oct 2015 04:39:33 -0400
parents 28c3695259c1
children b1ad6f515338
files RScript.r
diffstat 1 files changed, 6 insertions(+), 6 deletions(-) [+]
line wrap: on
line diff
--- a/RScript.r	Thu Oct 29 11:21:33 2015 -0400
+++ b/RScript.r	Fri Oct 30 04:39:33 2015 -0400
@@ -442,14 +442,14 @@
 interval = intervalFreq
 intervalOrder = data.frame("interval"=paste(">", interval, sep=""), "intervalOrder"=1:length(interval))
 product = data.frame("Titles"=rep(Titles, each=length(interval)), "interval"=rep(interval, times=10), "V_Segments"=rep(V_Segments, each=length(interval)), "J_Segments"=rep(J_Segments, each=length(interval)))
-lapply(patients[4], FUN=patientCountOnColumn, product = product, interval=interval, on="Frequency", appendtxt=T)
+lapply(patients, FUN=patientCountOnColumn, product = product, interval=interval, on="Frequency", appendtxt=T)
 
 cat("<tr><td>Starting Cell Count analysis</td></tr>", file=logfile, append=T)
 
 interval = intervalReads
 intervalOrder = data.frame("interval"=paste(">", interval, sep=""), "intervalOrder"=1:length(interval))
 product = data.frame("Titles"=rep(Titles, each=length(interval)), "interval"=rep(interval, times=10), "V_Segments"=rep(V_Segments, each=length(interval)), "J_Segments"=rep(J_Segments, each=length(interval)))
-lapply(patients[4], FUN=patientCountOnColumn, product = product, interval=interval, on="normalized_read_count")
+lapply(patients, FUN=patientCountOnColumn, product = product, interval=interval, on="normalized_read_count")
 
 if(nrow(single_patients) > 0){
 	scales = 10^(0:6) #(0:ceiling(log10(max(scatterplot_locus_data$normalized_read_count))))
@@ -971,12 +971,12 @@
   one = triplets[triplets$Sample == "14696_reg_BM",]
   two = triplets[triplets$Sample == "24536_reg_BM",]
   three = triplets[triplets$Sample == "24062_reg_BM",]
-  #tripletAnalysis(one, "14696_1_Trio", two, "14696_2_Trio", three, "14696_3_Trio", product=product, interval=interval, on="normalized_read_count", T)
+  tripletAnalysis(one, "14696_1_Trio", two, "14696_2_Trio", three, "14696_3_Trio", product=product, interval=interval, on="normalized_read_count", T)
   
   one = triplets[triplets$Sample == "16278_Left",]
   two = triplets[triplets$Sample == "26402_Left",]
   three = triplets[triplets$Sample == "26759_Left",]
-  #tripletAnalysis(one, "16278_Left_Trio", two, "26402_Left_Trio", three, "26759_Left_Trio", product=product, interval=interval, on="normalized_read_count", T)
+  tripletAnalysis(one, "16278_Left_Trio", two, "26402_Left_Trio", three, "26759_Left_Trio", product=product, interval=interval, on="normalized_read_count", T)
   
   one = triplets[triplets$Sample == "16278_Right",]
   two = triplets[triplets$Sample == "26402_Right",]
@@ -992,12 +992,12 @@
   one = triplets[triplets$Sample == "14696_reg_BM",]
   two = triplets[triplets$Sample == "24536_reg_BM",]
   three = triplets[triplets$Sample == "24062_reg_BM",]
-  #tripletAnalysis(one, "14696_1_Trio", two, "14696_2_Trio", three, "14696_3_Trio", product=product, interval=interval, on="Frequency", F)
+  tripletAnalysis(one, "14696_1_Trio", two, "14696_2_Trio", three, "14696_3_Trio", product=product, interval=interval, on="Frequency", F)
   
   one = triplets[triplets$Sample == "16278_Left",]
   two = triplets[triplets$Sample == "26402_Left",]
   three = triplets[triplets$Sample == "26759_Left",]
-  #tripletAnalysis(one, "16278_Left_Trio", two, "26402_Left_Trio", three, "26759_Left_Trio", product=product, interval=interval, on="Frequency", F)
+  tripletAnalysis(one, "16278_Left_Trio", two, "26402_Left_Trio", three, "26759_Left_Trio", product=product, interval=interval, on="Frequency", F)
   
   one = triplets[triplets$Sample == "16278_Right",]
   two = triplets[triplets$Sample == "26402_Right",]