Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

comparison RScript.r @ 38:f5b242a5337f draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Fri, 11 Sep 2015 08:31:59 -0400
parents 623bbe972363
children 9fdcf6bfe024
comparison
equal deleted inserted replaced
37:623bbe972363 38:f5b242a5337f
315 } 315 }
316 p = NULL 316 p = NULL
317 if(nrow(scatterplot_locus_data) != 0){ 317 if(nrow(scatterplot_locus_data) != 0){
318 if(on == "normalized_read_count"){ 318 if(on == "normalized_read_count"){
319 scales = 10^(0:6) #(0:ceiling(log10(max(scatterplot_locus_data$normalized_read_count)))) 319 scales = 10^(0:6) #(0:ceiling(log10(max(scatterplot_locus_data$normalized_read_count))))
320 p = ggplot(scatterplot_locus_data, aes(type, normalized_read_count)) + scale_y_log10(breaks=scales,labels=scales) 320 p = ggplot(scatterplot_locus_data, aes(type, normalized_read_count)) + scale_y_log10(breaks=scales,labels=scales) + expand_limits(y=10^6)
321 } else { 321 } else {
322 p = ggplot(scatterplot_locus_data, aes(type, Frequency)) + scale_y_continuous(limits = c(0, 100)) + expand_limits(y=c(0,100)) 322 p = ggplot(scatterplot_locus_data, aes(type, Frequency)) + scale_y_continuous(limits = c(0, 100)) + expand_limits(y=c(0,100))
323 } 323 }
324 p = p + geom_point(aes(colour=type), position="jitter") 324 p = p + geom_point(aes(colour=type), position="jitter")
325 p = p + xlab("In one or both samples") + ylab(onShort) + ggtitle(paste(patient1[1,patientIndex], patient1[1,sampleIndex], patient2[1,sampleIndex], onShort, product[iter, titleIndex])) 325 p = p + xlab("In one or both samples") + ylab(onShort) + ggtitle(paste(patient1[1,patientIndex], patient1[1,sampleIndex], patient2[1,sampleIndex], onShort, product[iter, titleIndex]))