Mercurial > repos > davidvanzessen > argalaxy_tools

comparison report_clonality/RScript.r @ 29:1f83e14f173b draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Wed, 01 Mar 2017 10:54:30 -0500
parents 28fbbdfd7a87
children b50965edac24
comparison
equal deleted inserted replaced
28:798b62942b4b 29:1f83e14f173b
893 D.REGION.reading.frame$percentage = round(D.REGION.reading.frame$Freq / D.REGION.reading.frame$sample.sum * 100, 1) 893 D.REGION.reading.frame$percentage = round(D.REGION.reading.frame$Freq / D.REGION.reading.frame$sample.sum * 100, 1)
894 894
895 write.table(D.REGION.reading.frame, "DReadingFrame.txt" , sep="\t",quote=F,row.names=F,col.names=T) 895 write.table(D.REGION.reading.frame, "DReadingFrame.txt" , sep="\t",quote=F,row.names=F,col.names=T)
896 896
897 D.REGION.reading.frame = ggplot(D.REGION.reading.frame) 897 D.REGION.reading.frame = ggplot(D.REGION.reading.frame)
898 D.REGION.reading.frame = D.REGION.reading.frame + geom_bar(aes( x = D.REGION.reading.frame, y = percentage, fill=Sample), stat='identity', position='dodge' ) + ggtitle("D reading frame") + xlab("Frequency") + ylab("Frame") 898 D.REGION.reading.frame = D.REGION.reading.frame + geom_bar(aes( x = D.REGION.reading.frame, y = percentage, fill=Sample), stat='identity', position='dodge' ) + ggtitle("D reading frame") + xlab("Frame") + ylab("Frequency")
899 D.REGION.reading.frame = D.REGION.reading.frame + scale_fill_manual(values=sample.colors) 899 D.REGION.reading.frame = D.REGION.reading.frame + scale_fill_manual(values=sample.colors)
900 D.REGION.reading.frame = D.REGION.reading.frame + theme(panel.background = element_rect(fill = "white", colour="black"),text = element_text(size=15, colour="black"), axis.text.x = element_text(angle = 45, hjust = 1), panel.grid.major.y = element_line(colour = "black"), panel.grid.major.x = element_blank()) 900 D.REGION.reading.frame = D.REGION.reading.frame + theme(panel.background = element_rect(fill = "white", colour="black"),text = element_text(size=15, colour="black"), axis.text.x = element_text(angle = 45, hjust = 1), panel.grid.major.y = element_line(colour = "black"), panel.grid.major.x = element_blank())
901 901
902 png("DReadingFrame.png") 902 png("DReadingFrame.png")
903 D.REGION.reading.frame 903 D.REGION.reading.frame