# HG changeset patch
# User davidvanzessen
# Date 1482240972 18000
# Node ID 02efa5764a0a10abdf5a6a9820f82976123d3482
# Parent  15961ca8d9ce3e4e275677859f646b1147e1cb83
Uploaded

diff -r 15961ca8d9ce -r 02efa5764a0a report_clonality/RScript.r
--- a/report_clonality/RScript.r	Tue Dec 20 06:26:44 2016 -0500
+++ b/report_clonality/RScript.r	Tue Dec 20 08:36:12 2016 -0500
@@ -374,7 +374,7 @@
 CDR3Length = merge(CDR3Length, TotalPerSample, by="Sample")
 CDR3Length$Frequency = CDR3Length$Count * 100 / CDR3Length$total
 CDR3LengthPlot = ggplot(CDR3Length)
-CDR3LengthPlot = CDR3LengthPlot + geom_bar(aes( x = CDR3.Length, y = Frequency, fill = Sample), stat='identity', position='dodge' ) + theme(axis.text.x = element_text(angle = 90, hjust = 1)) + 
+CDR3LengthPlot = CDR3LengthPlot + geom_bar(aes( x = factor(reorder(CDR3.Length, as.numeric(CDR3.Length))), y = Frequency, fill = Sample), stat='identity', position='dodge' ) + theme(axis.text.x = element_text(angle = 90, hjust = 1)) + 
   ggtitle("Length distribution of CDR3") + 
   xlab("CDR3 Length") + 
   ylab("Percentage of sequences") +